
NASA TN D-7518 


on 


case file 



THE EFFECT OF WIND-TUNNEL 

WALL INTERFERENCE ON THE PERFORMANCE 

OF A FAN-IN-WING VTOL MODEL 


by Harry H. Heyson 

Langley Research Center 
Hampton , Va. 23665 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION • WASHINGTON, D. C. • FEBRUARY 1974 




1. Report No. 2. Government Accession No. 

NASA TN D-7518 

3. Recipient's Catalog No. 

4. Title and Subtitle 

THE EFFECT OF WIND-TUNNEL WALL INTERFERENCE ON 
THE PERFORMANCE OF A FAN-IN- WING VTOL MODEL 

5. Report Date 

February 1974 

6. Performing Organization Code 

7. Author(s) 

Harry H. Heyson 

8. Performing Organization Report No. 

L-9202 

10. Work Unit No. 

760-61-02-01 

9. Performing Organization Name and Address 

NASA Langley Research Center 
Hampton, Va. 23665 

1 1 . Contract or Grant No. 

13. Type of Report and Period Covered 

Technical Note 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 20546 

14. Sponsoring Agency Code 

15. Supplementary Notes 

16. Abstract 

A fan-in-wing model with a 1.07-m (42-in.) span was tested in seven different test sec- 
tions with cross-sectional areas ranging from 2.2 m2 to 265 m2 (24 ft^ to 2857 ft2). The data 


from the different test sections are compared both with and without correction for wall inter- 
ference. The results demonstrate that extreme care must be used in interpreting uncorrected 
VTOL data since the wall interference may be so large as to invalidate even trends in the data. 
The wall interference is particularly large at the tail, a result which is in agreement with 
recently published comparisons of flight and large-scale wind-tunnel data (NASA CR-2135) for 
a propeller-driven deflected- slipstream configuration. The data of the present investigation 
verify the wall-interference theory of NASA TR R-124 even under conditions of extreme inter- 
ference. A method given by Tyler and Williamson in AGARD CP-91-71 yields reasonable esti- 
mates for the onset of Rae T s minimum-speed limit. 


The present investigation shows that the rules for choosing model sizes to produce negli- 
gible wall effects, as given by Cook and Hickey in NASA SP-116, are considerably in error and 
permit the use of excessively large models. Even simple momentum theory appears to yield 
more nearly correct performance estimates in transition flight than uncorrected wind-tunnel 


data when the model span approaches one-half of the tunnel width. The "fan- induced” lift indi- 
cated by a number of previous studies in which the model was of similar relative span appears 
to be largely the result of wall interference that was not accounted for in reducing the data. 


17. Key Words (Suggested by Author(s) ) 

Wind tunnels 
Wall interference 
Fan in wing 
VTOL 


18. Distribution Statement 
Unclassified 

- Unlimited 

STAR Category 11 

19. Security Classif. (of this report) 

20. Security Classif. (of this page) 

21. No. of Pages 

22. Price* 

Unclassified 

Unclassified 

230 

$5.75 


For sale by the National Technical Information Service, Springfield, Virginia 22151 
























CONTENTS 


Page 

SUMMARY 1 

INTRODUCTION : 1 

SYMBOLS 3 

APPARATUS AND TESTS 6 

Model 6 

Model Instrumentation 7 

Wind Tunnels 8 

Procedure 13 

Precision of Measurement 14 

RESULTS AND DISCUSSION OF DATA FROM MODEL WITH FANS COVERED ... 15 

Uncorrected Data 15 

Considerations in Correcting Data 16 

Procedure in Correcting Data 18 

Corrections 20 

Corrected Data 21 

RESULTS AND DISCUSSION OF DATA FROM MODEL WITH FANS 

OPERATING 21 

Uncorrected Data 21 

Correcting the Data 33 

Interference in Uncorrected Data 36 

Corrected Data 37 

Interference in Corrected Data 37 

Maximum Model Size 37 

CONCLUSIONS 38 

APPENDIX A - JUSTIFICATION OF VORTEX-DENSITY CORRECTION FOR 

FANS AND JETS 40 

APPENDIX B - FORTRAN PROGRAM FOR CORRECTING DATA FROM A FAN- 

IN- WING MODEL TESTED IN FOUR DIFFERENT TEST SECTIONS 41 

REFERENCES 66 

TABLES 70 

FIGURES 72 

iii 



THE EFFECT OF WIND-TUNNEL WALL INTERFERENCE ON THE 
PERFORMANCE OF A FAN-IN-WING VTOL MODEL 


By Harry H. Hey son 
Langley Research Center 

SUMMARY 

A fan-in-wing model with a 1.07-m (42-in.) span was tested in seven different test 
sections with cross-sectional areas ranging from 2.2 m2 to 265 m2 (24 ft2 to 2857 ft2). 

The data from the different test sections are compared both with and without correction 
for wall interference. The results demonstrate that extreme care must be used in inter- 
preting uncorrected VTOL data since the wall interference may be so large as to invali- 
date even trends in the data. The wall interference is particularly large at the tail, a 
result which is in agreement with recently published comparisons of flight and large-scale 
wind-tunnel data (NASA CR-2135) for a propeller-driven deflected- slipstream configura- 
tion. The data of the present investigation verify the wall-interference theory of NASA 
TR R-124 even under conditions of extreme interference. A method given by Tyler and 
Williamson in AGARD CP-91-71 yields reasonable estimates for the onset of Rae's 
minimum- speed limit. 

The present investigation shows that the rules for choosing model sizes to produce 
negligible wall effects, as given by Cook and Hickey in NASA SP-116, are considerably in 
error and permit the use of excessively large models. Even simple momentum theory 
appears to yield more nearly correct performance estimates in transition flight than 
uncorrected wind-tunnel data when the model span approaches one-half of the tunnel width. 
The "fan-induced" lift indicated by a number of previous studies in which the model was of 
similar relative span appears to be largely the result of wall interference that was not 
accounted for in reducing the data. 


INTRODUCTION 

Despite considerable theoretical study (e.g., refs. 1 to 6) of wind-tunnel interference 
for VTOL and STOL aircraft, it is not a general practice to correct all such data for wall 
effects. This failure to correct is due in part to conflicting reports of the efficacy of such 
corrections (e.g., refs. 7 and 8); it is due in part to some confusion between the effects of 
corrections and of the minimum-speed limits proposed by Rae (ref. 9); and finally, it is 


due in part to the rather considerable effort required to program corrections for data 
reduction when the programing may be significantly different for different types of model. 

The magnitude of wall interference and the extent to which the data may be corrected 
for such interference become of paramount importance in the design of a new wind tunnel 
because the required test-section dimensions must be selected so that the data from the 
tunnel will be representative of the model operating in free air. Consequently, in connec- 
tion with the design of a new full-scale subsonic wind tunnel (refs. 10 to 12), a major 
experimental study of wall effects was undertaken. This program was a joint effort of 
the Langley and Ames Research Centers of NASA. The model chosen was a simplified 
fan-in-wing aircraft differing from the model of reference 13 only in the addition of a 
large tail and a slight increase in wing- section thickness ratio. This model was chosen 
because it is considerably more complex, from a wall-effects viewpoint, than the models 
which have heretofore been used in V/STOL wall-effects investigations (e.g. , refs. 8 and 
14 to 18), and because the general type of configuration was representative of the config- 
urations of reference 19 and therefore could provide an evaluation of the conclusions of 
reference 19. The model was tested, with and without smaller test-section inserts, in a 
2.13- by 3.05-m (7- by 10-ft) wind tunnel at Ames as well as briefly in the 12.2- by 
24.4-m (40- by 80-ft) wind tunnel at Ames. It was also tested with and without a test- 
section insert in the 9.14- by 18.3-m (30- by 60-ft) Langley full-scale wind tunnel. 

The immediate objectives of the test program were twofold: First, since the tests 

conducted by Rae (ref. 9), which defined the problem of minimum-speed limits, were all 
conducted using relatively large single rotors, it was desired to examine the differences 
in these limits which might be caused by distributing much of the lift into two discrete 
highly loaded fans. Second, it was desired to obtain some experimental indication of the 
magnitude and correctability of the wall interference engendered by a model of this type. 
The approach used was to correct all the data to the maximum extent possible, and then to 
examine the differences in the data (both corrected and uncorrected) from tests under con- 
trolled conditions in the various test sections. 

Examination of the data indicates that the first of the aforementioned objectives was 
only partially achieved. Some insight was obtained into the relative magnitude of the 
minimum-speed limits in different size test sections; however, the results are not ade- 
quate to distinguish any order of relative merit between the different cross-sectional 
shapes of the test sections. The second objective was met in a more satisfying manner. 
The data presented herein demonstrate the extent to which V/STOL data from different 
wind tunnels can be correlated, even in the face of extraordinarily large wall interference. 

This wall-interference study is of particular interest since it demonstrates that, for 
V/STOL flight conditions, the interference may be of such magnitude that even the trends 
of the data may be incorrect. An example of correlation of wind-tunnel and flight-test data 
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for an entirely different aircraft (ref. 20) is presented to demonstrate that this observa- 
tion is also true for configurations totally different from that of the present investigation. 
Comparisons are made between the present work and previously published theoretically 
and experimentally chosen limits for V/STOL wind-tunnel testing (refs. 6 and 7). 

In correcting wind-tunnel data of the nature of those presented herein, the biggest 
problem is the lack of uniformity over the model of the wall-induced interference (ref. 21). 
Some compensation must be made for the varying effective angles of attack and dynamic 
pressures over the different components of the model. Thus, in order to correct the data 
in a complete manner, it is necessary to have at least a rudimentary theoretical treat- 
ment of the performance of each component as affected by changes in velocity and angle 
of attack. In the present case, a simple momentum theory for the lifting fan in cross flow 
was used. This theoretical treatment, based largely on reference 22, is presented in its 
entirety in a separate paper (ref. 23). Throughout the present paper, the theoretical pre- 
dictions of reference 23 are compared with the measured model performance as obtained 
both with and without wall interference. 

Reference 13 has noted that a vortex- density correction is needed in applying the 
theory of references 2 and 3 to the correction of data obtained for fan-supported models. 

A justification of this vortex- density correction is presented in appendix A. A sample of 
the FORTRAN programs used in correcting the data obtained in the present investigation 
is presented as appendix B. 


SYMBOLS 

Because of the limited font of characters available in the automatic figure -plotting 
equipment, certain symbols may vary between the text and the figures. Where this varia- 
tion occurs, the symbol used in the figures is shown parenthetically at the beginning of the 
definition. 

A aspect ratio, b^/s^ 

Aj^ momentum area of VTOL elements 

7T 2 

Ajyj momentum area of wing, - b 

A>p cross-sectional area of wind-tunnel test section 

b span of wing 

C D drag coefficient, D/qS w or D j q c S ^ 
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lift coefficient, L/qS^ or L /q c s 


W 


L,j 


J a 


'N,T 


se 


c 

D 

d E 

D 

d e 

h 

L 

Nn 


q 

q c 

q F 

q T 

R 

S F 

Srp 

s w 


lift coefficient based on fan-area and fan-exit dynamic pressure, 

lift-curve slope of wing, dC-^jda, per rad 

tail normal-force coefficient, N^yqS^ or N T y / q c S T 
chord 

drag due to lift, total drag less drag at a = 0° 
sum of D and D se 

drag equivalent to shaft power, P s jv 

equivalent fan diameter, ^4Sjp/ir 

height of fan exit above test- section floor 
lift 

tail normal force 
shaft power 

(Q) dynamic pressure of test-section flow, - V 2 

dt 

corrected dynamic pressure at wing 
(Q-p) correc t ed dynamic pressure at fans 
corrected dynamic pressure at tail 
body radius 
fan area 
tail area 
wing area 


A p Vj 2 s F 
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Tg static thrust 

V test-section, or forward, velocity 

Vj fan-exit velocity in static thrust, ^Tg^pSp 

Wq vertical induced velocity in forward flight, positive upward 

a angle of attack, angle between relative wind axis and longitudinal axis of 

model, positive nose up, deg 

AD increment in fan external drag resulting from changes in a and V (see 

eq. (24)) 

Aip difference in Aa at wing and fans, (Aq!) f - (Aa)^, deg except rad in equa- 

tions (24) and (25) 

Ai-p difference in Aar at wing and tail, (Aar),p - (Aar)^, deg 

AL increment in fan lift resulting from changes in or and V (see eq. (25)) 

ALi so-called "fan-induced" lift, total lift less the independent lifts of the fans 

and the model with the fans covered 

Au longitudinal component of wall-induced interference velocity, positive 

rearward 

Aw vertical component of wall-induced interference, positive upward 

Aar change in angle of attack caused by wall interference, referred to wing angle 

of attack unless otherwise subscripted, positive nose up, deg except rad in 
equations (13), (15), and (16) 

S e elevator deflection angle, positive trailing- edge down, deg 

5 u Q interference factor for longitudinal interference due to drag 

6 U p interference factor for longitudinal interference due to lift 
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6 W D interference factor for vertical interference due to drag 

6,„ T interference factor for vertical interference due to lift 

W, ±j 

e downwash angle at tail, positive downward, deg 

p mass density of air 

X wake skew angle, angle measured from vertical axis of test section to center 

of wake, positive rearward, deg 

X e effective wake skew angle, deg 

Subscripts: 


c 

corrected 

F 

fans 

T 

tall 

u 

uncorrected 

W 

wing 


APPARATUS AND TESTS 
Model 

The model used in this investigation is shown in figure 1 and pertinent dimensions 
are further detailed in tables I and n. The model consisted of a symmetrical streamline 
body 2.13 m (84 in.) long with a maximum diameter of 0.2 m (8 in.). A symmetrical 
tapered wing with a 1.07-m (42-in.) span was mounted at the midpoint of the body. The 
airfoil section at the wing tip was NACA 16-015 and the section increased in thickness to 
NACA 16-017 at the centerline of the body; straight-line fairings were used between these 
two stations. 

Two commercially available 0.2-m (8-in.) tip- turbine-driven fans were mounted on 
centers spaced 0.56 m (22 in.) apart at the midchord position of the wing. The inlets to 
these fans were of the simple bellmouth type obtained by providing a reasonable radius at 
the intersection of the fan duct and the upper wing surface. 
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A slab tail with a 0.76-m (30-in.) span and a 0.32-m (12.5-in.) chord was mounted 
symmetrically so that its trailing edge was coincident with the rearmost end of the fuse- 
lage. This tail was installed during all tests for which the data are presented herein. 

The model was mounted on a pivot at the midpoint of the fuselage and 6.67 cm 
( 2 | ln ) below the centerline of th e model. A linear actuator, installed between the 

mounting strut and a point farther rearward on the model, provided remote control of 
angle of attack. 


Model Instrumentation 

The model was designed to be operated on the normal external mechanical balances 
of the wind tunnels; thus, it was not necessary to provide a sting balance for measurement 
of the overall forces and moments. The mechanical balances involved are all of the 
simple platform type and have relatively poor resolution of moments for model forces of 
the magnitude encountered during the tests. The expected balance accuracy, together with 
some anticipated difficulty in setting precisely the same powered- lift flight conditions, 
precluded the possibility of obtaining accurate measurements of the effect of the tail on 
the moments by comparison of tail-on and tail-off tests. Consequently, the tail and tail- 
cone were mounted on the body by means of a commercial 1.9-cm-diameter (3/4-in.) 
six-component strain-gage balance. The primary measurement desired was the tail nor- 
mal force, and the balance had a maximum load capability of 445 N (100 lb) for this com- 
ponent of force. 

Numerous pressure and temperature transducers were provided in the independent 
pneumatic systems powering the two fans. The only measurement pertinent to the final 
results was the rotational speeds of the fans, for which magnetic pickups were provided 
in the fan casings. Considerable difficulty was experienced with this system because of 
60-Hz pickup during the tests. The initial series of tests was actually conducted by set- 
ting the fan rotational speeds with a stroboscopic tachometer. For subsequent series of 
tests, a discriminator circuit was constructed to minimize the pickup problem, and mag- 
netic tachometers with higher output were used. This aspect of the testing is discussed 
more completely in a later section of this paper. 

Angle of attack was measured by an accelerometer-type transducer mounted within 
the model, except in the 12.2- by 24.4-m (40- by 80 — ft) tunnel, where a selsyn indicator 
mounted at the actuator strut was used. Differences in the data-acquisition systems of 
the other two tunnels required the use of different transducers in each tunnel. In the 
smaller two tunnels, the accuracy of the overall system was approximately the same. In 
the largest tunnel the overall accuracy was somewhat less. 
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Wind Tunnels 


2.13- by 3.05-m (7- by 10-ft) tunnel .- The smallest of the three tunnels used in 
this investigation was the Ames 7-foot by 10-foot Subsonic Wind Tunnel No. 2. This tun- 
nel is described on pages 1-32 and 1-33 of reference 24. The model was mounted in the 
tunnel on a single unfaired strut (fig. 2). The pivot point at the top of the strut was on the 
centerline of the tunnel; thus, the model aerodynamic center was slightly above the tunnel 
centerline. 

Air was supplied to the fans by means of two 5. 1-cm-diameter (2-in.) hoses which 
were dressed closely to the front and back of the strut by means of guide rings. Below 
the floor of the tunnel, and above the balance frame, some slack was provided in the air 
lines to provide for the motion which occurred as a result of changes in angle of attack. 

An elaborate trapeze connection was provided between the balance frame and the main air 
supply. Tests of this system under pressure, with the model hoses blocked, indicated no 
measurable effect on the loads as seen by the balance. 

Instrument leads were taped tightly to the sides of the strut and were connected to 
the data-acquisition system by means of a large hanging loop of wiring below the balance 
frame. The gap between the strut and the floor was closed to a minimum by specially 
trimmed sheet metal screwed to the floor of the tunnel. 

Tunnel airspeed was measured by means of a pitot-static tube mounted from the 
ceiling of the tunnel. Corrections for position error are discussed in a later section of 
this paper. The tube was mounted 0.254 m (10 in.) below the tunnel ceiling, and the 
static-pressure holes of the tube were 1.33 m (52.5 in.) ahead of the model pivot point. 

(At zero angle of attack this location is 0.267 m (10.5 in.) ahead of the nose of the model.) 
The dynamic pressure measured by this tube was passed through a pressure transducer 
and then to both the data-acquisition system and a digital indicator, which was used as a 
speed reference during the tests. 

Since this tunnel has continuous speed control, it had been hoped to maintain a close 
control over tunnel speed during each set of data points; however, this was not possible in 
practice. At high speed, blockage of the tunnel caused by the powerful variation of fan 
momentum drag with speed and angle of attack resulted in excessive time losses in 
attempting to set the tunnel speed precisely. At the lowest speeds, recirculation effects 
became so severe that the tunnel speed was found to lope; the pulsations in the tunnel flow 
were obvious even to the ear. Consequently, the tunnel speed was taken as the average of 
three readings, each of which in turn was averaged over a time of 1.25 seconds. 

Tuft boards were placed on the floor of the tunnel for visual observations of the flow 
when recirculation began (refs. 9 and 13). 
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Inserts in 2.13- by 3.05-m (7- by 10-ft) tunnel .- In order to simulate still smaller 
test-section sizes, the insert technique of references 9, 13, 18, and 25 was used. Two 
rectangular test sections were simulated by means of two vertical walls (one of plywood 
and the other of transparent plastic to permit observation of the tufts) between which two 
horizontal surfaces were suspended to simulate the floor and ceiling of the small test 
sections. The entire assembly in each case was generously braced to insure stability 
and dimensional constancy. 

The first of these simulated test sections had a width of 1.83 m (6 ft) and a height 
of 1.22 m (4 ft) providing a width-height ratio of 1.5. The second test section had a width 
of 2.24 m (88 in.) and a height of 1.12 m (44 in.) providing a width-height ratio of 2.0. A 
third test section was obtained by fitting the 2:1 insert internally with sheet metal ends 
which were rolled to a semicircular cross section; thus a flat-oval test section having a 
width-height ratio of 2.0 was provided. All these inserts were 3.66 m (12 ft) long. The 
model was centered longitudinally within each insert. Photographs of the model mounted 
in these test sections are given in figure 3. 

The cross-sectional areas of the 1.5:1 rectangular insert and the 2:1 flat-oval insert 
were essentially identical with each other at 2.23 m 2 (24 ft 2 ). The 2:1 rectangular test 
section, with a cross-sectional area of 2.50 m 2 (26.89 ft 2 ), was approximately 12 percent 
larger in cross-sectional area. The choice of these sizes was not accidental. The 
dimensions of the 2:1 flat-oval test section were specifically chosen to represent the 
wing-span to tunnel-width ratios used in several Ames full-scale tunnel tests of fan and 
fan-in-wing models (e.g., refs. 26 to 32). 

Speed measurement in the inserts was by means of the same pitot- static tube used 
in the basic wind tunnel. The longitudinal location of this pitot- static tube was constant, 
and in each insert the vertical location was adjusted so that the tube was 25.4 cm (10 in.) 
from the insert ceiling. 

Each insert was generously tufted for visual flow observations; however, the curved 
sheet metal walls of the flat-oval section severely limited the field of view. 

Compressed air to drive the fans was supplied from a large high-pressure storage 
tank. This air supply was adequate to drive the fans at nominal rotational speed of 10 000 
and 12 000 rpm. 

No evidence of any flow inclination was found in the data. Thus, the wind-tunnel 

* 

stream angle was zero irrespective of the presence or absence of the inserts. Under 
these conditions, it was possible to set the model angle of attack directly to the desired 
values throughout the tests. 

In any wind-tunnel wall-effects investigation the relative sizes of the test sections 
are of vital importance. A sketch illustrating the relative sizes is presented in figure 4. 
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12.2- by 24.4-m (40- by 80-ft) tunnel. - In order to obtain conditions essentially 
free of wall constraints, the model was tested briefly in the 12.2- by 24.4-m (40- by 
80-ft) Ames full-scale tunnel. This wind tunnel is described in reference 24. The 
external balance of this tunnel was not designed to measure loads as small as those which 
were produced by the present model. In order to gain some increase in precision, the 
model was mounted with its span vertical (fig. 5) so that the lift could be measured by the 
side-force scales, which have a greater sensitivity than the lift scales. 

The model was mounted on the same strut that was used in the tests conducted in 
the previously described tunnel; however, the mechanical arrangements did not allow the 
air hoses to be dressed closely to the strut. Instead, angled fittings were provided at the 
model and at the base of the strut. The required motion of the hoses with angle of attack 
was obtained by bending the supply hoses. As may be seen in figure 5, the resulting 
installation was substantially less clean than the installation in the smaller Ames tunnel. 

A different pitch actuator was installed, and angle of attack was measured as a func- 
tion of the actuator extension. The least division of the angle-of-attack indicator was 
0.25°, and this reading was manually inserted into the data-acquisition system. A strain- 
gage balance was inserted into the actuator linkage in order to measure pitching moments; 
however, these measurements were invalidated by the omission of a static tare accounting 
for the moments imposed on the model by bending the air-supply hoses. 

The discriminator circuit and large magnetic pickups were used in measuring the 
fan rotational speeds. This arrangement substantially reduced the amount of 60-Hz noise 
accepted by the counters; however, the static thrust measurements indicated that some 
spurious counts were still obtained. The counters were not connected directly to the data- 
acquisition system; the readings were manually inserted into the system. 

The air supply in this tunnel was not adequate for continuous operation of the fans at 
12 000 rpm. Therefore, the tests were conducted at 10 000 rpm and at the maximum 
available rotational speed, which tended to be on the order of 11 500 rpm. 

Three different systems of tunnel flow-velocity measurement were employed. 

There were substantial disagreements in the data measured by the three systems. The 
staff of the tunnel provided their best estimates of the true velocities, and these values 
were punched into the data cards at a later date. 

The tail-balance readings were recorded on a second data-acquisition system. This 
second system proved troublesome, with obviously mispunched cards being obtained even 
while recording zeros. It is believed that reasonably accurate readings were obtained 
during initial tests with the fans covered; however, the data obtained with the fans oper- 
ating were so different from the data obtained in all the other test sections that they were 
rejected. Approximately half of the powered phase of the testing was complete when this 
system failed completely and no further tail data were obtained. 
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Initial tests with the fans covered indicated a very large stream -angle correction. 
Inasmuch as this tunnel is not equipped with flow-survey apparatus, it was not possible to 
obtain direct measurements of the flow inclination. Subsequent tests, with the model 
removed and the air hoses taped tangentially to the top of the strut, indicated that a large 
lift tare was also present in the data. No complete sequence of tare tests were conducted 
to obtain the precise magnitude of the tare. In analyzing the data, the stream angle and 
the lift and drag tares were obtained by finding those values that yielded the same per- 
formance as in all the other test sections when the fans were covered. These values 
were assumed to be unaltered by fan operation. 

The stream angle obtained in the foregoing manner was significantly different from 
that presumed to exist during the conduct of the tests. Consequently, the maximum true 
angle of attack obtained in this tunnel was several degrees less than that obtained in the 
other tunnels. 

The tests in this tunnel were conducted under the direction of Kenneth W. Mort, of 
the NASA Ames Research Center. 

9.14- by 18.3-m (30- by 60-ft) tunnel . - The deficiencies inherent in the tests con- 
ducted in the Ames 12.2- by 24.4-m (40- by 80-ft) tunnel were such that the resulting 
data were too ambiguous to be accepted as defining the free-air characteristics of the 
model. Consequently, more complete tests were conducted in the 9.14- by 18.3-m (30- 
by 60-ft) Langley full-scale tunnel. This tunnel is described in reference 33. Some 
later information on the wind tunnel is presented in references 24 and 34. 

The ground board normally used in the Langley full-scale tunnel was in place during 
these tests. The upper surface of this ground board is approximately 0.61 m (2 ft) above 
the lower edge of the jet boundary and thus reduces the cross-sectional area of the test 
section to 141.8 m^ (1527 ft^). By comparison, the model is very small; its wing area is 
less than one-half of 1 percent of the test-section cross-sectional area. 

Because of the size of this tunnel, it was necessary to prepare a new mounting strut 
for the model. The new strut was designed so that the model was mounted vertically on 
the centerline of the active region of the tunnel (4.26 m (14 ft) above the ground board). 

As nearly as possible, the uppermost 1.07 m (3.5 ft) of the strut was identical with the 
strut used in the smaller tunnel. The end fitting on this strut, the hoses and their 
arrangement, and the angle-of-attack actuator were the same as those which were used in 
the smallest tunnel. A close-fitting fairing was installed around the strut starting 1.07 m 
(3.5 ft) below the model and continuing downward to meet the ground board. All hoses 
and electrical leads were dressed to the strut in, as closely as possible, the identical 
manner in which they were installed in the smaller tunnel. Photographs of the model 
installed in the tunnel are presented in figure 6. 
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The air-pressure lines were brought across the balance in a trapeze arrangement. 
Tests conducted under pressure with the hoses blocked at the model indicated no effect on 
the balance readings. The instrument leads were carried across the balance by means of 
a large hanging loop. 

Prior to mounting the model on the strut, the region occupied by the model was sur- 
veyed with a pitot-static-pitch-yaw head. The dynamic pressure measured by this survey 
instrument was used to calibrate the velocity at the model as a function of static depres- 
sion in the tunnel settling chamber; this static depression in turn was used to determine 
the tunnel velocity during the tests. The survey also disclosed the presence of a signifi- 
cant stream angle (approximately 0.7°) at the model location. The presence of this 
stream angle was confirmed later by the raw data from the symmetrical model when it 
was tested with the fans covered. The effects of this stream angle have been removed 
from all the data presented herein. 

The Langley full-scale tunnel does not have continuous speed control throughout the 
velocity range covered in these tests; instead, it has some 24 discrete power settings, or 
"points." A number of these points appropriate to the prior tests in the smaller tunnel 
were selected. The actual velocity presented herein was determined from the average of 
no fewer than 10 samplings, spaced 1 second apart, of the static pressure. 

In order to accommodate the different data-acquisition systems in this tunnel, it was 
necessary to use a different type of angle-of-attack transducer within the model. Again, 
the values presented result from the average of no fewer than 10 samplings of the trans- 
ducer output. 

Insert in 9.14- by 18.3-m (30- by 60-ft) tunnel .- It was desired to insure continuity 
of the test results between the tests conducted in the two wind-tunnel facilities. Conse- 
quently, a 2.13- by 3.05-m (7- by 10-ft) insert, 6.4 m (21 ft) long, was built up around 
the model in the Langley tunnel without disturbing the mounted model on the strut. The 
insert was fitted with a simple 15.2-cm-diameter (6-in.) semicircular sheet metal bell- 
mouth inlet to discourage separation of the flow at the inlet. 

The insert was constructed of 1.9-cm-thick (3/4-in.) plywood and was rigidly 
braced by angle iron to insure dimensional stability during the tests. It was supported by 
pipe columns and cable bracing so that the model pivot point was on the centerline of the 
insert, that is, in the same location as in the tests at the Ames Research Center, and so 
that it was centered longitudinally on the model. Photographs of this installation are pre- 
sented in figure 7. 

Within the insert, the gap at the floor of the tunnel was reduced to minimal size by 
means of closely trimmed sheet metal plates screwed to the floor. The fairing around the 
lower portion of the strut was sealed to the exterior of the insert. 
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The flow velocity within the insert was measured from the average of four sets of 
total- and static-pressure measurements. The probes for these measurements were 
located 45.7 cm (18 in.) behind the leading edge of the insert and 30.5 cm (12 in.) inward 
from the walls of the insert. Since no divergence was built into the insert, a small cor- 
rection (approximately 4 percent) was made to the velocity in order to account for the 
difference in boundary-layer displacement thickness between the probe and the model 
locations. 

It was not possible to survey the flow within the insert walls with the existing equip- 
ment in the Langley full-scale tunnel. Stream angle was determined by finding the angle 
which was required to reduce the lift of the symmetrical model with the fans covered to 
zero at an angle of attack of zero. In this regard, a root- mean- square average of such 
angles for all tunnel speeds was used. The resulting stream angle was approximately 
-0.2° and is accounted for in all data presented herein. 

Air to power the fans was provided by a permanent compressor in the tunnel. It 
would have been desirable to maintain the same rotational speeds as were used in the 
earlier Ames tests. Unfortunately, the compressor proved inadequate in capacity for 
continuous operation at 12 000 rpm. Consequently, the tests at Langley were conducted 
at lower rotational speeds, 8000 and 10 000 rpm, which overlapped those in the other test 
sections. 

All the tests in the Langley full-scale tunnel, as well as all the tests in the Ames 
7-foot by 10-foot Subsonic Wind-Tunnel No. 2, were conducted under the personal super- 
vision and direction of Frank A. Lazzeroni, of the U.S. Army Air Mobility R&D Labora- 
tory, Ames Directorate. 


Procedure 

The same test procedure was used in all the tunnels and test-section inserts. First, 
the fans were started and brought to the required rotational speed. Generally, static 
thrust was measured, usually throughout the same angle-of-attack range as in the subse- 
quent forward-flight tests. Then, the tunnel was started and brought to the desired veloc- 
ity. Data were recorded in the following angle-of-attack sequence: 0°, -10°, -5°, 0°, 5°, 
10°, and 16°. The tunnel speed was then altered to the next desired speed. Although the 
angle-of-attack sequence was constant, the progression of tunnel speeds was not constant. 
The sequence of speeds was often reversed so that the test commenced with the highest 
speed and ended with static thrust. Even more erratic velocity sequences were used in 
the Langley tunnel, where, because of a pole change in the motor-control system at 
approximately 48 knots, it was often more convenient to descend in velocity to that speed, 
drop to the smallest velocity, and then increase tunnel speed to obtain velocities up to 
48 knots. 
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Data recording procedures differed in the three tunnel facilities because of differ- 
ences in the data- acquisition systems. In the smallest tunnel, the data were obtained as 
three sets of time-averaged data (over a 1.25-second period) and were punched on cards 
for off-line reduction. In this tunnel, angle of attack was set manually and "dialed" into 
the data system manually. A similar set of two independent systems was used in the 
largest tunnel. In the Langley full-scale tunnel the data were obtained as at least 10 (and 
often 25) sets of samplings (with essentially no time averaging on each of the sets); the 
data were stored on magnetic tape for off-line processing. In this latter system, angle of 
attack was included as one of the directly recorded variables. 

In each case, essentially no data were available during the actual testing. While 
this "blindness" may be a disadvantage during tests of a specific configuration, it is an 
advantage during tests of the present type because it eliminates any tendency to tinker 
with the model in order to obtain a preconceived result. 

Precision of Measurement 


Detailed examination of the data, together with the known capabilities of the external 
balances, indicates that the forces should be accurate to within the values shown in the 
following table: 


Tunnel facility 

Force 

Lift 

Drag 

Tail normal 
force 

N 

lb 

N 

lb 

N 

lb 

2.13 by 3.05 m (7 by 10 ft) 

±8.9 

±2.0 

±2.2 

±0.5 

±4.4 

±1.0 

9.14 by 18.3 m (30 by 60 ft) 

±13.3 

±3.0 

±4.4 

±1.0 

±4.4 

±1.0 

12.2 by 24.4 m (40 by 80 ft) 

±22.2 

±5.0 

±22.2 

±5.0 

±4.4 

±1.0 


The values given for the 12.2- by 24.4-m (40- by 80-ft) tunnel include an allowance 
for the ambiguous nature of the stream angle and the tares. 


It will be observed that these accuracies vary percentagewise according to the over- 
all level of forces observed, and further that they will be reflected in the zeros for the 
data as well as in the data points themselves. As a proportionate point of reference for 
those figures in which the data have been nondimensionalized with respect to static thrust, 
it should be noted that the static thrust for the complete model ranges from about 196 N 
(44 lb) at a nominal speed of 8000 rpm to about 480 N (108 lb) at a nominal speed of 
12 000 rpm. 
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In the Ames full-scale tunnel the dynamic pressure is believed accurate within 
5 percent; angle of attack, to within 2.0°. In all the other test sections, dynamic pressure 
is believed accurate to within 1 percent and angle of attack to within 0. 1°. 

All data presented herein have been corrected for the effects of stream angle on the 
forces and the angle of attack, where such correction is appropriate. Where adequate 
rotational- speed data were recorded, the quantities Vj and T g used in nondimension- 
alizing much of the data have been corrected for the actual rotational speed. Forces, 
where presented directly rather than as coefficients, have been corrected to standard 
density from the density at which the data were obtained. 

Corrections for wall effects are discussed separately at appropriate points in the 
discussion of the results. 

RESULTS AND DISCUSSION OF DATA FROM 
MODEL WITH FANS COVERED 

Uncorrected Data 

The uncorrected data for the model operating with the fans covered are presented 
in terms of lift, drag, and tail normal-force coefficients as a function of angle of attack 
in figure 8. Because of the small loads and the coarse sensitivity of the external balances 
employed, only the data from the highest dynamic pressure run in each test section are 
presented. 

The strut used to mount the model during these tests was not faired, and a different 
length of this strut was exposed to the full dynamic pressure of the tunnel in each test 
section. No series of tare runs was made to determine the tare loads in the data; how- 
ever, as noted earlier, the mounting arrangements near the model were as identical as 
possible during most of the tests. Consequently, an amount of drag equal to the entire 
drag of the model with fans covered at zero angle of attack has been removed from the 
data in this figure and in all subsequent figures. The resulting values of drag and drag 
coefficients may be considered to be approximately those due to lift. 

The data for each coefficient, as obtained in the 9.14- by 18.3-m (30- by 60-ft) tun- 
nel (where boundary interference is negligible due to the extremely small size of the 
model compared with the test section), were subjected to least- squares analysis. The 
resulting expressions for a quartic fit to the data are displayed as a curve on each figure. 

It will be observed that even though the model was symmetrical, the data do not 
quite possess the expected symmetries and antisymmetries with angle of attack. This 
result is rational, for the rearmost portions of the model at positive angle of attack were 
immersed in a region of lowered dynamic pressure behind the mounting strut and were 
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free of this region when at negative angle of attack. It is clear that under such conditions 
the emphasis placed on maintaining the mounting conditions as identical as possible in 
most of the test sections was entirely justified and necessary. 

Considerations in Correcting Data 

The data of figure 8 contain several types of boundary -induced interference. First, 
there is solid blockage. This interference is easily evaluated to a sufficient degree of 
accuracy from the compilation of studies presented in section 6:10 of reference 35. Next, 
there are the boundary -induced effects due to the presence of the lifting model within the 
tunnel. In the present case this last-named effect was obtained using the method of ref- 
erences 2 and 3 as implemented by the FORTRAN programs given in reference 36. 

The use of references 2 and 3 presents two problems when the theory is applied to 
a model for which the lift may be zero. First, the momentum theory (refs. 2 and 22) 
used to obtain the wake skew angle appears to fail when the lift is negative. This diffi- 
culty is resolved by calculating the skew angle using the absolute value of the lift in the 
equations and subsequently choosing the proper quadrant for the wake according to 
whether the lift is positive or negative. The second problem is that the computer pro- 
grams of reference 36 are arranged in such a manner that they yield the correct inter- 
ference factors only when the wake skew angle is greater than -90° and less than or equal 
to 90° (that is, the wake cannot pass upward as it passes rearward). Some rules for 
treating the calculation by symmetries are presented in reference 5; however, in the 
present case, where the wake skew angles are only slightly greater than 90° (slightly 
upward) , it is more convenient merely to extrapolate from the values calculated for the 
first quadrant. 

The model was somewhat unusual in that the wing was closely coupled to an extra- 
ordinarily large tail. Furthermore, the tail had a greater aspect ratio than the wing (2.4 
compared with 1.6), and thus would be expected to have a higher lift-curve slope than the 
wing. Under such conditions, it would be expected that the model would behave more 
nearly as a tandem-wing system than as a simple wing-tail combination; this expectation 
is confirmed by the nonlinear character of the lift-curve slope (fig. 8(a)). It is important 
to consider this tandem-wing-like character of the model in the corrections; that is, the 
effect of the interference at the tail must be considered not only with respect to tail nor- 
mal force, but also with respect to the overall lift and drag of the model. 

The appropriate interference factors for the wing due to its own presence may be 
obtained from the FORTRAN program given as appendix B of reference 36. It was 
assumed that the wing had an elliptic load distribution. Since the quarter chord of the 
model is displaced from the pivot point, both vertically and longitudinally, these interfer- 
ence factors will be a function of angle of attack by virtue of the different vertical location 
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of the wing within the tunnel at each angle of attack. (See eq. (58) of ref. 3.) The effect 
of the presence of the tail and its loads on the tail itself is also significant and can be 
obtained from the same program; it is imperative that the location of the tail as a func- 
tion of angle of attack be considered. Observe that this effect would have been difficult 
to consider if it had not been for the use of a tail balance to measure the tail loads. The 
interference factors at the tail due to the presence of the wing may be obtained from 
appendix D of reference 36. The interference at the wing due to the presence of the tail 
could be obtained from the same program (by considering the wing to be a canard tail); 
however, the rapid decrease of interference with distance upstream from the causative 
lifting element precludes any significant effect from this source and it may be ignored 
safely. 

In the tests of the model installed in the inserts in the 2.13- by 3.05-m (7- by 10-ft) 
wind tunnel at the Ames Research Center, one other feature must be considered. This 
feature is the tunnel velocity measurement by means of a pitot- static tube near the nose 
of the model. At this location, the pitot- static tube is affected by the direct field of the 
model (both due to the body shape and to the lifting system) as well as by wall effects 
caused by the presence of the model. These effects must be evaluated in order to obtain 
the proper tunnel velocity to use in the interference calculations and in forming the cor- 
rected force coefficients. 

The solid blockage at the pitot locations is caused primarily by the body because the 
body is the portion of the model closest to the pitot tube. The blockage is not the same 
as a classical blockage correction (ref. 35), since the classical blockage calculation is 
for the model location. In the present analysis, this blockage effect was approximated by 
setting up a calculation based on the use of a source and a sink to represent a Rankine 
ovoid (ref. 37, p. 208), and using the technique of reference 38 to obtain the strengths and 
spacing of these elements to produce an ovoid which, in free-air, has the same length and 
diameter as the fuselage. The ovoid was then reflected both horizontally and vertically 
to produce a pattern which represents the boundary conditions at the walls. The interfer- 
ence velocities at the pitot-static tube location can be obtained from this field of elemental 
sources and sinks. 

It is the usual practice in such wall-effects calculations to omit the central image 
which represents the model itself on the basis that this is the portion to be measured and 
corrected. In the present calculations, the central image is retained since it is desired 
to include the direct field of the model as well as the blockage interference. The level of 
this correction is approximately 1 percent of the free-stream velocity. 

The foregoing treatment was used in the present analysis; however, it does contain 
certain inadequacies. First, the existence of the images representing the boundary condi- 
tions at the wall results in an overall velocity at the model which is somewhat greater than 
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the free-air condition for which the source and sink were chosen. Thus, the body for 
which the interference is obtained will be somewhat more slender than the desired body 
shape, and, furthermore, it will be slightly different in shape than a Rankine ovoid. These 
effects probably result in an underestimate of the actual interference. Second, the cal- 
culation method will produce a streamlined symmetrical body only at an angle of attack of 
zero. Thus, it is not possible to examine the effect of the angle of attack of the body on 
the calculated interference. Such effects would be expected to be large at positive angles 
of attack, where the nose approached more closely the pitot-static tube location; however, 
one would expect only smaller changes at negative angles of attack, where the model nose 
moved farther away from the tube. On an overall average basis, the actual effects of 
solid blockage at the pitot location probably are underestimated from the omission of 
angle -of- attack effects. 

The wall interference at the pitot location can be obtained from appendix D of ref- 
erence 36 by considering the pitot to be a canard tail of zero span. Considerable care 
must be exercised in choosing the tail length and height as a function of model angle of 
attack in order to retain the correct pitot location. The direct field of the lifting model 
is obtained by retaining the central image. This is accomplished most simply by altering 
the subroutine DLTAS given in appendix Q of reference 36. (Delete lines (Q13) and (Q67) 
through (Q105).) 


Procedure in Correcting Data 

The first step in correcting the data is to divide the loads between the wing and the 
tail. This is possible only because the present model was fitted with a tail balance. 

Then, the loads assigned to the wing are used to solve the momentum quartic (ref. 2 or 22) 
for v/wq and the wake skew angle x- 1° those test sections where the tunnel velocity 
was measured near the model, the measured tunnel velocity is then corrected, and the 
corrected value is used to recompute v/wq and \- 

The value of x obtained in this manner is the momentum-theory value and, as 
pointed out by reference 8, is not the value that should be used in wall-interference cal- 
culations. Because of wake roll-up, the wake vorticity will be concentrated at some 
higher location in the tunnel given by an effective average value of the skew angle x e - 
As discussed in references 6 and 39, the most appropriate choice for a simple wing is 
that given by 


tan *e 



( 1 ) 


The values of the interference factors (previously obtained from ref. 36) are then 
interpolated to obtain the values corresponding to this value of x e - In this range of skew 
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angles, the effect of x e on the interference factors for the effect on any element due to 
its own presence is small; thus, x may be assumed to be 90° when considering the 
effect of the tail on itself. 

At this point the lift and drag of the wing and tail may be used to compute the indi- 
vidual vertical and horizontal increments of interference velocity separately at the wing 
and at the tail (eqs. (40) to (43) of ref. 2). At the wing, the total values of Aw and Au 
are simply the sums of the respective components occasioned by the lift and drag of the 
wing; however, at the tail, the contributions of both the wing (which are different at the 
tail than at the wing) and the tail itself must be summed to obtain the total components of 
interference. 

The total values of Aw and Au are then used to obtain separately, at the wing 
and at the tail, the values of Aa and q c jq by use of equations (48b) and (49b) of refer- 
ence 2, which are 


A a = tan 


-1 Aw/V 

, Au 
1 + — 


V 


q 



( 2 ) 


( 3 ) 


The values of Aa and q c /q at the wing are used as a first correction to the data; 
however, it is also necessary to account for the differences in Aa and q c jq at the wing 
and the tail. These differences are conveniently expressed as 


Ai T = (Aa) T - (Aa) w (4) 



It will be noted that the difference in the two values of Aa is effectively a change 
in tail incidence (ref. 40), and q-j^q c is an alteration in the effective dynamic- pressure 
ratio at the tail. If these effects were not removed from the data, the model would not be 
aerodynamically equivalent to the model under test. 


19 



The procedure used herein was to resolve the forces of the entire model around a 
new effective stream direction given by A a, where 


(a) c = («) u + 



(c L ) u cos A “ - ( c d) u 


« c/l 


( c D ) u cos ia + (C L ' 
%/'> 


sin Aa? 


sin A a 


(6) 

( 7 ) 

( 8 ) 


Next the tail forces as measured were resolved, and then the tail lift and drag were 
adjusted for Ai t and q T /q c . This adjustment requires a knowledge of the lift-curve 
slope of the tail and the free-air dynamic-pressure ratio at the tail. It would be desirable 
to have test results for the tail in the presence of the body, but without the wing, as a guide 
in estimating these values. Unfortunately, the construction of the model did not allow for 
such tests; thus, the lift-curve slope was taken as 0.03 per degree (approximately the 
value given in fig. 5-5 of ref. 41) and the dynamic-pressure ratio at the tail was rather 
arbitrarily selected to be 0.9. A small correction to the induced drag of the tail was 
made to account for the difference in the measured and adjusted lift. A correction to the 
profile drag would be appropriate; however, insufficient data were available to make such 
an adjustment. In any event, such a profile-drag adjustment probably would be significant 
only if the tail were to stall during the test. 

The foregoing adjustments were sufficient to provide the corrected values of C N.T- 
As a final step, the differences in the lift and drag of the tail were applied as adjustment 
to the overall lift and drag of the model in order to obtain the final corrected values of 
C L and C D . 


Corrections 

The corrections obtained for the model with the fans covered are shown in figure 9. 
The dynamic pressure ratios differ from unity only by 2 or 3 percent and thus have only a 
comparatively small effect on the data. However, A a and Airp assume significant 
proportions in the smallest test sections at large angles of attack. 
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One common rule of thumb in wind-tunnel testing (e.g., ref. 42) is that Aa should 
not exceed 2°. It is evident from figure 9(a) that Aa has assumed almost this value in 
the smaller inserts even with the fans covered, and that Ai,p (fig. 9(c)) is well in excess 
of 2°, yielding a total correction angle at the tail on the order of 4° (eq. (4)). 

Corrected Data 

After the application of corrections, the data for the model with the fans covered 
appear as shown in figure 10. The solid line shown in figure 10 is again a least-squares 
quartic faired through the data obtained in the 9.14- by 18.3-m (30- by 60-ft) tunnel. 

On the basis of the force accuracies previously given and the dynamic pressures of 
the tests, the anticipated agreement should be on the order of 0.03 for C^, 0.008 for C^, 
and 0.03 for C N.T- Examination of figure 10 indicates that the correlation between C L 
and Cp is generally within these limits but that rp appears to be overcorrected to 
a somewhat greater extent than would be anticipated by a simple examination of the mea- 
surement accuracy at the highest values of lift. 

One possible cause of the poorer correlation in the case of the tail normal-force 
coefficient could lie in the required estimates of the tail lift- curve slope and tail dynamic- 
pressure ratio. These estimates are far more critical in correcting tail normal force 
than in correcting the overall lift and drag of the model. Another possible cause could be 
the effect of the wall-induced velocities in relocating the wake to a higher position in the 
small test sections than in the large tunnel. References 43 and 44 have examined this 
latter effect theoretically. The maximum ratios of Cj^A encountered in the present 
test are on the order of 0.5. For such values of references 43 and 44 indicate 

that the correction to the tail should increase when the tail moves with the model. Such 
an effect would further degrade the present correlation. Finally, the assessment of test 
accuracy may be excessively optimistic since the value quoted represents only 1 percent 
of the full normal-force capability of the tail balance and considerable vibration and buf- 
feting of the tail was obvious during the tests. 

RESULTS AND DISCUSSION OF DATA FROM 
MODEL WITH FANS OPERATING 

Uncorrected Data 

Presentation of data .- In view of the difficulties experienced with the measurement 
of fan rotational speed during the tests, it is not possible to present the data directly for 
constant rotational speeds. Instead, the forward velocity has been nondimensionalized 
with respect to Vj, which is the fan efflux velocity in static thrust, defined from simple 
incompressible momentum theory as 
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( 9 ) 



Similarly, forces are presented only in nondimensional quantities, generally referenced 
either to the static thrust or to each other. The static thrust used in these nondimension- 
alizations is always the value obtained in the largest test section used in each series of 
tests. It is also measured at zero angle of attack. These two conditions insure that the 
value of static thrust used is the best available from the viewpoint of minimum flow 
recirculation in the tunnel during the measurement. Indeed, comparisons of the static- 
thrust data, with and without the insert, in the Langley full-scale tunnel, indicate that any 
errors caused by recirculation in the 2.13- by 3.05-m (7- by 10- ft) tunnel are within the 
accuracy of the data. 

The uncorrected data in the form of lift, drag, and normal-force coefficients are 
presented in figures 11 to 13. They are presented in the form of the ratios of l/Tq, 
D/Tg, and D/L in figures 14 to 16. Finally, the ratio l/d e ^where Dg is the sum 
of the external drag and a drag equivalent to the power supplied to the fans (ref. 23)^ is 
presented in figure 17. 

For lift coefficients and the ratios of lift to static thrust, a line on the figures indi- 
cates the values which would be obtained if the lift were simply the direct sum of the ver- 
tical component of the fan static thrust and the lift of the wing with the fans covered. In a 
similar manner, the momentum-theory values of all the other parameters ^with the excep- 
tion of Cj^ rp, for which momentum theory is inappropriate^ have been computed by means 
of the equations of reference 23, and these calculated values are compared with the cor- 
rected data. 

When examining figure 17 it should be noted that no measurements adequate for the 
calculation of the power supplied directly to the fans were actually made. Indeed, in view 
of the small size and fairly low efficiency of the model fan turbines, such measurements 
of power would be meaningless in relation to flight hardware. Instead, the momentum- 
theory value of shaft power, as computed from reference 23, has been converted into an 
effective drag by means of the relationship 

P s 

D = — (10) 

se y v ' 

The values of D se obtained from equation (10) have been added to both the experi- 
mental data and the theoretical curve. Although figure 17 presents no measured data that 
were not available in the preceding figures, it does serve the purpose of illustrating the 
effects of wall interference on the efficiency of this type of aircraft in transition. 
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Effect of wall interference on lift .- Figure 11 indicates clearly that at any constant 
angle of attack, the measured lift coefficient increases as the cross-sectional area of the 
test section decreases. The magnitude of this effect is disguised somewhat by the log- 
arithmic scales and by the effect of the variation in dynamic pressure in computing the 
coefficient when a large part of the lift (from the fans) is essentially independent of for- 
ward speed. Figure 14 presents a truer picture of the influence of the walls by presenting 
the lift in the form of a ratio to the static thrust. Reference 23 has shown that the ratio 
l/t s is proportional to a lift coefficient based on fan area and fan- exit velocity; that is 


L 

T S 


= ic 


2 L,j 


( 11 ) 


In figure 14 the differences in the data, as measured in the various test sections, 
are demonstrated to represent very significant differences in lift. For example, at an 
angle of attack of zero (fig. 14(b)) and a speed of y/Vj = 0.4 (which would represent a 
speed near the high-speed end of transition), the data from the small inserts indicate that 
a lift of about 25 percent more than the static thrust would be obtained; the data from the 
moderately larger 2.13- by 3.05-m (7- by 10-ft) test section would indicate that the gain 
in lift would be only 10 percent; and the data from the largest test section indicate that a 
small loss in lift would be encountered. Indeed, the data from the largest test section 
indicate that this model would have a loss of lift (from that expected from a simple addi- 
tion of lift components) for all angles of attack and for all forward speeds less than 
y/Vj = 0.5. This speed range encompasses the entire feasible transition range of lifting 
fans with modern pressure ratios. 

Reference 7 presents a set of charts which define relative proportions between 
model and test section which were believed to yield negligible wall effects at a speed of 
30 knots. Figure 18 shows the degree to which the present tests meet these size limits. 
Only the highest and the lowest disk loadings encountered are shown. The test conditions 
include points between these two disk loadings. In particular, the present tests in the 
1.12- by 2.24-m (44- by 88-in.) insert meet these limits at least as well as many of the 
tests reported in references 26 to 32. In the 2.13- by 3.05-m (7- by 10-ft) tunnel, the 
size of the present model falls well within the size limitations of reference 7. In con- 
trast, the data in figure 14 clearly indicate unacceptably large overestimates of "fan- 
induced" lift in these test sections. Therefore, it must be concluded that the size limits 
proposed by reference 7 are not valid for configurations such as that of the present inves- 
tigation; indeed, since there is nothing very unusual about this configuration except the 
relative size of the tail, it would be presumed that these size limits are equally inapplica- 
ble to other configurations as well. 
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Reference 7 attempts to limit its conclusions to conditions for which the overall 
drag of the model is trimmed. This limitation to zero net drag is based upon the the- 
oretical results of reference 2, which the authors of reference 7 claim to be incorrect 
and even in the wrong direction. In fact, reference 2 was misinterpreted in arriving at 
the limitation to zero net drag. As will become clear in the subsequent discussion of 
correcting the present fan-in-wing data, each element of the aircraft must be considered 
individually. Thus, the fans, except for a few isolated conditions, always have a drag; 
also, the wing always has an induced drag. The addition of a centered jet exhausting 
directly rearward (y = 90°) could balance the drag of the model under any condition; how- 
ever, the thrusting jet would contribute nothing to the interference at the model (from 
ref. 2, 6 , T = b = 6„ r» = 0, and 6,„ T has no effect since the lift of the thrusting 

U j Jj W j U VX W j Xj 

jet would be zeroj. Consequently, the limitation to trimmed drag in reference 7 is 
meaningless. 

The comparisons between flight and wind-tunnel data given in reference 7 have 
already been discussed in reference 39, which shows that the conclusions of reference 7 
were based upon faulty comparisons between flight and wind tunnel. Such an error should 
have been anticipated, since one of the conclusions was that both Glauert's corrections 
and those of reference 2 were in the wrong direction. Since both Glauert and reference 2 
predict upwash interference in a closed tunnel, this result of reference 7 could only be 
obtained if a downwash interference was produced by the walls. Such a result is physi- 
cally impossible for an overall correction in a closed tunnel. Indeed, references 4 and 39 
have already demonstrated that the calculated flow of reference 2 is in the correct 
direction. 

Since the only guide in choosing model sizes for the fan-in-wing tests of refer- 
ences 26 to 32 has been the set of limits given in reference 7, the data shown in figures 14 
and 18 should lead to serious concern regarding the highly favorable "fan-induced" lift 
reported as one of the main advantages to the fan-in-wing configuration in those studies 
which have produced and correlated uncorrected wind-tunnel data (e.g., refs. 19 and 26 
to 32). One such correlation (from ref. 19) is presented in figure 19, where the "fan- 
induced" lift is correlated as a function of the ratio of fan area to wing area. 

The present model has a ratio of fan area to wing area of 0.094 and, as may be seen 
in figure 14, has a "fan-induced" lift in the 1.12- by 2.24-m (44- by 88-in.) flat-oval 
insert which lies very near the lower boundary of the correlation region. This value does 
not really correspond to the other data in figure 19 because all those data were obtained 
with models having either no tail or a small tail; whereas the present model has an 
extremely large tail which carries a significant download (fig. 13) under almost all con- 
ditions. It is easily shown from the definition of ^ and Vj that 
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The value of C N T at a = 0° and v/Vj = 0.4 is obtained from figure 13(b) as 
-0.53; Srp^/Sp is 3.73;’ thus, from equation (12) the ratio N t /t s is equal to -0.16 for 
the conditions of figure 19. Removing the tail load from the data of figure 14 increases 
AL^ (= L - Tgj from 0.25 to 0.41, which is near the upper edge of the correlation band of 
reference 19. (See fig. 19.) On the other hand, in the largest test section, AL^Tg is 
negative for the complete model at this value of V/Vj. 

The slope of the correlation band of reference 19 which is reproduced in figure 19 
deserves some comment. It is obvious in examining the data points of figure 19 that a 
band of the same width, drawn parallel to the abscissa, and thus indicating total indepen- 
dence from the ratio Sjy/s^ would have encompassed a larger number of data points 
than the band which was drawn. In either event, the major exceptions to the correlation 
band are those configurations in which the fans are displaced far from the center of pres- 
sure of the wing. Such aircraft would be unflyable as VTOL configurations without the 
provision of additional fans to provide moment balance. 

Since reference 7 presents several different criteria upon which to scale wind- 
tunnel tests for wall interference, it is advisable to perform a first-order analysis in 
order to determine which parameters really are significant. For this first-order anal- 
ysis, examine the zero-angle-of-attack case, for which the wing of the present model 
would have no lift in free air. Then assume that the horizontal components of wall- 
induced interference have only a second- order effect and that Aa is sufficiently small 
to let Aa ~tan Aa. Under these assumptions, following references 1 to 3, 
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where Aa is in radians and where S„ r T and 6,„ are calculated for the fans. 
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From reference 23, momentum theory shows that for the fans at a = 0, 
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Substitution of equations (14) into equation (13) yields 


A a = - 





(15) 


Observe that both 6 T and 6 ^ are negative in a closed tunnel; thus, Aa 

WjIj W jD 

will be positive (upwash). Reference 23 shows that the lift of the fan is virtually insensi- 
tive to angle of attack for angles near zero; therefore, the increase in lift will be essen- 
tially all on the wing. This increase in lift may be written as 



(16) 


Substitute equation (16) into equation (15) to obtain 


AL = - ( 6 w,L + Vj 6 w,d)| C l/ 2s W ^ ± 


(17) 


Divide both sides of equation (17) by Tg = PSpVj to yield 
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Consider the product Cj^Sw in equation (18). Since Cl^ = this product may 
be rewritten as 
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Finally, substitute equation (19) into equation (18) to yield 
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which is the wall-induced lift. 
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Observe that the only term of equation (20) which explicitly involves the model 
dimensions is b^A-p For test sections having approximately the same width-height 
ratios (in the present case, from 1.4 to 2.0), b^/Arp will be approximately proportional 
to the square of the ratio of the wing span to the test- section width. The ratio of fan area 
to test-section cross-sectional area is completely immaterial. 

One more significant point is evident in the preceding analysis. It is generally 
believed that wall effects are greatest at low speeds. For a constant model configuration 
and for data presented in terms of coefficients based on free-stream dynamic pressure, 
wall effects are greatest at low speed. ^Note that ACl = AaC^ and, from eq. (15), 
that A a has a 1/V component.^ On the other hand, again for a constant model config- 
uration, when the data are presented in terms of forces or force ratios, equation (20) 
clearly shows that the greatest effect of wall interference will be at high speed. This 
conclusion is confirmed by the data presented in figure 14. 

It is obvious that the present results from the small insert lead to a gross overes- 
timate of "fan-induced" lift. The correlation with the data presented in reference 19 
indicates that the data presented therein also include substantial overestimates of "fan- 
induced" lift which would not be obtained in flight. (Observe that, with two exceptions, 
the models of ref. 19 have essentially the same span-to-width ratio as the present model 
in the smallest inserts. Of the two exceptions, one is anomalous because of its thin delta 
wing; the second was notable for producing the smallest "fan-induced" lift of any of the 
models of ref. 19.) Further, the model of reference 19 in which the fans are behind the 
trailing edge of the wing would be expected to show a far smaller "fan-induced" lift than 
indicated in figure 19, and, similarly, the model with the fans well forward of the wing 
would be expected to show far greater "fan-induced" losses than indicated. In either of 
these two cases, the results would be affected substantially by provision of the additional 
fans required for moment control in the VTOL mode. This latter effect is evident in fig- 
ure 19 when these two configurations are combined into one. (See the appendix of ref. 23 
for a further discussion of the effect of fan location on mutual interference.) Irrespective, 
however, of whether or not VTOL moment control is feasible for the configurations of ref- 
erence 19, it is obvious that all the data of that paper contain a large increment of wall- 
induced, rather than "fan-induced," interference. The "good" configurations will be far 
less "good" in free air; the "poor" configurations will be even worse in free air. 

It will be observed that there are differences in notation between the present paper 
and reference 19. In the present paper, ALj is defined (at a - 0°) as L - Tg since 
Tg is equal to the thrust in forward flight according to ideal momentum theory (ref. 23). 
Similarly, Vj is defined (see Symbols) as the fan efflux velocity in static thrust. In ref- 
erence 19, ALj is defined as total lift, less any wing lift (which is zero at a = 0° in 
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the present tests) , less the thrust in forward flight as measured by rakes in the fan exit; 
and Vj is defined in relation to this measured thrust. 

Because of the square root involved in determining Vj from the thrust, as well as 
the relatively flat character of l/ts near v/Vj = 0.4 and 0 = 0° (the conditions 
chosen by ref. 19; see fig. 14(b)), there will be little effect of the difference in definition 
of Vj. The difference in definition of ALj has a more serious effect. In practice, 
because of inlet efficiency, the actual value of thrust in forward flight will be somewhat 
less than the theoretical value of Tg. Figure 6 of reference 32 indicates that at v/Vj 
of 0.4, a loss of 10 to 15 percent of Tg may be expected for a typical lift-fan model. 

For complete comparability, this loss should be added to the present results; that is, in 
figure 19 the values of ALj^Tg should be about 0.1 to 0.15 greater than indicated therein 
for the present model. Thus, the effect of wall interference on the data of reference 19 
may be even greater than indicated by the data shown in figure 19. 

Effect of wall interference on drag .- The drag of the model also increases as the 
tunnel size decreases (figs. 12 and 15); however, the increases in drag are not commen- 
surate with the increases in lift. Indeed, at the lower angles of attack, the increases in 
drag are minimal. The disparity between the increases in lift and drag may be seen 
more clearly in figure 16, which presents the external drag-lift ratio for the model. At 
all angles of attack, D/L is greater in the larger test sections, thus indicating poorer 
efficiency. 

The apparent gain in efficiency in the smallest test sections retained even when 
the data are presented in terms of L (fig. 17). For example, at V^Vj = 0.4 and 
at a = 0° (fig. 17(b)), L/D E as measured in the smallest test sections is approxi- 
mately 25 percent greater than the same values measured in the largest test sections. 

Thus, wall effects are sufficient to indicate a 25-percent decrease in the power required 
to fly in the transition speed range. 

The values of l/Dj; shown in figure 17 appear at first glance to be remarkably 
small. They are confirmed however by the momentum theory presented in reference 23. 
This confirmation is demonstrated by the theoretical curves (from ref. 23) given in fig- 
ure 17. They are further confirmed by calculations made using the measured shaft powers 
given in reference 45, as well as by the extraordinary fuel consumption in low-speed flight 
found in design studies of fan-supported aircraft (ref. 46). 

Effect of wall interference on tail normal force .- The uncorrected measurements of 
tail normal force, as a function of V/Vj, are shown in figure 13. At low speed, the trends 
shown for the various test sections are observed to scatter. This effect is probably due 
to Rae’s limit (refs. 4, 9, 13, and 18), and it will be discussed in the next section. 
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At the higher speeds and for angles of attack less than 10° (figs. 13(a) to 13(c)), the 
observed tail normal-force coefficient is essentially independent of the test- section size 
or shape. As the angle of attack becomes greater, the data from the various test sections 
show greater differences (fig. 13(e)), with the tail normal force becoming more positive 
as the test- section size decreases. 

Even constancy in tail normal-force coefficient would indicate a serious degree of 
wall interference since, in free air, the increased lift (shown in fig. 14) in the small sec- 
tions would increase the downwash at the tail, reduce the tail angle of attack, and result 
in a more negative tail normal-force coefficient. However, the data of figure 13 indicate 
that the wall-induced interference at the tail is of sufficient magnitude to negate, or even 
to reverse, the trend that would be expected with increased lift. 

Wall-induced effects at the tail, of course, are not confined to this configuration. 
Large wall effects have also been noted in comparing large-scale wind-tunnel and flight 
data; for example, consider the comparison of flight- test data and uncorrected wind-tunnel 
test data (Ames full-scale tunnel test 388) presented in reference 20 for a YOV-10 air- 
craft fitted with a rotating- cylinder flap. Serious differences were found in maximum lift 
and the angle of attack at which it was obtained; however, by far the greatest disagree- 
ment between wind tunnel and flight was with regard to the effects at the tail. 

Figure 20 shows these differences (as presented in ref. 20) in terms of the elevator 
angle required to trim the aircraft as a function of forward speed. The uncorrected wind- 
tunnel data indicate positive speed stability with the stick moving rearward (the elevator 
moving trailing edge upward) as the speed decreases; the elevator is 20° trailing edge up 
when 55 knots is reached. In contrast, the flight data indicate a speed instability; the ele- 
vator angle is always in the opposite sense (trailing edge down); and at 55 knots the eleva- 
tor angle is 13° trailing edge down. The total disagreement between tunnel and flight at 
55 knots is 33°, and this disagreement is in the same direction as that indicated in 
figure 13. 

The trends shown in figure 20 are given further import by the flight-test data when 
extended to slightly lower speeds (fig. 21). Here the speed instability became more 
dramatic, and the minimum speed in many cases was determined by the speed at which 
the elevator contacted the limit of travel in the trailing-edge-down direction and not by 
maximum lift. Needless to say, under such circumstances, the pilot finds himself in 
somewhat compromised circumstances because he has no control left for any unanticipated 
maneuvering requirement. The point here, of course, is that not even full-scale wind- 
tunnel tests of the actual aircraft gave any indication that the pilot would find himself in 
these circumstances because wall effects were not properly accounted for in the data 
reduction. 
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It is noted that the YOV-IO with the rotating- cylinder flap, as tested in the tunnel, 
also fell well within the boundaries of reference 7, which, according to that paper, would 
indicate negligible wall effects (fig. 22). This evidence confirms the previous conclusion 
that the testing boundaries of reference 7 are erroneous. 

The opinion is sometimes voiced that wind-tunnel interference does not affect the 
trends shown by the data, or, as expressed in reference 47 (p. 7-1): "Informative results, 

even when the model lifting system spans 2/3 to 3/4 of the wind tunnel test section width, 
will still be obtained." Neither the model of the present investigation, nor the aircraft of 
reference 20, approached so great a size relative to the test section; and, in each case, 
the wall interference was so great that even the trends shown by the data were in the 
opposite direction from "free air" results at low speed. Such results clearly show that 
extreme caution must be used when interpreting uncorrected wind-tunnel data. 

Effect of test-section size on Rae's limit .- The separately instrumented tail was 
installed on the model in the hope that measurements of tail normal force would provide 
a sensitive indication of the onset of the recirculation which results in Rae's minimum- 
speed limit (ref. 9). This procedure was chosen because of the dramatic alterations in 
tail lift which were observed behind a rotor in reference 18. 

Although the tail normal-force- coefficient data presented in figure 13 do show 
marked effects as a function of tunnel configuration at low forward speed, effects as 
definitive as those of reference 18 were not always observed. One reason may be the 
magnitude of the wall interference in the present tests. This aspect of the problem will 
be discussed in subsequent sections of the present paper. A second reason is that the 
tail on this model, as may be seen by comparisons between the individual parts of fig- 
ure 13, has very nearly zero tail effectiveness ^that is, ~ 1, so that 1 - ~ 0 ) until 

the highest angle of attack (16°) is reached. At a = 16°, the tail normal-force coeffi- 
cient suddenly turns upward as the speed is reduced in the small test sections. In the 
small inserts, C N T departs from the trends shown in the data from the 9.14- by 18.3-m 
(30- by 60-ft) tunnel at a value of ^/Vj below 0.38; a similar departure may be observed 
in the data from the 2.13- by 3.05-m (7- by 10-ft) tunnel at a value of v/Vj below about 
0.2. The values correspond approximately to conditions at which visual tuft observations 
indicated substantial flow reversal on the floor; furthermore, the values are roughly in 
proportion to the height of the various test sections, as might be expected from the cor- 
relation rules presented in references 6, 9, and 21. Unfortunately, those rules are 
expressed in terms of the momentum wake angle. Since the momentum wake angle for 
the fan is always along its axis and is actually negative for the data of figure 13(e), those 
rules cannot apply to the present case. 

Tyler and Williamson (refs. 48 and 49) have conducted a systematic program to 
determine minimum-speed test limits for jet lifting systems. Their results indicate 
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incipient stagnation (near a = 0°) on the floor of the test section when 

^r =L59 TT (21) 

Vj h 

for single and tandem-paired jets, and when 

— = 1.31 — (22) 

V i h 

for a laterally paired system of two jets spaced 4.3 nozzle diameters apart. The spacing 
of the two fan nozzles in the present model is considerably closer (2.75 diameters); nev- 
ertheless, using equation (22) yields v/Vj = 0.67 for the two 1.12- by 2.24-m (44- by 
88-in.) test sections; V/Vj = 0.62 for the 1.22- by 1.83-m (4- by 6-ft) test section; and 
V/Vj = 0.39 for the 2.13- by 3.05-m (7- by 10-ft) test section. The corresponding val- 
ues for the largest test section are below the smallest velocities at which tests were run. 

To define the point of incipient stagnation and to define the point at which the data 
will be affected are two different things, as is noted in reference 49. Tyler and 
Williamson suggest that test speeds as small as 55 percent of the speed for incipient 
stagnation may be acceptable for single jets and 65 percent of this speed may be accept- 
able for widely spaced lateral pairs of jets. If the values obtained in the preceding para- 
graph are reduced by a multiplying factor of 0.6 (an average of 0.55 and 0.65), they will 
be observed to agree closely with the previously noted points of figure 13(e). Therefore, 
it would appear that the Tyler and Williamson relations (eqs. (21) and (22)) provide a rea- 
sonable means of estimating the minimum speed for wind-tunnel testing of jet- and fan- 
supported models. 

The value observed for the degradation of data due to recirculation in the 1.12- by 
2.24-in. (44- by 88-in.) flat-oval test section was about V/Vj = 0.38, and that obtained 
from equation (22) reduced by 40 percent was 0.40. Not only the correlation between 
these values is of interest; their magnitude is significant in itself. Observe that the cor- 
relation of "fan-induced" lift in reference 19 was obtained at v/Vj = 0.4. It is entirely 
possible that some of the data upon which reference 19 is based are suspect because of 
recirculation effects, since the model to tunnel-size ratios in those data are comparable 
with those obtained in the present small flat-oval insert. Furthermore, at 30 knots, it is 
clear that much of the data used to prepare the testing limits defined in reference 7 were 
obtained for flow conditions which were unrepresentative of flight in free air because of 
flow breakdown induced by the model in the wind tunnel. 
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Two primary requirements exist in planning wind-tunnel tests. One is simulation 
of the aircraft, and, given a drawing of the aircraft, it is simple to produce a reasonable 
model of it. Equally important, however, is that the basic free-air flow must also be 
simulated. At speeds less than Rae's limit, a powerful cylindrical sheet of vorticity is 
formed ahead of the intersection of the wake on the floor (refs. 4, 13, 18, and 39). This 
sheet ultimately extends across the floor and up the sides of the test section. Except in 
ground effect, no equivalent vortex formation exists in actual flight. Under such condi- 
tions, the flow in the test section does not simulate free air and almost any result may be 
obtained. 

It is particularly important to realize that the existence of the basic alteration of 
the flow does not depend upon the presence, or the absence, of a tail on the model. The 
flow alteration is caused by the presence within the walls of the main lifting system. 
Indeed, the models used by Rae (ref. 9) when he discovered this effect had no tails; neither 
did the models used by Tyler and Williamson (refs. 48 and 49). 

Comparison of simple momentum theory with experimental results .- Reference 23 
develops a simple incompressible-flow momentum theory for the fan-in-wing configura- 
tion based upon the assumption that there is no mutual interference between the fans and 
the wing. Momentum theory, by itself, is incapable of calculating the actual lift of the 
model, since the lift depends intimately on the local angles of attack of the wing and of the 
fan blades. However, once the lift is given, momentum theory is capable of estimating 
the remaining performance items. Momentum theory, obviously, also is incapable of 
predicting the tail normal-force coefficient because this coefficient depends upon a 
detailed calculation of the flow field in the vicinity of the tail. 

In the present case, it is assumed (following ref. 23) that the thrust of the fan is 
unaltered by forward speed or angle of attack. This assumption is verified by figure 6 
of reference 32, which shows that the actual thrust for a typical lift-fan model ^at v/Vj 
as great as 0.6^ is only 10 or 15 percent less than the static thrust. When the normal 
component of static thrust is added to the lift of the wing with covered, inoperative fans, 
the results previously presented in figures 11 and 14 are obtained. Evidently, at high 
speed, significant fan-wing interaction effects are present; however, throughout the usable 
transition speed range ^0 = V/Vj < 0.5^, the assumption of zero interaction yields values 
close to the observed total lift. The differences in notation between reference 19 and the 
present paper have no effect herein, since both theory and experiment are presented in 
the identical manner. 

All the remaining curves in figures 12 and 15 to 17 follow directly from the equa- 
tions of reference 23 once the lifts are assumed. It will be observed that, for transition 
speeds, the obsei’ved performance is predicted more closely by even this simple momen- 
tum theory than by the data from the small inserts. Note that the model in the present 
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investigation spanned only a little less than half the width of the smallest inserts. This 
relative size is essentially the same as that used in references 26 to 32, and a similar 
result may be implied to be true for those tests as well. 

Correcting the Data 

Considerations in correcting the data .- Correcting the data with the fans operating 
follows the same general procedure described earlier for the data with the fans covered. 
Obviously, the procedure is complicated to a degree in accommodating the presence of 
the fans. In this case, the interference factors are obtained from appendices O and P of 
reference 36.* The previously discussed modifications to subroutine DLTAS (appendix Q 
of ref. 36) were used to obtain the interference factors at the pitot-static tube location. 

The solid blockage factors are identical with those used when the fans were covered. 

Since there was no independent measurement of the thrust of each of the fans, there is no 
alternative but to deal with them simultaneously. Therefore, the appropriate interference 
factors for the pair of fans are the average of those for the fan due to its own presence 
and those for the fan due to the presence of the other fan. In all cases, slight changes to 
the programs of reference 36 allowed data cards containing the interference factors to be 
punched automatically as they were calculated. These cards were used as input data to 
the data correction program to be discussed shortly. This procedure eliminates the pos- 
sibility of errors in transcription when preparing the input to the correction program. 

Reference 36 offers choices of wing load distribution and rotor-disk load distribu- 
tion. In the absence of definitive measurements to the contrary, an elliptic load distribu- 
tion was chosen for the wing. In order to ascertain the degree to which this choice might 
affect the corrections, the data were also reduced using the interference factors for a 
uniform load distribution; no significant effect was found for this model, perhaps because 
of the magnitude of the corrections. Because of the large central boss in the fans, the 
disk load distribution over the faces of the fans is not uniform. Consequently, the trian- 
gular disk load distribution was taken as being more nearly representative of the actual 
fan load distribution. In any event, the fans were so small compared with the test-section 
dimensions that little effect of this choice should be evident. 

*Three known errors exist in the programs given in reference 36. Two of these 
affect the work contained herein. The following lines should be corrected to read as 
follows: 


805 XDELTA(Ll) =XDfcLTA<Ll)+UELTA< L1)*XL0AU(N1) 
SUML*0. 063052 
SUML s 0» 252208 


(E 79) 
(P 113) 
(P 136) 
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The first problem is to divide the measured loads between the elements which pro- 
duce them. The tail presents no problem since it was mounted on its own strain-gage 
balance; however, there was no balance to separate the independent forces of the wing 
and fans when they were operating in unison. In the absence of specific information, the 
wing was assumed to produce the same lift and drag as it did when the fans were covered. 
Thus, the fan lift and drag are assumed to be the main balance readings, less the mea- 
sured tail loads, and less the aforementioned assumed lift and drag of the wing. Then the 
fan lift and drag thus obtained were used to solve the momentum quartic and to calculate 
the fan wake skew angle x and the fan velocity ratio v/wg. The resulting values of x 
were within 2° or 3° of being equal to -a, as they should be (ref. 23). For the fans, it is 
more appropriate to use the effective skew angle 


*e 


_ 90° + X 
2 


(23) 


as given in references 8 and 39. For the wing and for the tail, the effective skew angle 
was assumed to be 90°. In the face of the powerful downwash field generated by the fans, 
the use of the momentum quartic given in references 2 and 22 is not strictly applicable 
because it would be necessary to include the local effective downwash angle in the vector 
diagram defining x- 1° any event, the use of x = 90° for the wing and the tail is a great 
simplification in the calculation. 

At this point, everything is in hand to compute the average interference velocity 
components over the wing, over the tail, and at the pitot tube. Observe that the fans now 
contribute substantially to the interference velocities at each location. 

Reference 13 has noted that it is necessary to apply a vortex-density correction to 
the theoretical interference factors when applying them to ducted fans. This correction 
was used in the present calculations. A justification of the vortex-density correction is 
presented in appendix A. 

The first step in applying corrections to the data from the smaller Ames tunnel is 
to correct the pitot-static tube reading of tunnel velocity. The forces that were charged 
to the wing are directly dependent on dynamic pressure; thus, it is necessary at this point 
to return to the original division of loads and redo that division with the corrected dynamic 
pressure, and then repeat all the steps to this point. 

The calculation now proceeds as before, correcting the overall performance to the 
corrected flight condition at the wing and then adjusting the tail loads to account for the 
substantially different wall-induced interference at the tail location. Despite the fact that 
the fans are mounted within the wing planform, there is a significant difference in the 
average wall-induced interference over the full span of the wing and the similar average 
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over the faces of the fans. It is necessary to remove this difference by adjustments to 
the fan lift and drag. The adjustment is accomplished by the use of the following equa- 
tions from reference 23: 


AD = 



+ Aip 


cos a 


(24) 


AL = -Tg Aip sin a 


(25) 


where Aip is in radians. It will be observed from figure 15 that an a-dependent multi- 
plying factor, generally greater than 1.0, could have been applied legitimately to equa- 
tion (24) (that is, the actual drag of the fans is generally greater than the momentum- 
theory value). This was not done; equation (24) was used directly as given above. 

Corrections at zero tunnel velocity are particularly suspect. As noted in refer- 
ence 8, a hovering condition in the tunnel leads to an interference which is a pure upwash. 
Proceeding in a formal manner, this upwash is equivalent to change in angle of attack of 
90°. While true in a sense, it is more rational to consider the model to be at the same 
angle of attack, but with a rate of sink equal to the upwash velocity. The corrections to 
the data then would depend upon the effect of a rate of sink. Unfortunately, sufficient data 
are not present to make such a correction for this model. Furthermore, at zero speed, 
the test conditions always violate Rae's minimum- speed limit (ref. 9). 

In view of the foregoing observations, no attempt has been made herein to correct 
data obtained at zero velocity. When such points are shown in the corrected data, they 
are identical with the uncorrected data, and they are presented only to preserve the con- 
tinuity of the data set. 

Computer program.- Appendix B presents one of the computer programs used in 
correcting the data obtained in the present investigation. Because of differences in the 
measurement of fan rotational speed, it was necessary to write slightly different pro- 
grams for the two separate sets of tests in the Ames tunnel. The absence of a stream- 
angle correction in that tunnel provides one simplification in that it was possible to com- 
pute the interference factors specifically for the angles of attack which were used; thus, 
only a single interpolation against x e i s required to obtain the proper factors for each 
data point. The presence of a small stream angle in the Langley 2.13- by 3.05-m (7- by 
10-ft) insert was accommodated by means of a double interpolation against both x an d 
a in a third modified version of the computer program. 
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The program of appendix B is substantially more complex than would normally be 
required because it simultaneously treats four different test sections. This feature 
requires considerably more storage and additional steps (which increase running time) 
than would a program written for a single test section. Nevertheless, in the Langley 
Computer Complex, the program requires only 54000g (=22 500iq) storage locations for 
compilation, 460008 (~19 500io) storage locations to run, and completely corrects more 
than 360 data points in about 30 seconds (including compilation time) at a cost of only $7. 
The storage lengths, the time, and the cost obviously would be different in almost any 
other computer; however, it also is obvious that only minimal costs and computer capa- 
bilities are required to fully correct data for wall effects even for a fairly complex model. 

The program of appendix B produces several files of output data and a sequence of 
punched-card sets for subsequent plotting of the data. In sequence, the written files pre- 
sent the interference factors used in the correction routines, the uncorrected data 
together with a preliminary breakdown of the division of loads (in the Langley insert, the 
presence of a stream angle requires interpolation in order to obtain this file), a point-by- 
point listing of the corrected data, a listing of the corrected data at fixed angles of attack 
obtained by interpolation of the previous listing, an interpolated listing of the corrections 
themselves at a series of fixed corrected angles of attack, and finally a listing of the cor- 
rections according to the uncorrected angles of attack. Punched-card decks of the last 
four listings are provided for subsequent automatic plotting of the data. 

The data herein are presented as a function of forward speed. If it is desired to 
obtain polar plots of the performance at fixed speeds, an interpolation against v/Vj 
would be required. The addition of one more interpolation should not present any signifi- 
cant difficulty. 


Interference in Uncorrected Data 

The corrections in the uncorrected data of figures 13 to 17, as calculated from the 
foregoing considerations, are presented in figures 23 to 28. The corrections are dis- 
tinguished by their enormous magnitude, which far exceeds the more reasonable values 
suggested as the maximum practical limits in references 6 and 42. Depending upon a 

and V/Vi, the average interference angle Aa at the wing varies from about 2^° to over 
J 2 
14 in the smaller inserts (fig. 23). Similarly, the effective dynamic pressure at the wing 

is reduced by 5 to 22 percent (fig. 24). The effective tail incidence in the smaller inserts 

is increased by from 5° to 12° (fig. 25) and the dynamic pressure at the tail varies from 

1.15 to almost 3 times that at the wing (fig. 26). In the small inserts, even the fans are 

operating at effective angles of attack as much as 7° more than the wing in which they 

were mounted (fig. 27). Only the ratios of the dynamic pressures at the fans to those at 

the wing remain relatively small (fig. 28). The wall interference in the 2.13- by 3.05-m 
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(7- by 10-ft) test section is generally of a lesser magnitude, although, even in that test 
section, Aar and Ai-j. (figs. 23 and 25) tend to be larger than would be desirable. 

It is noticeable that the values for Rae's limit, which were obtained earlier, were 
for speeds so low that a prudent investigator would have discontinued testing long before 
this limit was a serious concern. This result is in accordance with the results of refer- 
ence 6, which show that Rae’s limit is of primary concern only for those models which 
are very small with respect to the test-section size. In all other cases, the magnitude 
of the wall-induced distortions of flow over the model are the controlling factors in 
choosing the maximum permissible model size. 

Corrected Data 

It is obvious from the magnitude of the indicated corrections (figs. 23 to 28) that 
perfect correlation of the data from all the test sections cannot be expected. Neverthe- 
less, the corrected data show remarkably improved agreement (figs. 29 to 35). This 
agreement is poorest at those speeds previously determined to be less than Rae’s limit 
(v/Vj ~ 0.38 in the smallest test sections^. It is also somewhat poorer in the smaller 
inserts than in the more moderately sized 2.13- by 3.05-m (7- by 10-ft) test section. 

(See particularly fig. 33(a).) On the other hand, considering the magnitude of the required 
corrections, the data presented in figures 29 to 35 are an impressive verification of wall- 
interference theory (refs. 1 to 6). 

Interference in Corrected Data 

It is obvious that the angular range of model settings in figures 29 to 35 differs from 
that in figures 11 to 17 because the geometric angle of attack in figures 11 to 17 has been 
decreased by Aa to obtain the corrected values of a in figures 29 to 35. Since the 
wall-induced interference increases as the lift and drag increase, and since these forces 
increase with a, the corrections are somewhat less in figures 29 to 35 than in figures 11 
to 17. For completeness, figures 36 to 41 have been prepared to indicate the magnitude of 
the corrections that are actually present in the corrected data. Although slightly less than 
the values presented earlier, the corrections in the final data are still extremely large. 

Maximum Model Size 

It is clear from the data presented in this paper that the model was so large in rela- 
tion to the small inserts that the corrections were excessive. At low speeds, the model 
was excessively large in the 2.14- by 3.05-m (7- by 10-ft) test section as well. It would 
appear that prudent model sizing would have led to a model having a span of about a quar- 
ter of the test- section width. 
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A similar conclusion as to an appropriate model size would have been reached by 
examining the charts of reference 6. The overall corrections (a. a and q c /q] and the 
wall-induced tail incidence ^Ai-pj are about the same in those charts as were found herein; 
however, the charts of reference 6 fail totally to indicate the magnitude of the wall- 
induced dynamic-pressure ratio ^q-jy/q c j at the tail. Such a discrepancy is understand- 
able. The charts of reference 6 are based on the assumption of a winglike model having 
a uniformly loaded span and a single, blended wake; these assumptions are grossly vio- 
lated by the present fan-in-wing model. Thus, in using reference 6 to size models of 
unusual VTOL configurations which do not approximate the assumptions of that paper, it 
is best to err on the safe side by selecting a model size even smaller than indicated 
therein. 

It will be observed that nonuniformities in wall-induced interference decrease more 
rapidly (ref. 6) with decreases in span (in a given tunnel) than do the overall corrections, 
which are roughly proportional to the square of the span. Thus, a small model does not 
require the same rigor in applying corrections as is required for a large-span model. 
Furthermore, considerably more confidence in the final results is justified when the 
model is small enough to require only minimal correction to the data. 

While the present model should have been about half its present size, in the smallest 
insert, this conclusion should not be extended to indicate that all VTOL models should 
span about one-quarter of the test-section width. The allowable size of a VTOL or STOL 
model will depend upon the configuration, upon the minimum speed for which useful data 
are required, and upon the degree of correction applied to the data. Reference 6 should 
be some help in this regard; however, the only real safety will be found in correcting, as 
fully as possible, all wind-tunnel data as a standard practice. As noted earlier, the addi- 
tional computing cost is minor in comparison with the total cost of a wind-tunnel 
investigation. 


CONCLUSIONS 

The results of this investigation of wind-tunnel wall interference on the performance 
of a fan-in-wing model are as follows: 

1. Extreme caution must be used in interpreting uncorrected wind-tunnel data 
obtained at low speeds. Unless the model is extremely small in relation to the test- 
section size, the wall interference can be so large that even the trends in the data may be 
opposite to those which would be obtained in flight. 

2. Wall-induced interference is particularly large at the model tail. This result 
confirms recently published (NASA CR-2135) conclusions based on the correlation of flight 
and wind-tunnel data for a YOV-IO aircraft. 
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3. The theory of wall interference for VTOL and STOL models, presented in NASA 
TR R- 124 and subsequent papers, has been verified under conditions of extreme wall 
interference. 

4. The rules for choosing model sizes to produce negligible wall effects, as given 
by Cook and Hickey in NASA SP-116 (also in AGARD Rep. 520), appear to be considerably 
in error and to permit the use of models which are significantly too large for the tunnel. 

5. The method presented by Tyler and Williamson in AGARD CP-91-71 yields rea- 
sonable estimates of the onset of Rae's minimum-speed limit for jet- and fan-supported 
models; however, for reasonably large models, wall interference becomes so great that 
testing should be discontinued at a speed significantly greater than Rae's limit. 

6. The "fan-induced" lift indicated by a number of previous investigations appears 
to be largely the result of wind-tunnel wall interference which was not accounted for in 
reducing the data. The uncorrected results obtained herein, when the model spanned 
almost half of the tunnel width, fall directly on a previously published correlation of "fan- 
induced ' lift (Hickey and Cook, AGARD CP-22, paper No. 15); however, the increase in 
lift for the model under conditions which approach testing in free air was small or nega- 
tive for the actual transition speed range. 

7. The simple incompressible-flow momentum theory presented in NASA TN D-7498 
appears to yield reasonable estimates of fan-in-wing performance in transition flight; 
indeed, the theoretical predictions are more accurate than uncorrected wind-tunnel test 
data in which the model span is approximately half of the tunnel width. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., December 12, 1973. 
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APPENDIX A 


JUSTIFICATION OF VORTEX-DENSITY CORRECTION 
FOR FANS AND JETS 

The theoretical treatment of wall interference for VTOL-STOL aircraft (refs. 1, 2, 
and 5) sets up the inclined wake of the aircraft in free air as a doublet string extending 
from the aircraft to infinity. If the model is small, this doublet string might represent a 
rotor or a lifting fan or a jet. Indeed, the relationship between the doublet strength and 
the induced velocity was obtained directly from an earlier analysis of a rotor wake in the 
wind tunnel (ref. 50). The basis of the relationship was the doublet strength required to 
match the vorticity along the edge of the vortex cylinder comprising the wake. 

Now, if one considers a vortex cylinder cutting through an otherwise unbounded 
flow, and then takes the line integral of V • ds (where V is the total vectorial velocity 
and s the path length), the eventual result is that the vorticity along the edge of the cyl- 
inder is precisely equal to the velocity jump across the cylinder. Unfortunately, integra- 
tion, by means of the Biot-Savart law, of all the vorticity in the wake, leads to a velocity 
only one-half this great at the origin of the cylinder. In order to obtain the correct veloc- 
ity at the end of the cylinder, it is necessary to double the vorticity. Since the correc- 
tions for wall interference depend upon the velocity Wq, which is the mean vertical 
induced velocity, it appears appropriate to take the corrections due to the fan as being 
twice as great as those of references 1 to 3. The changes need not be made in the inter- 
ference factors, but can be made most simply by doubling the wall-induced interference 
components caused by the presence of the fans. 

This effect was first noted in reference 13 and the vortex-density correction was 
used both in that paper and in the present analysis. The results of both papers appear to 
justify its use. 

It will be observed that no similar correction is required for rotors or propellers. 
In those configurations, the induced velocity at the origin of the wake should be equal to 
one-half the vortex density. 
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APPENDIX B 


FORTRAN PROGRAM FOR CORRECTING DATA 
TESTED IN FOUR DIFFERENT 


FROM A FAN-IN-WING MODEL 
TEST SECTIONS 


THIS PROGRAM WAS WRITTEN IN CDC FORTRAN, VERSION 2.1, TO RUN ON THE CDC 
6000 SERIES COMPUTERS IN THE LANGLEY RESEARCH CENTER COMPUTER COMPLEX. MINOR 
MODIFICATIONS MAY BE NECESSARY PRIOR TO USE IN OTHER COMPUTERS. A DESCRIPTION 
OF THIS PROGRAM IS GIVEN IN THE TEXT OF THIS PAPER. A COMPLETE LISTING OF THE 
INTERFERENCE FACTORS USED IS INCLUDED. EACH LINE IS CODED AT THE END BY THE 
ANGLE OF ATTACK FOR WHICH THE FACTORS WERE COMPUTED, AN INTEGER CODE SPECIFYING 
THE TEST SECTION, AND THE CODE WORD DESCRIBED WITHIN THE PROGRAM. THESE INTER- 
FERENCE FACTORS WERE OBTAINED USING THE COMPUTER PROGRAMS OF NASA TM X-1740 IREF 
361. CERTAIN ERRORS IN THAT REFERENCE ARE OISCUSSED IN THE TEXT. 


THE SUBROUTINE OISCOT IS INCLUOED FOR COMPLETENESS. THIS IS A RELATIVELY 
STANDARD SINGLE OR DOUBLE INTERPOLATION ROUTINE. IT, OR ITS EQUIVALENT, WILL 
BE FOUND IN MOST COMPUTER SYSTEM LIBRARIES. 


PROGRAM AARL6A (I NPUT , OUTPUT ,TAPE5= INPUT , TAPE6=0UTPUT , T APE 1 , TAPE3 , (B 1) 
1 TAPE4, PUNCH) (B 2) 

C 


c 


TUNNEL 

CODE 




c 

I TUN* 1 

IS 

44X88 

INCH 

WITH 

ROUND ENDS 


c 

I TUN= 2 

IS 

44X88 

INCH 

WITH 

RECTANGULAR 

ENDS 

c 

I TUN=3 

IS 

48X72 

INCH 

WITH 

RECTANGULAR 

ENOS 

c 

I TUN = 4 

IS 

7X10 

FOOT 

TUNNEL 



C 

C FORTRAN WORDS REPRESENTING INTERFERENCE FACTORS ARE ALL CODED BY 

C THE LAST FOUR CHARACTERS OF THE WORD. STARTING FROM THE RIGHT- 

C HAND SIDE OF THE WORD, THE FIRST CHARACTER REPRESENTS THE ELEMENT 

C ACTED UPON AND THE SECOND CHARACTER THE ELEMENT WHICH CAUSES THE 

C WALL INTERFERENCE, WHERE: W=WING; F = FANS; T= TA I L ; AND P = PITOT. 

C THE NEXT TWO CHARACTERS ARE THE SUBSCRIPTS OF THE INTERFERENCE 

C FACTORS AS DEFINED IN NASA TR R-124. VARIOUS PREFIXES ARE 

C APPENDED TO THESE CODE LETTERS TO DISTINGUISH SPECFIC CHOSEN 

C VALUES. 

C 

REAL LOTS 

DIMENSION WLFF(4,6,8),ULFF(4,6,8),WDFF(4,6,8) ,UDFF(4,6,8 ), 

1 WLWF(4,6),ULWF14,6) , WLWM4.6) ,ULWW(4,6> , WLFW I 4 , 6 , 8 ) , ULFW (4 , 6 ,8 ) , 

2 WDFWC4,6,8> ,UDFW<4,6,8) , WLF T ( 4, 6 , 8 ) , UL FT I 4, 6 , 8 ) , WDFT ( 4, 6 , 8 ) , 

3 UDFT<4,6,8) ,WLWT(4,6) ,ULWT(4,6) ,WLTT(4,6> , ULTT ( 4 ,6 ) , ANAMFI 28 ) 
DIMENSION CWLFFI8) , CUL FF ( 8 ) , CWDF F ( 8 ) , CUDFF ( 8 ) , CWLFW ( 8 ) ,CULFW(8I, 

1 CWDFW ( 8 ) , CUDFW I 8) ,CWLFT(8) ,CULFT( 8) ,C WDFT 18) ,CUDFT!S) 

DIMENSION TDALPFI6) ,TDALPT (6) * TQOQF I 6) , TQ0QTI6) ,TDVOVJ( 6 ) ,TALP(6) , 

1 TO ( 6 » , TDALPW ( 6 » , TQOQWI 6 ) ,TQC<6) ,TALPC( 6) , TCNC ( 6 ) , TL I FT ( 6 ) , SLT ( 6 ) , 

2 TVOVJI 6) »TLQTS(6 » , TD0H6) , A ( 5 ) , TQOQ J( 6 ) , EPS I LON ( 4 ) , BDT ( 6 » , BLTI 6 ) , 

3 CHI (8) ,AT(4),SLAT(4,6) , SO AT ( 4 , 6 ) , JTUN I 4 ) , J TY PE ( 3 ) ,SLAT1 (6) , 

4 SLAT2I6) , SLAT3(6),SLAT4(6),SDAT1(6),SDAT2(6),SDAT3I6),SDAT4(6), 

5 UOVPIT (4),W0VPIT(4) 


IB 3) 
IB 4) 
IB 5) 
IB 6) 
IB 7) 
IB 8) 
IB 9) 
(B 10) 
IB 11) 
IB 12) 
IB 13) 
(B 14) 
IB 15) 
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APPENDIX B 


C 

C 

C 

C 

C 

C 


C 

C 

C 

C 

C 

c 

c 

c 


c 

c 

c 


c 

c 

c 


DIMENSION WLF P 14,6*8) , ULFP(4«6,8), WOFP 14, 6, 8 ) *UDFP (4*6 ,8 ) t 

1 WLWP (4*6) tULWP (4,6 ) t WDWP( A. 6) ,UDWP (4,6)*CWLFP(8 ) ,CULFPI 81 . 

2 CWOFP( 8) tCUDFPl 8) 

NOTE THAT WLFP,ULFP,WDFP,UDFP,WLWP,ULWP,WOWP, AND UDWP ARE ALL 
COMPUTED RETAINING THE N=M=0 TERMS IN THE WALL-EFFECTS CALCULA- 
TION. THUS, THEY CONTAIN THE DIRECT EFFECT OF THE MOOEL FLOW 
FIELD AS WELL AS THE WALL EFFECTS AT THE PITOT LOCATION. 

OAT A ( AID, 1 = 1, 0. *5., 10. , 16./ 

DATA (ATI II ,1=1, 41/24. 004 *26. 889, 24. ,70./ 

DATA (CHI(I), 1 = 1, 81/20. ,30., AO. , 50. ,60. ,70. ,80. ,90./ 

DATA ( JTUNUI , I =1 , A » / 10H2-1 ROUND ,10H2-I RECT. , 10H 1.5-1 , 

1 10H 7X10 / , 

DATA ( JTYPEI I I , I=1*3)/10H WT TARE *10H AERO TARE*10H OATA / 

DATA ( E PS I LON ( I ), 1*1, A)/. 0 27 05, .02AA5,.02565* .00716/ 

UOVPIT AND WOVPIT ARE THE SOLIO BLOCKAGE EFFECTS (INCLUDING THE 
DIRECT MODEL FIELO) OF THE BODY AT THE PITOT LOCATION. THESE ARE 
COMPUTED FROM A SOURCE AND SINK SPACED SO AS TO PRODUCE IN FREE 
AIR A RANKINE OVOIO OF THE SAME LENGTH AND DIAMETER AS THE MODEL 
FUSELAGE. IN THE TUNNEL THE BODY WILL BE SOMEWHAT SLENDERER 
AND THE CALCULATED BLOCKAGE WILL BE SOMEWHAT UNDERESTIMATED. 

DATA ( UOVPIT ( 1 1,1 = 1 ,4 I / -.0 3672,-. 00707,-. 00443, 0.00035/ 

DATA ( WOVP IT ( I I, 1 = 1 , AI/2 *0.00685, 0.00622, 0.CC191/ 

STATIC WEIGHT TARE DATA 

OATA (BOT< I ) , 1=1,6)/- 1.7 38 3, -0.8929,-0.02 38,0. 8546,1.7515,2.7357/ 
DATA (BLT( I 1 ,1 = 1, 61/0.15 66, 0.0113, 0.0081, 0.0197, 0.0835, 0.3 73 2/ 

DATA (SLT(I), 1=1, 61/0.1999, -0.07 50, 0.22 50,0. 1000,0.0333,0.0374/ 

AERODYNAMIC TARE DATA 

DATA ( SLAT 1( II , 1=1,6 )/-. 477,-. 200, .026, . 258,. 518. . 897/ 

DATA ( SLAT2( I ) , 1 = 1 , 6 I /- . 482 , -. 207 , .02 1 , . 2 51 , . 523, . 902 / 

DATA (SLAT3I I), I =1,6)/-. 433,-. 179,. 036,. 2 56,. 508,. 882/ 

DATA (SLAT4( I), 1=1, 61/-. 390, -.185, .010, .205, .420, .697/ 

DATA (S DA TUI), 1 = 1, 6)/. 198, .152, .137, .144, .178, .273/ 

DATA ( S DAT 2 ( I), 1 = 1, 6)/. 1 97, .153, .138, .146, .179, .274/ 

DATA (S DA T3(I), 1=1, 6)/. 182, .140, .128, .135, .169, .258/ 

DATA (S DAT 4(1), 1=1, 6)/. 2 08, .164, .151, .155, .182, .246/ 

DO 6 1=1,6 

SLAT(1,I)=SLAT1( I) 

SL AT ( 2 , I ) =SL AT2 ( I ) 

SL AT ( 3 , I ) =SL AT3 ( I ) 

SLATC4, I )=SLAT4( I ) 

SD AT ( 1 , I ) =SDAT 1 ( I ) 

SDAT (2,1) =S0AT2( I ) 

SD AT ( 3 , I ) =S0AT3 ( I ) 

6 SDAT(4,I)=S0AT4(I ) 

PI=3. 14159265258979 
PAD=PI /180. 

SW=7. 41125 
ST =2 . 6042 
SF=2.*PI/9. 

REWIND 1 
REWIND 3 
REWIND 4 


B 16) 
B 17) 
B 18) 


B 19) 
B 20) 
B 21) 
B 22) 
B 23) 
B 24) 
B 25) 


IB 26) 
(B 27) 


(B 28) 
(B 29) 
IB 30) 


IB 31) 
(B 32) 
IB 33) 
IB 34) 
(B 35) 
(B 36) 
IB 37) 
(B 38) 
(B 39) 
(B 40i 
IB 41) 
IB 42) 
IB 43) 
IB 44) 
(B 45) 
IB 46) 
IB 47) 
IB 48) 
IB 49) 
IB 50) 
IB 51) 
IB 52) 
IB 531 
IB 54) 
IB 55) 
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no o o o o 


APPENDIX B 


REAO IN INTERFERENCE FACTORS 


WRITE 

( 6* 1 20 ) 


(B 

56) 

READ 

( 5,2024) 

WINGLDG 

(B 

57) 

WRITE 

(6,2025) WINGLDG 

(B 

58) 

DO 20 

I TUN= 1 

,4 

(B 

59) 

WRITE 

(6,127) 

JTUN( I TUN ) 

(B 

60) 

DO 20 

I ALPHA 

*1,6 

(B 

61) 

READ 

(5,121) 

( WLFF ( I TUN , I ALPHA, ICHI ) , ICH I = 1 , 8 ) , ANAME ( 1 ) 

(B 

62) 

IF (EOF, 5) 999,22 

( B 

63) 

22 REAO 

(5,121) 

( ULFF ( ITUN, (ALPHA, ICHI ) , ICH I *1 , 8 ) , ANAME ( 2 ) 

(B 

64) 

READ 

(5,121) 

( WOFF ( I TUN, I ALPHA, ICHI 1 , ICH I = 1 , 8 » , ANAME ( 3 ) 

(B 

65) 

READ 

(5,121) 

( UOFF ( ITUN , I ALPHA , ICHI ), ICHI = 1,8), ANAME (4) 

(B 

66) 

READ 

(5,122) 

WLWF ( ITUN, I ALPHA ) , ANAME ( 5 ) 

(B 

67) 

READ 

(5,122) 

ULWF ( ITUN, (ALPHA) , ANAME 1 6 ) 

(B 

68) 

READ 

(5,2022 ) 


(B 

69) 

READ 

(5,2022) 


(B 

70) 

READ 

(5,122) 

WLWW ( I TUN , I ALPHA), ANAME (7) 

(B 

71) 

READ 

(5,122) 

ULWW( ITUN, I ALPHA) ,ANAME(8> 

(B 

72) 

READ 

(5,2022) 


(B 

73) 

READ 

(5,2022) 


(B 

74) 

READ 

(5,121) 

(WLFWI ITUN, I ALPHA, ICHI ) , ICH 1 = 1, 8 ) , ANAME ( 9 ) 

(B 

75) 

READ 

(5,121) 

(ULFWI I TUN, I ALPHA, ICHI ) , I CH I = 1 , 8 ) , ANAME ( 10 ) 

(B 

76) 

READ 

(5,121) 

( WDFW( ITUN, I ALPHA, ICH I ) , ICH 1= 1 , 8 ) , AN AME( 1 1 ) 

(B 

77) 

READ 

(5,121) 

( UOF W ( I TUN, I ALPHA, ICHI I , ICH 1 = 1 , 8 ) , ANAME ( 12 ) 

(B 

78) 

READ 

(5,121) 

( WLF T ( ITUN, I ALPHA, ICHI 1 , ICH 1 = 1, 8 ) , ANAME! 13 ) 

( B 

79) 

READ 

(5,121) 

( ULF T ( I TUN, I ALPHA, ICHI ) , ICHI =1 , 8 ) , ANAME ( 14 ) 

( B 

80) 

READ 

(5,121) 

(WOFTI ITUN, I ALPHA, ICHI ) , I CH I = 1 , 8 ) , AN AME ( 1 5 ) 

(B 

81) 

READ 

(5,121) 

( UDFT ( ITUN , I ALPHA, ICHI), ICH 1 = 1 , 8 ), ANAME ( 16 ) 

(B 

82) 

READ 

(5,122) 

( WLWT( ITUN, I ALPHA) , ANAME ( 17) ) 

( B 

83) 

READ 

(5,122) 

(ULWT( I TUN, I ALPHA) , ANAME ( 18) ) 

(B 

84) 

REAO 

(5,2022) 


(B 

35) 

READ 

(5,2022) 


(B 

86) 

READ 

(5,122) 

( WLTT( ITUN, I ALPHA) .ANAME ( 19) ) 

(B 

87) 

REAO 

(5,122) 

(ULTT( I TUN, I ALPHA) , ANAME ( 20) ) 

(B 

88) 

READ 

(5,2022) 


(B 

89) 

READ 

(5,2022) 


(B 

90) 

READ 

(5,121) 

( WLFP ( ITUN, I ALPHA, ICHI ), ICH 1 = 1 , 8 ), ANAME < 21 ) 

(B 

91) 

READ 

(5,121) 

( ULFP( I TUN, I ALPHA, ICHI ) , I CH I =1 , 8 ) , ANAME ( 22 ) 

(B 

92) 

READ 

(5,121 ) 

(WDFPI ITUN, I ALPHA, ICHI ) , I CH I = 1 , 8 ) , ANAME ( 23 ) 

(B 

93) 

REAO 

(5,121) 

(UDFPl I TUN, I ALPHA, ICHI ) , ICH 1 = 1 , 8 ) , ANAME ( 24 ) 

(B 

94) 

READ 

(5,122) 

( WLWP ( I TUN, I ALPHA) .ANAME (25) ) 

(B 

95) 

READ 

(5,122) 

( ULWP ( I TUN, I ALPHA) , ANAME ( 26) ) 

(B 

96) 

READ 

(5,122) 

( WDWP ( ITUN, I ALPHA) .ANAME (27) ) 

(B 

97) 

READ 

(5,122) 

(UOWP( ITUN, I ALPHA) , ANAME ( 28) ) 

(B 

98) 

WRITE OUT INTERFERENCE FACTORS 



WRITE 

(6,1231 

(CHIII), 1=1, 8) 

(B 

99) 

WRITE 

(6,1231 

( WLFF ( ITUN, I ALPHA, ICHI ), ICHI= 1, 8 ), ANAME ( 1 ) 

(B 

100) 

WRITE 

(6,123) 

( ULFF ( ITUN, I ALPHA, ICHI ) , ICHI = l , 8 ) , ANAME ( 2 ) 

(B 

101) 

WRITE 

(6,123) 

(WOFF ( ITUN, I ALPHA, ICHI ) , ICH 1= 1 , 8 ) , ANAME ( 3 ) 

(B 

102) 

WRITE 

(6,123) 

(UOFF! ITUN, I ALPHA, ICHI ) , ICH 1 = 1 , 8 ) , ANAME ( 4 ) 

(B 

103) 

WRITE 

(6,124) 

( WLWF ( ITUN, I ALPHA), ANAME! 5) ) 

(B 

104) 

WRITE 

(6,124) 

( ULWF ( ITUN, I ALPHA) ,ANAME(6) ) 

(B 

105) 

WRITE 

(6,124) 

(WLWW( ITUN,! ALPHA) , ANAME! 7) ) 

(B 

106) 

WRITE 

(6,124) 

(ULWW! I TUN, I ALPHA) , ANAME ( 8 ) ) 

(B 

107) 

WRITE 

(6,123) 

(WLFW( ITUN, I ALPHA, ICHI ) , ICH 1 = 1 , 8 ) , ANAME ( 9 J 

(B 

108) 

WRITE 

(6,123) 

(ULFW( I TUN, I ALPHA, ICHI ) , ICHI = 1 , 8 ) , ANAME ! 10 ) 

(B 

109) 
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appendix b 


c 

c 

c 


c 

c 

c 

c 

c 

c 


c 

c 

c 


c 

c 

c 


WRITE (6,123) ( WDE W ( ITUN,I ALPHA f I CHI ) , ICHI=1 , 8 ) , ANAME ( 11 ) (B 

WRITE (6,123) (UDFWUTUN, IALPHA, ICHI) , ICH 1= 1, 8 ) , ANAME ( 12 ) (B 

WRITE (6,123) (WLFTdTUNtl ALPHA, ICHI), ICHI=1 , 8 ) , ANAME ( 13 ) (B 

WRITE (6,123) (ULFT( ITUN , I ALPHA, ICHI), ICHI = 1 , 8 ) , ANAME ( 14 ) (B 

WRITE (6,123 ) ( WOET( ITUN, I ALPHA, ICHI ) , I CH 1= 1 , 8 ) , AN AME ( 1 5 ) ( B 

WRITE (6,123) (UDFT( ITUN, IALPHA, ICHI), ICHI*1, 8), ANAMEU6) (B 

WRITE (6,1241 WLWT( ITUN, I ALPHA) ,ANAME( 17) (B 

WRITE (6,124) ULWT( ITUN, IALPHA), ANAME(18) (B 

WRITE (6,124) WLTK ITUN, I ALPHA) ,ANAME< 19) (B 

WRITE (6,124) ULTT( ITUN, I ALPHA) ,ANAME( 20) (B 

WRITE (6,123) (WLFP( I TUN, I ALPHA, ICHI ) , I CH 1 = 1 , 8 ) , AN AME ( 2 1 ) (B 

WRITE (6,123) < ULFP ( ITUN, I ALPHA, ICHI ) , IC H I = 1 , 8 ) , ANAME ( 22 ) ( B 

WRITE (6,123) ( WDFP ( I TUN, I ALPHA, ICHI ) , I CH 1= 1 , 8 ) , AN AME ( 23 ) (B 

WRITE (6,123) ( UDFP ( ITUN, I ALPHA, ICHI ) , I CH 1= 1 , 8 ) , ANAME ( 24 ) ( B 

WRITE (6,124) (WLWP( ITUN, I ALPHA) ,ANAME(25) ) (B 

WRITF (6,124) (ULWP( ITUN, IALPHA), AN AM E< 26)) (B 

WRITE (6,124) (WDWP( ITUN, l ALPHA) ,ANAME(27) ) IB 

WRITE (6,124) (UDWP( ITUN, IALPHA) ,ANAME(28) ) (B 

WRITE (6,2022) (B 

20 CONTINUE (B 

wRITE (6,110) (B 

L I NE =0 (B 

READ IN TEST DATA 

1 RF AD (5,100) IRUN, IRPM, I TUN, ITYPE, l FR A ME , Q , B ALDR AG , BALLI FT, SCLDRAG (B 
1 , SCLLIFT, ALPHA, RPM, DENS ITY <B 

IF ( EOF , 5 ) 998,2 (B 


NOTE THAT PITOT MEASUREMENT OF TUNNEL VELOCITY IS CORRECTED ON 
THE FIRST PASS THROUGH THE INTERFERENCE CALCULATIONS* ON THE 
SECOND PASS, THE LOADS ARE REDIVIDED ACCORDING TO THE CORRECTED 
PITOT READING AND ALL THE CORRECTIONS ARE CALCULATED. 


2 I TR I P= 1 (B 

AQ=Q (B 

1=1 ( B 

IF ( ALPHA. GT. -6.) 1=2 (B 

IF ( ALPHA. GT.-l.) 1=3 (B 

IF ( ALPHA. GT. 4. ) 1=4 (B 

IF { ALPHA. GT. 9. ) 1=5 (B 

IF ( ALPHA. GT. 15. ) 1=6 <B 

L I NE = L l NE + 3 <B 

IF (LINE. LE. 40) GO TO 75 (B 

L I NE = 0 (B 

WRITE (6,110) (B 

SUBTRACT OUT WEIGHT TARES 

75 BALDR AG=B ALDRAG-BDT ( I ) (B 

BALLIFT=BALLIFT-BLT( I ) (B 

SCLL I FT = SCLL I F T-SLT ( I ) (B 

DIVISION OF FORCES BETWEEN WING, FANS, AND TAIL 

IF ( I TYPE-2 ) 1,1,5 <B 

5 TL=SCLLIFT-BALLIFT*COS(ALPHA*RAD)+ BALDRAG*SIN( ALPHA*RAD) (B 

TD = SCLDRAG-BALLIFT*‘SIN< ALPHA*RAD>- BALDRAG*COS ( ALPHA*RAD) ( B 

FL = SCLLIFT-SLAT( ITUN, I)*Q*SW ( B 

FD=SCL DRAG-SOAK ITUN, I )*Q*$W ( B 


110 ) 

111 ) 

112 ) 

113) 

114) 

115) 

116) 

117) 

118) 

119) 

120 ) 
121 ) 
122 ) 

123) 

124) 

125) 

126) 

127) 

128) 

129) 

130) 

131) 


132) 

133) 

134) 


135) 

136) 

137) 

138) 

139) 

140) 
141 ) 

142) 

143) 

144) 

145) 

146) 


147) 

148) 

149) 


150) 

151) 

152) 

153) 

154) 
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CN=BALLIFT/(Q*ST) 

(B 

155) 

CA=BALDRAG/(Q«ST) 

( B 

156) 

IF ( I TR I P • EQ * 2 ) GO TO 42 

(B 

157) 


WRITE OUT UNCORRECTED DATA 


WRITE 16,1111 JTUNI I TUN) , I RUN, I FRAME , IRPM , AL PHA , Q , SCI L I F T , SCLORAG, 

(B 

158) 

1 CNt CA, TL,TD,FL,FD 


( B 

159) 

IF ( IRPM.EQ.O ) GO TO 1000 


(B 

160) 

IF (IRPM.EQ.8 ) T ST AT 10 = 40 • 2 


(B 

161) 

IF (IRPM. EQ. 10) TST AT IC =66. 3 


(B 

162) 

IF (IRPM. EQ. 12) TSTATIO95.0 


( B 

163) 

QJ=TSTATIC/(2.*SF) 


( B 

164) 

IF (Q.LE.O.) Q= 10E-10 


( B 

165) 

QOQJ = Q/ Q J 


(B 

166) 

VOVJ=SQRT(QOQJ) 


(B 

167) 

VO V JUN= VQ V J 


(B 

168) 

T0TS=( SCLL I FT/T STATIC) ♦(0.002378 /DENSITY) 


(B 

169) 

DOL=(SCLDRAG-SDAT( I TUN , 3 ) «Q^ SW ) / SCL L I FT 


(B 

170) 

WRITE (6,1001) QOQJ, VOVJ, TOTS, DCL 


( B 

171 ) 

RPM=10)0.*FLOAT( IRPM) 


(B 

172) 

CL=SCLL IFT / ( Q*SW ) 


(B 

173) 

I S I Z E — 3 


( B 

174) 

IF (Q. IT. 0.8) Q=CL=CN-VOVJ=Q0QJ=0. 


(B 

175) 

PUNCH CARDS FOR SUBSEQUENT PLOTTING OF UNCORRECTED 

DATA. 



PUNCH 2023, IRUN, IFRAME , RPM, ALPHA, Q, SCLL IFT, CL ,CN, 

QOQJ, VOVJ, DOL, 

(B 

176) 

1 TOTS , ITUN, I SI ZE 


( B 

177) 

CONTINUE 


( B 

178) 

ELIMINATE STATIC THRUST POINTS FROM THE DATA TO BE 

CORRECTED. 



IF (Q.LE.O. 8) GO TO 1 


( B 

179) 

CORRECT FOR SOLID BLOCKAGE 




QFAC=( 1«0+EPSIL0N( ITUN) )**2 


(B 

180) 

Q=AQ*OFAC 


(B 

181) 

CL W = SL A T ( ITUN, I )/QFAC 


(B 

182) 

CD W= SDAT ( ITUN, I )/QFAC 


(B 

183) 

C 00= SDA T ( I TUN , 3 ) /QF AC 


( B 

184) 

CN=CN/ QFAC 


(B 

185) 

C A = C A/QFAC 


(B 

186) 

SOLVE MOMENTUM QUARTIC FOR SKEW ANGLE AND VELOCITY 

RATIO. 



XT^O. 


(B 

187) 

DE LTX = 0 .01 


(B 

188) 

AX=— 1 . 


( B 

189) 

IF (CL) 1*44,45 


(B 

190) 

DOL=FD/FL 


(B 

191) 

VWH=-SQRT (2.^Q*SF/FL) 


(B 

192) 

i X= XT +DE L T X 


(B 

193) 

XT = X 


(B 

194) 

IF (XT. GT. 1.02) GO TO 51 


(B 

195) 

IF (XT. LT. -0.01 ) GO TO 52 


(B 

196) 

X=( 1.0+DOL*DOl)*X*X*X*X«-2.0*DOL*VWH*X*X*X + VWH*VWH*X*X-1.0 

(B 

197) 

IF ( ABS (X ) .LT .0.000001 I GO TO 47 


(B 

198) 

IF IX/AX) 48,47,49 


< B 

199) 


45 
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48 

DELTX=-0.5*DELTX 

(B 

200) 

49 

AX = X 

(B 

201) 


GO TO 46 

(B 

202) 

47 

WWH=XT 

( B 

203) 


VWO=VWH/WWH 

( B 

204) 


T ANCH I =-VWO-DOL 

(B 

205) 


CHIM-ATANITANCHI ) /RAD 

(B 

206) 


CHIEFF=45.0*CHIM/2. 

(B 

207) 


GO TO 53 

(B 

208) 

44 

CHIM=CHIEFF=90. 

(B 

209) 


VW0=10E 10 

(B 

210) 


GO TO 53 

(B 

211) 

51 

WRITE (6,151) IRUN, I FRAME 

(B 

212) 


GO TO 1 

(B 

213) 

52 

WRITE (6,152) IRUN, I FRAME 

( B 

214) 


GO TO 1 

(B 

215) 


SELECT TABLES OF INTERFERENCE FACTORS FOR APPROPRIATE TUNNEL. 


53 DO 60 I CH I = 1 , 8 

C WLFP { I CHI I = WLF P ( I TUN, I , ICHI > 

CULFP ( ICHI )=ULFP ( I TUN, I, ICHI ) 

CWDFP( ICHI ) =WDFP( I TUN, I , ICHI) 

CUDFP ( ICHI ) = UDFP ( ITUN, I, ICHI ) 

IF (ITRIP.EQ.il GO TO 60 
CWLF F ( ICHI )=WLFF ( ITUN, I, ICHI) 

CULFF ( ICHI ) =ULFF ( ITUN, I , ICHI) 

CWDFF ( ICHI )=W0FF ( ITUN, I, ICHI ) 

CUOFFI ICHI ) =UDFF ( ITUN, I , ICHI ) 

CWLFWI ICHI ) = WLFW ( ITUN, I, ICHI ) 

CULFWl ICHI )=ULFW( ITUN, I , ICHI ) 

C WDF W ( ICHI )- WDFW ( ITUN, I, ICHI ) 

CUDFW ( ICHI ) =UDF W I ITUN, I , ICHI ) 

CWLFTI ICHI )=WLFT( ITUN, I, ICHI ) 

CIJLFTI ICHI ) -ULFT ( ITUN, I , ICHI ) 

CWDFTI ICHI ) =WDFT ( ITUN, I, ICHI) 

60 CUDFTI ICHI )=IJ0FT( ITUN, I , ICHI ) 

C 

C INTERPOLATE FOR CHI EFF IN TABLES OF INTERFERENCE FACTORS 

C 

CALL D I SCOT ( CHI EFF , CH I EFF , CH I , C WL F P , C WLFP , -0 30 1 8 , 0 , OWL FP ) 
CALL D I SCOT ( CH I EFF, CH I EF F, CH I , CUL FP , CULFP , -030, 8 , 0 , OUL FP ) 
CALL 01 SCOT (CHIEFF,CHIEFF,CHI,CWOFP,CWDFP,-O30, 8,0,DWDFP) 
CALL DI SCOT ( CHI EFF, CH I EFF, CH I , CUDFP , CUDFP ,-0 30, 8, 0, DUDFP ) 
IF (ITRIP.EQ.l) GO TO 39 

CALL DISCOT ( CH I EFF, CH I E FF , CH I , CWLFF, CWLFF , -030, 8, 0, DWLFF ) 
CALL DISCOT (CHIEFF, CHI EFF, CHI, CULFF, CULFF, -030, 8,0, DULFF) 
CALL DISCOT ( CH I EFF , CH I E FF , CH I , CWDFF , CWDF F , -0 30, 8, 0 , D WDF F ) 
CALL DI SCOT ( CH I E FF , CH I E FF ,CH I , CUOFF , CUDF F , -0 30 , 8 , 0 , DUDF F ) 
CALL DISCOT {CHIEFF, CHIEFF, CHI, CWLFW,CWLFW, -030,8,0,0 WLFW) 
CALL DISCOT ( CH l E FF , CH I E FF , CH I , CULF W , CULFW , -0 30 , 8 , 0 , DULFW ) 
CALL DISCOT ( CHI EFF, CHIEFF ,CH I ,CWDFW, CWDFW ,-03C, 8, 0,DWDFW ) 
CALL DI SCOT ( CH I EFF , CH I E FF , CH I ,CUDF W, CUDFW , -030 , 6 , 0 , DUOFW ) 
CALL DISCOT ( CH I EFF , CH I E FF , CH I , C WLF T , C WLF T ,-0 30 , 8 , 0 , D WLF T ) 
CALL DISCOT { CHI EFF ,CHI EFF ,CHI »CULFT, CULFT , -0 30, 8,0 , DULFT ) 
CALL DISCOT ( CH I E FF , CH I E FF , CH I , C WDFT , CWDFT ,-030, 8 , 0 , D WDF T ) 
CALL 01 SCOT (CHIEFF, CHIEFF, CHI, CUDFT,CUDFT, -030, 8, 0,DUDFT) 


(B 216) 
(B 217) 
(B 218) 
(B 219) 
(B 220) 
(B 221) 
(B 222) 
(B 223) 
(B 224) 
(B 225) 
(B 226) 
(B 227) 
(B 228) 
(B 229) 
(B 2301 
(B 231) 
(B 232) 
(B 233) 


(B 234) 
(B 235) 
(B 236) 
(B 237) 
(B 238) 
(B 239) 
(B 240) 
(B 241) 
(B 242) 
(B 243) 
(B 244) 
(B 245) 
(B 246) 
(B 247) 
(B 248) 
(B 249) 
(B 250) 
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C 

C CORRECT TUNNEL PITOT MEASUREMENT FOR WALL EFFECTS 

C 

3 9 OWP=2.*DWLFP*SF/( AT ( ITUN ) *VWO ) +2 . *OWOf P*S F*DOL / ( AT{ I TUN J *VWO) 
1 +WLWP( I TUN f I )*CLW*(-0.25)*$W/AT<ITUN) 

DUP=2.*DULFP*SF/< AT t I TUN ) *VWO ) +2 . *OUDFP*SF*DOL / ( A T (I TUN ) *VWO I 
1 +ULWPI I TUN f I )*CLW*(-0.25)*SW/AT(ITUN) 

DWP=DWP*WOVPIT( ITUN) 

DUP = DUP + UO VP I T ( ITUN) 

QFACT=DWP*DWP*( 1 • +DUP ) * ( 1 • +DUP ) 

IF (ITRIP.EQ.l) QF ACTOR-QFACT 
0 = AQ /OF ACT 

IF ( ITRIP.EQ.2) GO TO 40 
I TR I P=2 
GO TO 5 

40 VWO = VWO*$QRT(QFACTOR/QFACT ) 

CL W-CLW+QFACT/QF ACTOR 
CDO=CDO*QFACT/QFACTOR 
CN = CN *0 FACT /OF ACTOR 
C 

C FIND INTERFERENCE VELOCITY RATIOS 

C 

DWF=2«*DWLFF*SF/(AT( I T UN ) * VWO ) +2. *DWDFF*S F*DOL / ( AT (ITUN) * VWO ) 
1 +WLWF( ITUNt I )*CLW*SW*(-0.25)/AT ( ITUN) 

DUF=2»*DULFF*SF/(AT( ITUN) * VWO )+2.*DUDFF*SF*D0L/(AT(ITUN)*VW0) 
I + ULWF ( ITUNt I )*CLW*SW*(-0.25)/ATIITUN) 

DWW = 2.*DWLFW*SF/< AT( ITUN )* VWO ) +2 . *D WOFW*S F*DOL / ( AT (I T UN ) * VWO ) 
1 +WL WW ( ITUNt I ) *CLW*SW*< -0.25)/ AT ( ITUN) 

DUW=2.*DULFW*SF/(AT( ITUN ) *VWO ) + 2* *DUDFW*$F*DOL / ( AT ( ITUN)*VWO) 
1 + UL WW ( I TUN 1 1 )*CLW*SW*l-0.25)/AT(ITUN) 

DWT=2-*DWLFT*SF/(AT( ITUN )♦ VWO ) + 2 . *DWDFT*SF*DOL / ( A T (I TUN ) * VWO ) 

1 *WLWT( ITUNtI )*CLW*SW*(-0#25)/AT(ITUN) 

2 +WLTT ( ITUN, I ) *CN*ST*< -0. 25) / AT < I TUN ) 

DUT-2**DULFT*SF/<AT( I TUN ) *VWO )«-2 . *OUCFT*SF*OOL / ( AT ( I TUN ) *VWO ) 

1 +ULWT( ITUNt I )*CLW*SW*(-0.25)/AT(ITUN) 

2 +ULTT ( ITUN, I ) *CN*ST* I -0 . 25 ) / AT ( I TUN ) 


(B 251) 
IB 252) 
(B 253) 
(B 254) 
(B 255) 
(B 256) 
(B 257) 
(B 258) 
(B 259) 
(B 260) 
(B 261) 
(B 262) 
(B 263) 
(B 264) 
(B 265) 
(B 266) 


(B 267) 
(B 268) 
(B 269) 
(B 270) 
(B 271) 
(B 272) 
(B 273) 
(B 274) 
(B 275) 
(B 276) 
IB 277) 
(B 278) 
(B 279) 
(B 280) 


CALCULATE CORRECTIONS 


TANALPF*DWF/f l.*DUFI 
T AN AL P W= DW W/ ( l.+DUW) 

TANALPT = DWT/(1.*DUT ) 

DALPF=ATAN(TANALPF)/RAD 

OALPW=ATAN(TANALPW)/RAD 

DALPT=ATAN(TANALPT)/RAD 

QOQF = DWF*DWFM l.+DUF 1**2 

QOQW-DWW*DWW*( l.+DUW)**2 

QOOT=DWT* DWT M1.*DUT)**2 

AL PCF= AL PHA + OALPF 

ALPCW=ALPHA*DALPW 

ALPTC=ALPHA+DALPT 

OCF=Q*OOOF 

QC W=Q*QOQW 

QCT=Q*QOOT 

SLIFTC=SCLL IFT*COS( DALPW *R AO )- ( SCLDRAG-CD 0*Q*SW ) *S I N( OALPW*R AD) 
SDR AGC= SCLL I FT*S I Nl DAL PW*R AD ) «• I SCLQRAG-CDO*Q*S W ) *COS ( DALPW*R AD ) 
CNTC=CN/QOQT 


(B 281) 
(B 282) 
(B 283) 
(B 284) 
(B 285) 
(B 286) 
(B 287) 
(B 288) 
(B 289) 
(B 290) 
(B 291) 
(B 292) 
(B 293) 
IB 294) 
(B 295) 
CB 296) 
(B 297) 
IB 298) 
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C 

C ADJUST FOR DIFFERENCES IN CORRECTIONS AT WING, FANS, AND TAIL 

C ASSUME A TAIL EFFICIENCY FACTOR t QT/Q 1 OF 0.9 

C 

DDALPT=DALPT-DALPW 

CNTCA=(CNTC-0.030*D0ALPT*0.9) 

TOTL=SLI FTC -CN*ST*Q *COS ( AL PH A*RA D ) 

1 *CNTCA*QCW*ST*COS( AL PH A*RAD ) 

CDIUN=CN*CN/ (PI*2.40) 

CDIC=CNTCA*CNTCA/< PI *2.40) 

SORAGC=SDR AGC-CD IUN*Q*ST +0010*00 W*ST 

DDAFW=DALPF-DALPW 

TST AT I C =40 • 2 

IF (IRPM.EQ.10) TSTAT I C=66. 3 
IF (IRPM. EQ. 12) TST AT IC = 9 5.0 
QJ=TSTATIC/( 2.*SF ) 

QOQ J -OC W/ 0 J 

OOCJF=QCF/OJ 

VOVJ=SORT(QOQJ) 

VOVJF*SQRT I QOQJF ) 

DVOVJ=VOVJF-VOVJ 
SA=S IN( ALPCW+RAD ) 

C A = COS ( ALPCW*RAD ) 

TOTl=TOTL+TSTATIC*DDAFW*RAO*SA 
$DRAGC=SDRAGC-TSTATIC*(DVOVJ+DDAFW*RAO*CA ) 

TOTDOL-SDRAGC/TOTL 

C 

C STORE CORRECTED VALUES CN TAPE 1 ANC CORRECTIONS ON TAPE 4 

C 

WRITE (I) I TUN ,IRUN, IFRAME, IRPM, ALPHA, Q , CH l M , DALPF , QOQF , 

1 0 ALPW,QOOW,DALPT,OOQT,CCW,AL PC W,CNTCA,TOTL,TOTDOL, DENSITY 

WRITE (4) I TUN, I RUN , A L PHA , AL PC W , D ALP W , OALP F , DALPT , QOQW , QOQF , 

1 OOQT, VOVJ,DVOVJ , VOVJUN 

GO TO 1 
998 ENCFILE I 
RFWIND I 
ENDFILE 4 
REWIND 4 
WRITE (6,125) 

C 

C WRITE OUT CORRECTED VALUES FROM TAPE 1 

C 

L I NE =0 

DO 61 K = 1 , 1000 

READ (1) ITUN ,IRUN, IFRAME, IRPM, ALPHA, 0 * CH I M , D ALP F , OOQF , 

1 DALPW,QOQW,DALPT,QOQT,QCW,ALPCW,CNTCA,TQTL,TOTDOL* DENSITY 

IF ( EOF , 1 ) 999,62 

62 WRITE (6,126) J TUN ( I TUN ) , I RUN , I FR AME , I RPM , A LPH A , Q , CH I M , DALPF , 
1 QOQF , DAL PW,QOOW, DALPT, OOQT, QCW, ALPCW, CNTC A, TOTL 


(B 299) 
(B 300) 
(B 301) 
(B 302) 
(B 303) 
(B 304) 
IB 305) 
( B 306) 
(B 307) 
(B 308) 
(B 309) 
(B 310) 
(B 311) 
(B 312) 
(B 313) 
(B 314) 
(B 315) 
(B 316) 
(B 317) 
(B 318) 
(B 319) 
(B 320) 


(B 321) 
(B 322) 
(B 323) 
(B 324) 
(B 325) 
(B 326) 
(B 327) 
(B 328) 
(B 329) 
(B 330) 


(B 331) 
(B 332) 
IB 333) 
(B 334) 
(B 335) 
(B 336) 
(B 337) 


C 

C 

c 


NON DIMENSION ALIZE CORRECTED VALUES. 


IF ( IRPM.EQ.O ) 

GO TO 1002 

(B 

338) 

IF ( IRPM. EQ. 8 ) 

TST AT IC *40 • 2 

(B 

339) 

IF ( IRPM.EQ.10) 

TST AT I C *66 • 3 

(B 

340) 

IF (IRPM. EQ. 12) 

TSTAT IC *9 5. 0 

(B 

341) 

QJ=TSTATIC/(2.*SF» 

(B 

342) 

QOCJ=QCW/QJ 


(B 

343) 

VnVJ=SQRT(QOQJ» 


(B 

344) 

TOTS= (TOTL /TSTATIOM 0.002 378/DENS I TY) 

(B 

345) 
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WRITE (6,10011 Q0QJ,V0VJ,T0TS, TOT DO L (B 346) 

STORE CORRECTED N0N0 1 MENS I0NAL VALUES ON TAPE 3 


WRITE (3) I TUN, l RUN, I RPM , ALPHA , Q , DALP W, QOQW , QC W , ALPCW , CN TC A , TOTL , 
1 OOQJ, VOVJ,TOTS,TOTDOL 
L I NE=LI NE+3 
1002 CONTINUE 

IF (LINE. LT. 40) GO TO 61 
L I NE = Q 

WRITE (6,125) 

61 CONTINUE 
999 ENOFILE 3 
REWIND 3 


(B 347) 
(B 348) 
(B 349) 
(B 350) 
(B 351) 
(B 352) 
(B 3 53) 
(B 354) 
(B 355) 
(B 356) 


INTERPOLATE CATA ON TAPE 3 
CORRECTEC ANGLES OF ATTACK 


TO OBTAIN CORRECTEO VALUES AT FIXED 


WRITE (6,2020) 

L I NE=0 

2003 DO 2001 1=1,6 

READ (3) I TUN, I RUN, I RPM, T ALP (I ) , TQ (I ) ,T0 ALPW ( I ) , TQOQW ( I ) ,TQC( I ) , 

1 T ALPC ( I ) , T CNC ( I ) , TL I FT ( I ) , T QOQ J ( I ) , TVOVJ( I), TLOTSC I ) ,TDOL( I ) 

IF (EOF, 3) 9999,2001 
2001 CONTINUE 

00 2002 1=1,5 

IF ( I.EQ.1.AND.TALPC( t I.GT.-4.5) GO TO 2002 
IF ( I.EQ.1.AND.TALPC( 1I.GT.-4.5) L INE = L I NE* 1 
AA= A ( I ) 

CALL DI SCOT ( AA , A A , T AL PC , T AL P , T ALP * -0 10 , 6 , 0, AUN ) 

CALL 0 I SCOT (AA, AA, TALPC ,TG, TQ, -010 ,6,0, QUN ) 

CALL DISCOT ( A A , A A , TALPC , TDAL P W, TDALPW ,-0 1 0 , 6 , 0, 0 AL W ) 

CALL DISCOT ( A A , A A , TALPC , TQOQW , T COQW , -010 , 6 , 0, QOQ ) 

CALL D I SCOT ( AA , AA, T AL PC , T QC , TQC ,-0 10 ,6, 0 , QC ) 

CALL DISCOT (AA, AA, TAL PC ,TCNC,TCNC, -010,6,0, CNC) 

CALL DISCOT ( A A , A A , T ALPC , TL I F T , TL I F T , -010 , 6 , 0 , TLFT ) 

CALL DISCOT ( A A, A A, T ALPC , TQOQ J , TQOQ J ♦ -01 0 , 6 , 0 , QOQ J ) 

CALL D I SCOT ( A A , A A , T ALPC , T VOV J , T VOV J , -010 , 6 , 0, VCV J ) 

CALL D I SCOT ( A A , A A , T ALPC , TLOT S , TLOT S , -0 10 , 6 , 0, LOT S ) 

CALL DI SCOT ( AA , AA, TALPC , TOOL , TOOL , -010, 6, 0, CL ) 

DOTS=DL*LOTS 


(B 357) 
(B 358) 
(B 359) 
(B 360) 
(B 361) 
(B 362) 
(B 363) 
(B 364) 
(B 365) 
(B 366) 
(B 367) 
(B 368) 
(B 369) 
(B 370) 
(B 371) 
(B 372) 
(B 373) 
(B 374) 
(B 375) 
(B 376) 
(B 377) 
(6 378) 
(B 379) 


WRITE OUT, AND PUNCH FOR SUBSEQUENT PLOTTING, THE CORRECTED 
NONDI MENS IONAL VALUES. 


WRITE (6,2021) 

JTUN( I TUN ) , l RUN, I RPM, AA,QC, QOQ J, VOV J, TLFT, LOTS, 

(B 

380) 

1 CNCtOL.OOTSt AUN,QUN, OALW.QOO 

(B 

381) 

RPM=1000.*FLOAT< IRPM) 

( B 

382) 

CL=TLFT / ( QC*SW ) 


(B 

383) 

I S I Z E=1 


(B 

384) 

PUNCH 2023, IRUN, I RUN , RPM, A A, OC, TLFT, CL , CNC , QOQ J , VOV J , DL , LOTS, 

(B 

385) 

1 I TUN , ISIZE 


(B 

386) 

LINE = LINEU 


( B 

387) 

IF (LINE. LT. 35) 

GO TO 2002 

(B 

388) 

L I NE=0 


( B 

389) 

WRITE (6,2020) 


( B 

390) 

2002 CONTINUE 


( B 

391) 

WRITE (6,2022) 


( B 

392) 

GO TO 2003 


(B 

393) 

9999 WRITE (6,3000) 


( B 

394) 
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WRITE (6,3008) 
I C0R = 1 
L I NE = 0 


(B 395) 
(B 396) 
(B 397) 


INTERPOLATE DATA ON TAPE 4 TO OBTAIN VALUES OF CORRECTION ANGLES 
AND VELOCITY RATIOS AT FIXED CORRECTEC ANGLES OF ATTACK. 


3003 DO 3001 1=1,6 

READ (4) I TUN, IRUN,TALP< I ) , T AL PC ( I ) , T 0 AL PW < I ) ,TDALPF( I ) , TDAL PT ( I ) , 
1 T COQW ( l ) , TOOQF ( I) , TOOQT (I ),TVOVJ( I ) , T DVO V J ( I ) , VOV JUN 
IF (EOF, 4) 3002,3001 

3001 CONTINUE 

DO 3004 1 = 1,5 

A A = A ( I ) 

CALL D I SCOT ( AA, A A, TALPC , T ALP , TALP ,-010 , 6, 0, AUN > 

CALL D I SCOT < AA, AA, TAL PC ♦ T DAL PW , TDALP W , -0 10, 6 , 0 , DAL W ) 

CALL 01 SCOT (AA,AA,TALPC,TDALPF,TDALPF,-010,6,0,DALF) 

CALL 0 I SCOT <AA,AA, TALPC, T OALPT , T DALPT , -010, 6,0, DALT) 

CALL 0 1 SCOT ( AA,AA, TALPC, TQOQW,TQOQW, -010, 6,0, TOW) 

CALL DISCOT ( AA , AA, TALPC , TQOQF , TQOQF , -0 10, 6, 0, TQF ) 

CALL DISCOT ( AA, AA, TALPC ♦ TQOQT,TQOQT ,-010 ,6,0, TQT ) 

CALL D I SCOT ( A A , A A , T ALPC , T VO V J , T VOV J , -0 10 , 6 , 0, VOVJ ) 

CALL DISCOT C AA , A A , T ALPC , T 0 VO V J , TD VOV J ,-01 0 , 6, 0 , DVOV J ) 

DI T=DALT-DALW 
D I F=DAL F-DAL W 
QTOQC= TQT/TQW 
QFOQC=T OF/ TQW 
L I NE = L I NE + 1 


(B 398) 
(B 399) 
(B 400) 
(B 401) 
(B 402) 
(B 403) 
(B 404) 
(B 405) 
(B 406) 
(B 407) 
(B 408) 
(B 409) 
(B 410) 
(B 411) 
(B 412) 
(B 413) 
(B 414) 
( B 415) 
(B 416) 
(B 417) 
(B 418) 


WRITE OUT, AND PUNCH FOR SUBSEQUENT PLOTTING, THE CORRECTION 
ANGLES AND VELOCITY RATIOS. 


WRITE (6,300 5) JTUN( I TUN) , I RUN , A A , AUN , VO V J ,DALW , D I T , DI F , TQW, 
1 QTOQC , QFOQC , D VO V J 

PUNCH 3006, I RUN, AA, VOV J , DALW , D I T , D I F , TQW , QTOQC , QFOQC , DVOV J , 
1 I TUN , I COR 
3004 CONTINUE 

WRITE (6,3007) 

L INF = L INE + 1 

IF (LINE. LT. 36) GO TO 3003 

GO TO 9999 
C 

C WRITE OUT, AND PUNCH FOR SUBSEQUENT PLOTTING, THE CORRECTION 

C ANGLES AND VELOCITY RATIOS IN THE UNCORRECTED DATA. 

C 

3002 REWIND 4 
l C 0R= 3 

3011 WRITE (6,3009) 

WRITE (6,3008) 

L I NE = 0 

3014 DO 3010 1 = 1,6 

READ (4) I TUN, I RUN, ALPHA , ALPC W , DALPW , D ALP F, DALPT, QOQW, QOQF , 

1 QOQT, VOVJ , DVOV J, VOV JUN 
IF (EOF, 4) 3012,3013 

3013 DIT=DALPT-DALPW 
OIF=DALPF -DALPW 
CTOQC=UOQT/QOQW 
QF OQC = QOQF /QOQW 


(B 419) 
(B 420) 
(B 421) 
(B 422) 
(B 423) 
(B 424) 
(B 425) 
(B 426) 
(B 427) 


(B 428) 
(B 429) 
(B 430) 
(B 431) 
(B 432) 
(B 433) 
(B 434) 
(B 435) 
(B 436) 
(B 437) 
(B 438) 
(B 439) 
(B 440) 
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LINE=LINE*1 

WRITE (6,3005) JTUN( ITUN I , l RUN, ALPCW, ALPHA , VO VJ , DALP W, D I T , D I F , 
1 QOQW,QTOQC,QFOQC,DVOVJ 

PUNCH 3006, IRUN, ALPHA, VO V JUN , D AL PW , 0 I T , DIF, QOQW, QTOQC , QFOQC , 

1 DVOVJ,ITUN,ICOR 
3010 CONTINUE 

WRITE (6,3007) 

LINE^LINEU 

IF (LINE. LT. 35) GO TO 3014 
GO TO 3011 
3012 STOP 


( 8 441 ) 
(B 442) 
(B 443) 
(B 444) 
(B 445) 
(B 446) 
(B 447) 
( B 448) 
(B 449) 
(B 450) 
(B 451) 


FORMATS 


100 FORMAT (2I2,2I1,I4,F5.2,5F6.2,F7.1,F8.6) 

110 FORMAT ( 1 HI / 37X*FAN- I N-W I NG TESTS IN AMES 7X10 TUNNEL (AARL TEST 

1 NO. 6) *//4X*TUNNEL RUN FRAME RPM*5X*ALPHA*<,X*Q*6X*L*8X*D* 

2 5X*CN(T) CA(T)*4X*L(WB)*4X*D(WB)*5X*L(F1*5X*D(F)*2X/) 

111 FORMAT (2X,A10»2X,I2, 16, 14*, 000*F7 . 1 , F7. 2 , 2F9. 3, 2F 8. 4, 4F9. 3 , F8 . 2 ) 

120 FORMAT ( 1 HI / /30X * I NTERFERE NCE FACTORS*) 

121 FORMAT (8F7.4,16X,A8) 

122 FORMAT (49X,F7.4, 16X»A8) 

123 FORMAT ( 8F 10 . 4 , 5X , A 10 ) 

124 FORMAT ( 70X , F10 . 4 , 5 X , A 10 ) 

125 FORMAT ( 1H1//50X*C0RRECTED OAT A* //4X* TUNNEL RUN FRAME RPM*5X 

1 *ALPHA*4X*Q*3X*CHI DAL PF OOOF DALPW OOOW OALPT QCQT* 

2 4X*0CW ALPC(W)*4X*CNTC TOT L*/» 

126 FORMAT (2X.A10, 14, 1 6 , 1 4* , 000*F 7 . 1 , 2F 7 . 2 , 3 ( F 7 . 2 , F 7 . 4 ) 

1 ,F7.2,F8.2,F9.4,F8.2 ) 

127 FORMAT (//35X* A10// ) 

151 FORMAT ( 1 OX*R UN *12*, FRAME *14* FAILS, 

152 FORMAT (10X*RUN *12*, FRAME *14* FAILS, 

153 FORMAT (20X*CHI =*F6.2*, CHI E =*F6.2*, 

1001 FORMAT (10X*0/0J = *F8.4*, V/VJ = *F8.4 

1 F 6 . 3 / ) 

2020 FORMAT ( 1H1//49X* INTERPOLATED CORRECTED VALUES*//12X*TUNNEL*5X*RUN 
1*2X*RPM*2X*ALPHA C* 3X*QC *3 X*Q/ 0 J*2 X* V/ V J*2X*L I FT *4X*L/T S *3X*CNTC* 
23X*D/L*3X*D/TS*2X*ALPHA*4X*Q*4X*D ALP*2X*QC/0*/ ) 

2 021 FORMAT ( 1 OX , A10 , I 6 , 1 3*000*2F7 . 2 , 2F6.4 , F 7. 2 , 2F7 . 3 , F6 . 3 , F7 . 3 , F 6 . 2 , 

1 2F7.2,F8.4) 

2022 FORMAT ( ) 

2023 FORMAT ( I 3 , 1 4, F 5 . 0 , 3Ffc . 2 , 2F7 .4, 4F6 . 4, 2 1 1 ) 

2024 FORMAT (A8) 

2025 FORMAT (/30X.A8* WING LOADING*) 

3000 FORMAT ( 1H1//43X*C0RRECTICNS ACCORDING TO CORRECTED ALPHAS*/) 

3005 FORMAT (4X, A10, 17, 2F10 .2 ,F 10 .4 , 3F10 .2, 4F 10 . 4 ) 

3006 FORMAT ( I 3 , 9F 8. 4 , 2 1 1 ) 

3007 FORMAT () 

3008 FORMAT ( 6X*TUNNEL*6X*RUN*4X*ALPHAC*4X* ALPHAU*6X*V/V J*4X*0 ALPW* 

1 6X*D I T*6X*D IF*5X*QC/QW*5X*QT/QC*5X*QF/QC*3X*D VF/VJ*/) 

3009 FORMAT ( 1H1//42X*C0RRECTI0NS ACCORDING TO UNCORRECTED ALPHAS*/) 

END 


WWH TOO GREAT*) 
WWH TOO SMALL*) 
V/WO =*F8.2*» D/L 
*, L/TS = *F8 .4*, 


=*F8 • 4 ) 
D/L = * 


(6 

452) 

l B 

453) 

( B 

454) 

( B 

455) 

(B 

456) 

( B 

457) 

( B 

458) 

( B 

459) 

( B 

460) 

< B 

461) 

t B 

462) 

( B 

463) 

(B 

464) 

( B 

465) 

(B 

466) 

( B 

467) 

( B 

468) 

( B 

469) 

(B 

470) 

(B 

471) 

(B 

472) 

(B 

473) 

l B 

474) 

(B 

475) 

(B 

476) 

( B 

477) 

( B 

478) 

(B 

479) 

< 8 

480) 

(B 

481) 

(B 

482) 

(B 

483) 

(B 

484) 

<B 

485) 

<B 

486) 

(B 

487) 

(B 

488) 

(B 

489) 
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SUBROUTINE DISCOT ( X A , Z A , T ABX , TA BY , TA8Z ♦ NC , N Y , NZ , ANS ) 

(B 

490) 


DIMENSION T ABX (21 ,TABY(2) ,TABZ(2 ) 

( B 

491) 


DIMENSION NPX(8),NPY(8),YY(8) 

( B 

492) 


CALL UNS (NCflA.IDXtlOZ, IMSI 

(B 

493) 


IF C NZ- 1) 5,5,10 

( B 

494) 

5 

CALL DISSER ( X A , T ABX < 1 1 , 1 , NY , IDX , NN * 

( B 

495) 


NNN = I DX + 1 

( B 

496) 


CALL LAGRAN { X A , T ABX { NN ) , T AB Y ( NN ) , NNN , ANS ) 

( B 

497) 


GOTO 70 

(B 

498) 

10 

Z A P. G = Z A 

( B 

499) 


l P 1 X= l DX+ 1 

( B 

530) 


IP1Z=10Z+1 

( B 

501) 


IF (IA) 15,25,15 

(B 

502) 

15 

IF (ZARG-TABZ(NZ) > 25,25,20 

( B 

503) 

20 

7 ARG=T A 0Z ( NZ 1 

( B 

504) 

25 

CALL DISSER ( Z ARG , TABZ ( U , 1 , NZ, I DZ , NP Z > 

( B 

505) 


NX-NY/NZ 

(B 

506) 


NPZL=NPZ+IDZ 

( B 

507) 


1 = 1 

(B 

508) 


IF (IMS) 30,30,40 

(B 

509) 

30 

CALL DISSER ( X A, T A BX ( 1 ) , 1 , NX , I DX , NPX ( 1 ) ) 

(B 

510) 


DO 35 J J = NPZ , NPZL 

(B 

511) 


NPY ( I )= ( JJ-1 )*NX+NPX{ 1 ) 

( B 

512) 


NP X ( I ) =NPX ( 1 1 

(B 

513) 

35 

1 = 1+1 

( B 

514) 


GOTO 50 

(B 

515) 

40 

DO 45 JJ=NPZ,NPZL 

( B 

516) 


I S=( JJ-1)*NX+1 

(B 

517) 


CALL DISSFR (XA,TABX( 1) , IS, NX, IDX,NPX( I ) ) 

( B 

518) 


NP Y ( I ) = NPX ( I I 

( B 

519) 

45 

1 = 1+1 

( B 

520) 

50 

DO 55 LL=1, IP1Z 

( B 

521) 


NL OC =NPX (LL ) 

(B 

522) 


NL OC Y=NPY ( LL ) 

(B 

523) 

55 

CALL LAGRAN(XA,TABX(NLOC ) , T A BY ( NLOC Y ) , I P 1 X, Y Y ( LL ) > 

( B 

524) 


CALL LAGRAN ( Z ARG, T ABZ ( NPZ ) , YY ( 1 ) , I PI Z , ANS ) 

( B 

525) 

70 

RETURN 

(B 

526) 


FND 

( B 

527) 



SUBROUTINE UNS ( I C , I A , I DX , IDZ , I MS ) 

( B 

528) 


IF (ICI 5,5,10 

(B 

529) 

5 

I MS = 1 

(B 

530) 


NC = - 1 C 

(B 

531) 


GOTO 15 

( B 

532) 

10 

I M S= 0 

(B 

533) 


NC= I C 

( B 

534) 

15 

IF ( NC- 100 ) 20,25,25 

( B 

535) 

20 

I A = 0 

(B 

536) 


GOTO 30 

( B 

537) 

25 

I A = 1 

( B 

538) 


NC=NC- 100 

( B 

539) 

30 

I DX=NC / 10 

(B 

540) 


I0Z=NC-IDX*10 

(B 

541) 


RETURN 

(B 

542) 


END 

( B 

543) 
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SUBROUTINE DISSER ( XA , TA B , I , NX , I D , NPX ) 

{ B 

544) 


DIMENSION T AB ( 2 J 

(B 

545) 


NPT=ID+1 

(B 

546) 


NP B=NPT / 2 

(B 

547) 


NPU=NPT-NPB 

(B 

548) 


IE { NX-NPT ) 10,5,10 

(B 

549) 

5 

NP X= I 

IB 

550) 


RETURN 

(B 

551) 

10 

NL CW = I + NP B 

(B 

552) 


NUPP=I+NX-(NPU*1 ) 

(B 

553) 


DO 15 I I = NLOW , NUPP 

(B 

554) 


Nl 0C= I I 

(B 

555) 


IF (TABUI)-XAI 15 »20» 20 

(B 

556) 

15 

CONTINUE 

(B 

557) 


NPX=NUPP-NPB+1 

(B 

558) 


RETURN 

(B 

559) 

20 

NL=NLQC -NP B 

< B 

560) 


NU=NL+ID 

IB 

561) 


DO 25 J J= NL , NU 

(B 

562) 


ND I S = J J 

IB 

563) 


IF <TA8UJI-TAB< JJU1 > 25,30,25 

(B 

564) 

25 

CONTINUE 

(B 

565) 


NP X = NL 

IB 

566) 


RETURN 

( B 

567) 

30 

IF ( T AB ( ND I S ) — X A ) 40,35,35 

(B 

568) 

35 

NPX= ND I S — I D 

(B 

569) 


RETURN 

(B 

570) 

40 

NPX=NDI S+l 

( B 

571) 


RETURN 

I B 

572) 


END 

(B 

573 ) 



SUBROUTINE LAGRAN ( XA , X, Y , N, ANS ) 

I B 

574) 


DIMENSION X ( 2 ) , Y < 2 ) 

( B 

575 ) 


S U M = 0 • 0 

(B 

576) 


DO 3 1=1, N 

( B 

577 ) 


PROD= Y ( I ) 

( B 

578) 


DO 2 J= 1 , N 

l B 

579) 


A= X ( I ) - X I J ) 

( B 

580) 


IF (A) 1,2,1 

(B 

581 ) 

1 

B= ( XA- X ( J ) | /A 

t B 

582) 


PR0D=PR0D*6 

(B 

583) 

2 

CONTINUE 

( B 

584) 

3 

S UM = SUM + PR OD 

IB 

585) 


AN S= SUM 

(B 

586) 


RETURN 

I B 

587) 


END 

(B 

588) 
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. juOU 

.0116 
. 3681 
. 7516 
• 0556 


• 7751 
.4994 

• 5c>5 1 
.3149 
.4772 

• 2 0 o 7 

• 49 v 7 
.8694 


• 033& 
.026 

• C 9 1 0 
.8766 


.OuOO 

• 002 3 
.3865 
.7506 
.0537 


• 7o2 3 
.4752 
.^397 
.2954 

• 5010 
.2 785 

• 5 3o2 

• 9 3oC 


.0335 
.0403 
.C903 
. 6 7o4 


INTER! E RE NO t EACTOkS 
ELL IPT IC WING LCAUlNb 


2-1 RUUNC 
ALPHA = -10 


3 U . 00 00 

40.0000 

50.0000 

60.C00C 

70. OOuu 

80.0CO0 

9C.C000 

CHI 

-.8363 

-.6547 

-.5091 

- .4268 

-.4055 

-.4103 

-.4063 

-1C1WLFF 

.4843 

.4934 

.4245 

.3048 

.1654 

• 0271 

-.0776 

-101ULFF 

-.5460 

-.2678 

-.2319 

-.1354 

-.0619 

.0069 

.0776 

-1C1W0FE 

.1213 

.1216 

.0855 

.0416 

.0091 

-.CC46 

C.COOO 

-1C1UDEE 







-.6552 

-lOiwLwF 







-.0592 

-101ULWF 







-.5024 

-IOInLWW 







-.0383 

-1C1LLWW 

- .5925 

- • 444 7 

-.3449 

-.2951 

-.2643 

-.2873 

-.2821 

-lOlwLFW 

.5029 

• 44 9C 

.3578 

.2483 

. 1j5G 

.0260 

-.0508 

-101GLFW 

- .3865 

-.2553 

-.1644 

-.1011 

— . 051 u 

-.0027 

.0469 

-101WDPW 

. 3108 

.2749 

.2303 

. 1914 

. 1655 

.1542 

.1563 

-ICIUDFw 

- .6266 

-.7998 

-.9492 

-.9402 

-.8037 

-.8332 

-.9208 

-101WLFT 

-.1603 

-.0490 

.1622 

.421 7 

. 4299 

.1956 

- . C264 

-101ULPT 

-.5371 

-.5528 

-.5183 

-.39o7 

-.2341 

— • 1 C2 C 

.0437 

-lOlwGFT 

- • 81 67 

-.7219 

- .o06 6 

-.5216 

-. 5373 

-.5874 

-.6020 

-101UDFT 







-1.1440 

- lOlwL WT 







. 0035 

-1C1UL WT 







-.5791 

-1C1WLTT 







-.4450 

-1C1ULTT 

-.0327 

-.0322 

-.0324 

-.0329 

-.0335 

-.0335 

-.0329 

-lClwLFP 

-.0408 

-.0397 

-.0394 

-.0395 

-.0395 

-.0370 

-.0083 

-101ULFP 

-.0812 

-.0729 

- • Oc 5 7 

- • C59 1 

U52o 

- .0466 

-.0405 

-lQlwDFP 

.8740 

.87 15 

.8694 

.867 7 

• 666 8 

• 8666 

.8674 

-101UDFP 







-.0684 

-101 WL wP 







-.0408 

- 1 01 OL WP 







-.0493 

-lClwOWP 







1.1376 

-1C1U0WP 



ALPHA - 

-5 






30.0000 

40 . 00 CC 

50.0000 

60 .OOOC 

7C. OOUU 

80 .0000 

9C.CG00 

CHI 

-.8370 

-.6550 

-.5092 

- • 427 C 

4059 

- .4107 

-.4064 

-51WLFF 

.4835 

.4929 

.4238 

.3037 

. 1692 

.0258 

-.0788 

-51ULFF 

-.5452 

-.3667 

-.2307 

-.1342 

- • 0608 

.0081 

• C7d 8 

-51WDFF 

.1212 

.1208 

.0846 

• C41 C 

. 0G8 7 

-.0048 

C.CCCO 

-51U0FF 







-.6602 

-5 1 WL wF 







— . C758 

-51UL WE 







-.5048 

- 51 WL wW 







-.0688 

-51ULWW 

-.5879 

-.4451 

-.3475 

-.2578 

- • 2860 

- .2676 

-.2811 

-51WLFW 

.4811 

.4311 

.3435 

.2365 

. 1237 

• C 1 4b 

-.0620 

-51ULFw 

— . 36 76 

-.2399 

-.1504 

— • C87 5 

-.0376 

.0106 

.0600 

— 5 1 wDFW 

.2959 

.2647 

.2240 

.1879 

. 1641 

.1348 

.1588 

-51UUFW 

-.6778 

-.8915 

-1.0795 

-1.0384 

-.8222 

-.8*15 

-.8819 

-51WLFT 

-.2400 

-.1262 

.1221 

.444 7 

• 4443 

.2019 

-.0140 

-51ULFT 

-.5918 

—.6256 

-.5992 

-.4545 

-. 2621 

-.1215 

.0225 

- 5 IwOFT 

-.8614 

-.7557 

-.6181 

— . 5066 

-.5245 

-.5773 

-.5932 

-51UDWT 







-1.1188 

— 51wLwT 







.0084 

- 5 ILL wT 







-.4932 

-51WLTT 







-.2463 

-51CLTT 

-.0323 

- .0318 

-.0319 

-.0225 

- • C33 1 

- • C330 

-.0324 

-51WLFP 

-.0403 

-.0292 

-.0389 

- • 03^ 0 

- .0389 

-.0363 

-.0076 

-51ULFP 

-.0805 

-.0724 

-.0653 

- . 058 8 

-.0525 

-.0464 

-.0404 

- 5 1WL)F P 

.6737 

.6711 

.8689 

.8672 

• 6c6 1 

.8659 

. 8665 

-51U0FP 


-.0673 -51WLWP 
- . C394 -51ULWP 
-.0492 - 51»DhP 
1.1361 -51L0WP 
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ALPHA = 0 


2C.COQO 

30.0000 

40.0000 

5C.CC00 

60.0000 

70.0000 

80.0000 

90.0000 

CHi 

-.9887 

-.6382 

-.6555 

-.5093 

-.4270 

-.4061 

-.4108 

-.4065 

OlWLFF 

• . 3858 

.4837 

.49 33 

.4238 

.303 5 

.1636 

.0254 

-.0792 

C1ULFF 

1 -.7513 

- .5455 

-.3666 

-. 2304 

-.1336 

-.0604 

.0064 

.0792 

GUDFF 

.0529 

.1201 

. 1 2 C6 

.0847 

.0388 

• 0066 

-.0048 

C.COOO 

01UDFF 








-.6642 

ClftLfcF 








-.0917 

01ULWF 








-.5079 

ClHLWN 








-.0990 

OlULHW 

-. 752C 

-.5851 

-.4468 

-.3509 

-.3011 

-.2884 

- .2867 

-.2808 

ClfcLFh 

• <i 52 8 

• 4611 

.4150 

.3307 

.2255 

.1133 

.0040 

-.0727 

010LFM 

-.5171 

-.3507 

-.2259 

-.1376 

-.0749 

-.0249 

.0233 

.0727 

ClfcDFta 

.2763 

.2612 

.2545 

.2173 

.1836 

. 162 0 

.1544 

.1603 

01UDFW 

-.5261 

-.7369 

-1.0049 

-1.2490 

-1.1605 

-.8497 

-.8191 

-.e525 

OlWLFT 

-.3560 

-.3335 

-.21 87 

.0809 

.4898 

.4662 

.2106 

-. CC15 

C1ULFT 

-.5773 

-.6556 

-.7143 

-.7000 

-.5229 

-.2936 

-.1424 

.0015 

01KDFT 

-.5922 

-.9170 

-.7994 

-.6296 

-.4860 

-.5115 

- .5666 

-.5861 

01UDFT 








-1 . 0952 

OlwLhT 








.0132 

C1LLKT 








-.4489 

OltaLTT 








-.0934 

C1LLTT 

-.0330 

-.0319 

-.0314 

-.0315 

-.0320 

-.0326 

-.0325 

-.0328 

ClWLFP 

-.0417 

-.0397 

-.0367 

-.0383 

-.0384 

-.0363 

-.0357 

-.0070 

C1ULFP 

- • Co 95 

-.0799 

-.0715 

-.0646 

-.0584 

-.0523 

-.0462 

-.0404 

01WDFP 

.6761 

.87 33 

.8707 

.8683 

.8665 

• 6654 

.8651 

.8656 

01UDFP 


-.0662 OlHLWP 

- . C382 OHILwP 

-.CA92 OlHOhP 

1 • 134 5 C1U0HP 

ALPHA = 5 


20.0000 

30.0000 

40.0000 

50.0000 

60.CCOC 

70. 0000 

80 .OCOU 

9C.C000 

CHI 

-1.016? 

-.8400 

-.6562 

-.5093 

-.4570 

- . 4U6 0 

-.4108 

-.4064 

5 1 hL FF 

.3860 

.4849 

.4947 

.4248 

.304 0 

.1642 

• 025o 

-.0788 

51ULFF 

-.7536 

-.5471 

-.3675 

-.2308 

-.1343 

- .0606 

.0080 

.0788 

51*0FF 

.0532 

.1210 

.1217 

.0852 

.0411 

.0087 

-.0048 

C.COOO 

51UDFF 








6671 

51 WLtaF 








-. 1069 

51UL*F 








-.5117 

51kLhW 








-.1287 

5 1ULWW 

- .7442 

-.5840 

-.4496 

-.3552 

-.3051 

-.2914 

- .2905 

-.2815 

5 1 WL Fh 

* 32 0 

.4430 

.4005 

. 3194 

.2159 

. 1 04U 

-.0058 

-.0826 

5 1ULFW 

-.4972 

- .3359 

-.2135 

-.1259 

- .0633 

-.0131 

.0352 

• C64 7 

5 1 tuDF w 

.2574 

.2666 

.2441 

.2102 

.1790 

• 1 59u 

.1530 

.1606 

51U0FW 

-.5518 

-.8048 

-1.1464 

-1.4726 

-1.3066 

- • 6648 

-.8254 

-.8320 

51*LFT 

- .4506 

-.4467 

-.3341 

.0381 

.5656 

.4948 

.222c 

.0113 

51ULFT 

- .6235 

-.7301 

-.8227 

- .8280 

-.6030 

-. 32bo 

-.1652 

-.0198 

51WDFT 

-1 .0606 

-.98 72 

-.8567 

—.6434 

- .4570 

-.4979 

-.5610 

-.5807 

51U0FT 








-1.C729 

5 lrtLtnT 








.0181 

5 lOL *T 








-.4369 

5 1 bL TT 








.0409 

51ULTT 

-.0326 

-.0314 

-.0309 

-.0310 

-.0315 

-.0321 

-.0319 

-.0323 

5 lfcLFP 

-.0411 

-.03 91 

-.0381 

-.0377 

-.0378 

-.0376 

-.0351 

- . CC63 

51ULFP 

- .Ottdb 

-.0793 

-.0714 

-.0645 

-.0582 

-.0521 

-.0461 

— . C404 

5 1 HDFP 

.8758 

.8729 

.8701 

.8677 

.8658 

• 8646 

.8642 

• 6646 

51L0FP 


- . C652 SlHLKiP 

— . C369 51ULHP 

- . C493 5 1 HO HP 

1.1328 51U0HP 
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APPENDIX B 


ALPHA = 10 


2L. COOL 

30.0000 

40.0000 

50.0000 

60.0000 

-l .01*6 
.3871 
-.6830 
.0546 

-.8422 

.4872 

-.5499 

.1229 

-.6569 

.4971 

-.3652 

.1233 

-.5092 

.4266 

-.2321 

.0863 

-.4266 

.3053 

-.1355 

.0417 


-.7385 

-.5844 

-.4535 

-.3602 

-.3098 

.4129 

.4265 

.3877 

.3096 

• 2 C76 

-.4800 

-.3230 

-.2026 

-.1155 

- • C52 7 

.2369 

.2522 

.2336 

.2027 

.1736 

-.5771 

-.8825 

-1.3255 

-1.7749 

-1.4724 

-.5604 

-.5871 

-.4841 

-.0044 

.6840 

-.6750 

-.8173 

-.9576 

-.9946 

-.6948 

1.1445 

— 1 .0766 

-.9333 

-.6581 

-.4158 


-.0321 

-.0309 

-.0305 

-.0306 

- .C311 

- • 0404 

-.0385 

-.0374 

-.0371 

- . C3 72 

- .0881 

-.0787 

-.0709 

-.0641 

-.0579 

. 6 753 

.8723 

.8695 

.8670 

.8651 


ALPHA = 16 


20.00U0 

30.0000 

40.0000 

50.0000 

60.C00C 

-1.0249 

-.8454 

-.6581 

-.5090 

- • 426 C 

• 3 d96 

.4912 

.5011 

.4299 

.3080 

-.7639 

-.5547 

-.3726 

-.2348 

-.1375 

.0578 

.1265 

.1263 

.0883 

.0427 


-. 7344 

- .5869 

-.4596 

-.3672 

- .3161 

. 3*21 

.4091 

.3745 

.299 8 

.1994 

-.4629 

-.3104 

-.1518 

-.1049 

-.0416 

.2170 

.2350 

.2208 

.1931 

.1662 

-. 60 <tC 

-.9897 

- 1.6116 

- 2.3032 

- 1.6787 

-.7229 

-.8075 

-.7 369 

-.0439 

.9005 

-.7436 

-.9418 

- 1.1686 

- 1.2737 

-.8174 

1.2718 

- 1.2202 

- 1.0644 

-.6723 

-.3438 


- .0315 

-.03 04 

-.0299 

-.0300 

- • C305 

-.0396 

-.0377 

-.0367 

-.0363 

-.0364 

0873 

-.0781 

- • 07 C5 

-.0638 

-.0577 

. d74t> 

.8716 

.8686 

.8661 

.8640 


70 • 0000 

ec.occu 

9C.C000 

CHI 

-.4056 

- .41 05 

-.4063 

101WLFF 

.1654 

.0270 

-.0776 

1C1ULFF 

-.0620 

• CC68 

.0776 

IClwOFF 

.0091 

- • C C46 

c.cooo 

101UDFF 



-.6687 

1C1WLWF 



-.1213 

101ULWF 



-.5160 

lClwLWW 



-.1579 

101ULWW 

-. 2550 

-.2531 

-.2829 

lOlftlFh 

.0957 

-.0147 

-.0522 

1C1ULFW 

-.0023 

.0464 

.0562 

lOlwLFW 

. 1551 

.1506 

.1558 

10 1CDFW 

— • 5266 

-.8399 

-.8202 

1C1WLFT 

. 5301 

.2378 

.0247 

101ULFT 

-.3680 

-.1904 

-.0420 

1C1WDFT 

-.4837 

-.5545 

-.5769 

101UDFT 



-1.0523 

1 ClWLUT 



.0234 

1 0 IOLmT 



-.4547 

101WLTT 



.1767 

1Q10LTT 

— . 031 6 

-.0214 

-.0318 

1C1WLFP 

-.0370 

-.0344 

-.0057 

1C1ULFP 

-.0520 

-.0461 

-.0404 

1C1WDFP 

• 8638 

.8633 

.8636 

1C1UDFP 



-.0642 

10 1WL WP 



-.0358 

lOlULhP 



-.0495 

10UDWP 



1.1311 

lClUDwP 


70.0000 

tiO.GCUO 

90.0000 

CHI 

-.4050 

-.4100 

-.4060 

16ULFF 

.16 7 8 

.0296 

-.0751 

161ULFF 

-.0644 

.0044 

.0751 

161WDFF 

.0100 

-.0042 

C.COOO 

161UDFF 



-.6689 

leluLWF 



-.1374 

161ULWF 



-.5219 

16 1 WL Wirt 



-.1919 

lfclULtaW 

-.3002 

-.2971 

-.2858 

1 6 1 WL F W 

• C 87 3 

-.0241 

-. 1026 

1 61 ULFirt 

.0094 

• 0586 

.1089 

161WDFW 

.1492 

.1463 

.1572 

lfclUDFW 

-.9850 

-.8686 

-.8179 

161WLFT 

.5804 

.2609 

.0421 

161ULFT 

-.4212 

-.2246 

-.0707 

16 lwDFT 

-.4657 

-.5481 

-.5747 

161LDFT 



-1.0308 

1 6 1 WL wT 



.0306 

1 6 10L hT 



-.5208 

161WLTT 



.3715 

161ULTT 

-.0310 

—.0308 

- . C31 2 

161WLFP 

-.0362 

-.0337 

- • CC51 

161ULFP 

-.0519 

-.0461 

-.0405 

161WDFP 

. 8627 

.8621 

.8622 

161U0FP 



-.0630 

161WLWP 



-.0345 

161ULWP 



- . C498 

16 lWUWP 



1.1288 

161U0WP 
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APPENDIX B 


2-1 RECT. 
ALPHA = -10 


CO 00 

30.0000 

40.00CC 

50.0000 

60.C00C 

7G.C000 

8C.CCCC 

9C.C000 

CHI 

1143 

-.9164 

-.7110 

-.5449 

-.4485 

-.4196 

-.4207 

-.4160 

-102WLFF 

4331 

.5430 

.5535 

.4768 

.3431 

. lb7o 

.0330 

-.0862 

-102ULFF 

62 33 

-.5931 

-.3339 

-.2426 

-.1365 

-.0577 

.013o 

.0862 

— 10 2 rtDFF 

1121 

.1757 

.1666 

.1173 

• C59 8 

.0162 

-.0038 

C.CCCO 

-1 C2UUFF 








-.6718 

-10 2toL hF 








— • C650 

-1 02 GL *F 








-.5064 

-1 C2WL*W 








-.0430 

—1 C2UL Ink 

8468 

-.6433 

-.4761 

-.3620 

-.3030 

-.2872 

-.2677 

-.2820 

- 1 0 2 *L Fw 

56GC 

.5641 

.5039 

.4020 

.2795 

.1528 

• C306 

-.0574 

-102ULF* 

6131 

-.4137 

-.2679 

-.1678 

- .C993 

-. C4o9 

.0020 

.0515 

-1 02 UUFh 

3 b06 

.3654 

.3157 

. 2571 

• 2 C6 0 

.1709 

.1545 

.1561 

-1G2UDFW 

3264 

- .6925 

-.8845 

-1.0485 

-1.0325 

-.8689 

-.8874 

-.9808 

-1 02 hL FT 

2 36t 

- .1626 

-.0383 

.1779 

.4666 

.4794 

.2213 

-.0276 

-1C2ULFT 

5331 

-.5940 

-.6104 

-.5705 

-.4328 

-.2496 

-.1047 

.0504 

-102KUFT 

904 0 

-.8245 

-.7202 

-.5956 

-.5015 

-.5234 

-.5853 

-.6042 

-1G2UOFT 








-1.2103 

-1C2WLWT 








.0055 

-1C2UL*T 








-.6270 

- 1 0 2»L TT 








-.4956 

-102LLTT 

03o2 

-.0349 

-.0344 

-.0345 

- . C350 

-.0357 

- . 0356 

-.0363 

-1C2KLFP 

0465 

-.0443 

-.0431 

-.0426 

- . C42 7 

-.0425 

-.0399 

-.0093 

-1C20LFP 

0973 

— • 0666 

-.0777 

-.0699 

- .0629 

-. C561 

-.0495 

-.0431 

-1C2*DFP 

904 5 

.9003 

.8967 

.893 7 

.8915 

.6903 

• 8 9u4 

.8916 

-1C2UOFP 








-.0733 

-102wL*P 








-.0448 

— 1 G2UL l*P 








-.0524 

- 10 2wDkP 








1.175C 

-1020DWP 




ALPHA = 

-5 






CJOO 

3 C • 0000 

40.0000 

50.0000 

60 .0000 

7C.COOO 

60.0000 

9C.C000 

CHI 

1147 

-.91 71 

-.7112 

-.5450 

-.4467 

-.4284 

-.4210 

-.4161 

-52WLFF 

4333 

.5451 

.5529 

.4759 

.3419 

. 1863 

.0317 

-.0875 

-52ULFF 

8223 

-.59 22 

-.3926 

-.2413 

-.1352 

-.0564 

.0151 

.0875 

- 5 2 fcUFf* 

,1100 

.1745 

. 1658 

.1165 

.0591 

• Cl 5 o 

-.0041 

c.cccc 

-52U0FF 








-.6775 

- 5 2*L t%F 








-.0837 

-520L*F 








-.5096 

-52l*L 1 h)h 








-.0764 

- 52UL*»W 

8346 

-.6383 

-.4768 

-.3651 

-.3062 

2893 

-.2883 

-.2811 

- 5 2mLFw 

5332 

.5399 

.4842 

.3863 

. 2 66 5 

.1407 

.0185 

-.0693 

-52ULFW 

3847 

-.3926 

-.2507 

-.1522 

-.0842 

C319 

.0169 

.0663 

-52*OFh 

3565 

.3486 

.3042 

.2500 

.2020 

. 1693 

.1551 

.1588 

-52UDFW 

55 2 2 

- .7488 

-.9859 

-1. 1927 

-1.1399 

-.8856 

-.8685 

-.9314 

— 52 wL F T 

3140 

-.2710 

-.1437 

.1341 

.4955 

.4965 

.2277 

-.0150 

-52ULFT 

59 3 7 

-.6549 

-.6918 

-.6604 

- .496 6 

-.2602 

-.1259 

.0263 

— 5 2hDF T 

9553 

-.8737 

-.7572 

-.6053 

-.4835 

-.5085 

-.5738 

-.5940 

-52U0FT 








-1.1739 

- 52 WLnT 








• CC97 

-52ULWT 








-.5258 

-52*LTT 








-.2741 

-52LLTT 

0358 

- .0345 

-.0339 

-.0340 

-.0346 

-.0352 

-.0351 

-.0357 

- 5 2 fcLFP 

0459 

.0437 

-.0425 

-.0421 

-.0421 

-.0419 

-.0392 

-.0085 

-52ULFP 

0965 

-.0859 

-.0771 

-.0694 

-.0625 

-.0556 

-.0493 

-.0430 

-52WDFP 

9041 

.6999 

.8961 

.8930 

.8908 

.8696 

.8896 

.8906 

-52UDFP 








-.0721 

-52WLHP 








-.0433 

-52ULWP 








-.0523 

-52WUWP 








1.1733 

-52UDWP 
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appendix b 


ALPHA = 


2U.00CC 

30.0000 

40.0000 

50.0000 

60.0000 

70.0000 

60 .0 COO 

SC.COOG 

CHI 

-1.1164 

-.9185 

-.71 18 

-.5450 

- .4487 

42u0 

-.4211 

-.4162 

Q2WLFF 

. 4325 

.5424 

.5534 

.4760 

.3417 

.1659 

.0312 

-.0880 

C20LFF 

-.8231 

- .5926 

-.3925 

-.2409 

-.1347 

-.056 0 

.0156 

.0880 

C2MDFF 

.1091 

. 1744 

.1658 

.1164 

.0585 

.0156 

-.0041 

C.0000 

C2U0FF 







-.6823 

C2MLMF 








1016 

C2UL*F 








-.5137 

0 2WL Mto 








-.1094 

C 2UL MM 

- • d2 33 

-.6354 

-.4789 

-.3692 

-.3102 

-.2922 

-.2698 

-.2811 

02MLFM 

.5082 

.5177 

.4663 

. 3723 

.2546 

. 1296 

.0073 

-.0806 

C2LLFM 

-.5595 

-.3738 

-.2352 

-.1379 

-.0701 

-.0176 

.0312 

.0806 

02MDFM 

. 3 371 

.3321 

.2528 

.2426 

.1976 

.1670 

.1547 

.1603 

C2UDFM 

-.5797 

-.8142 

-1.1120 

-1.3813 

-1.2745 

-.9136 

-.8621 

-.8552 

02MLFT 

-.4021 

-.3748 

-.2462 

.0892 

.5473 

. 5215 

.2370 

-.0024 

C2ULFT 

- .6 346 

-.7262 

-.7906 

-.7730 

-.5727 

-.3148 

-.1491 

.0024 

C2mDFT 

-1.0178 

-.9352 

-.8052 

-.6174 

-.4596 

-.4536 

-.5642 

-.5861 

02UDFT 








-1.1432 

C2MLMT 








.0136 

02LLHT 








-.4725 

02 ML IT 








-.1045 

C2ULTT 

-.0353 

-.0340 

-.0335 

-.0336 

-.0341 

-.0346 

-.0345 

-.0351 

C2WLFP 

- . C45 3 

-.0431 

-.0419 

-.0414 

-.0414 

-.0412 

- .0385 

-.0078 

02ULFP 

-.0956 

-.0852 

-.0765 

-.0690 

-.0621 

- • 0556 

-.0451 

-.0429 

02MUFP 

.5036 

.8993 

.6955 

.8923 

.8900 

• 686 7 

• 6886 

.8896 

02UDFP 








-.0709 

02mLmP 








-.0419 

02UL MP 








-.0523 

C2M0MP 








1.1714 

C2UD*P 


ALPHA = 


20.0000 

30.0000 

40 • 00 CO 

50.0000 

60 .COOC 

7C.00UU 

60.0000 

9C.C000 

CHI 

-1.1191 

-.9204 

-.71 25 

-.5445 

-.4466 

-.4200 

-.4211 

-.4161 

52MLFF 

.4327 

.5438 

.5550 

.477 L 

. 3423 

.1863 

.0316 

-.0875 

52ULFF 

-.8256 

-.5943 

-.3935 

-.2414 

-.1352 

-.0564 

.0151 

.0875 

52MDFF 

.1055 

. 1754 

.1668 

.1170 

.0592 

.0156 

-.0041 

c.cooo 

52U0FF 








-.6860 

52MLMF 








-.1188 

52ULMF 








-.5187 

52MLMM 








-.1420 

52ULWM 

-.6147 

-.6343 

-.4823 

-.3742 

-.3150 

-.2559 

-.2920 

-.2820 

52MLFM 

.4651 

.4975 

.4503 

.3599 

.2444 

.1196 

- • C C31 

-.0913 

52ULFW 

-.5373 

-.3573 

-.2213 

-.1248 

-.0571 

-.0046 

.0446 

.0942 

52mDFM 

.3161 

.3160 

.2814 

.2349 

.1925 

. 1638 

.1533 

. 1608 

52U0FW 

-.6061 

-.8699 

-1.2700 

-1.6312 

-1.4376 

-.9516 

-.8677 

-.8714 

52MLFT 

-.5047 

-.5004 

-.3743 

.0426 

.6336 

. 5545 

.2495 

.0105 

52ULFT 

-.6914 

- . 8096 

-.91 20 

-.9163 

-.6622 

-.3540 

-.1747 

-.0218 

52MDFT 

-1.0936 

-1.0129 

-.8645 

-.6318 

-.4265 

-.4783 

-.5561 

-.5803 

52UOFT 








-1.1181 

52ML MT 








.0176 

52ULMT 








-.4583 

52MLTT 








.0439 

52ULTT 

-.0348 

-.0335 

-.0330 

-.0331 

-.0336 

-.0341 

-.0340 

-.0346 

52MLFP 

- .0446 

-.0424 

-.0412 

-.0408 

-.0408 

-.0405 

-.0376 

- • CC72 

52ULFP 

-.0946 

-.0846 

-.0760 

-.0686 

-.0618 

-.0554 

-.0490 

-.0429 

52MUFP 

.9031 

.8987 

.8948 

.8916 

.8892 

.8878 

.8876 

• 8884 

52LUFP 


-.0698 52*LwP 
-.0405 52ULWP 
-.0524 52WDMP 
1.1694 52UOKP 
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APPENDIX B 


ALPHA = 10 


20.0000 

30.0000 

40.0000 

50.0000 

60.0000 

70.0000 

80.CC00 

90.0000 

CHI 

- 1 • 1 23 C 

-.9229 

-.7134 

-.5449 

-.4462 

-.4196 

-.4208 

-.4160 

102MLFF 

.4340 

.5463 

. 5576 

.4791 

.343 8 

.1876 

.0330 

- • C862 

102OLFF 

-.8298 

-.5974 

-.3954 

-.2429 

-.1365 

-.0577 

.0138 

.0862 

102MDFF 

.1111 

.1775 

.1666 

.1182 

• C60 C 

.0162 

-.0038 

C.CCOO 

102UDFF 








-.6884 

1 02mL MF 








-.1351 

1 02UL MF 








-.5244 

1C2MLMM 








-.1740 

10 2ULMM 

- .6086 

-.6351 

-.4869 

-.3801 

-.3205 

-.3002 

-.2551 

-.2838 

1C2ML FM 

.4638 

.4753 

.4361 

.3491 

.2355 

• 1106 

-.0125 

-.1012 

102ULFW 

-.3182 

- .3430 

-.2092 

-.1133 

- • C453 

• 0076 

.0572 

.1071 

1 02 MDFM 

.2956 

.3002 

.2700 

.2269 

• 1866 

.1599 

.1509 

.1600 

1C2UDFM 

-.6364 

-.9770 

-1.47C6 

-1.9698 

-1.6234 

— • 9986 

-.8851 

-.8601 

102MLFT 

-.0268 

-.6566 

-.5412 

-.0037 

. 7676 

.5949 

.2657 

.0242 

1C2ULFT 

-.7493 

-.9074 

-1.0632 

-1.1031 

-. 7653 

-.3983 

-.2033 

-.0469 

1 02MDFT 

-1.1868 

-1.1125 

-.9536 

-.6476 

-• 380C 

- • 4625 

-.5495 

-.5766 

102UDFT 








-1.098b 

102WLMT 








.0221 

1 C2UL mT 








-.4795 

102MLTT 








.1941 

1C2ULTT 

-.0343 

-.0330 

-.0325 

-.0326 

- • C33 C 

- • C33o 

-.0334 

-.0340 

1 C2*LFP 

-.0439 

-.0417 

-.0405 

-.0401 

-.0401 

-.0398 

-.0371 

-.0065 

102ULFP 

-.0941 

-.0840 

-.0755 

-.0682 

-.0615 

- • 055 c 

- .0 450 

-.0429 

1C2M0FP 

.9024 

• 89 80 

.8940 

.8907 

.8882 

. 8868 

.8 865 

.8872 

102UCFP 








-.0687 

102MCMP 








-.0393 

1 C2UL mP 








-.0526 

1G2MDMP 








1.1673 

1C2UDMP 




ALPHA = 

16 






20. OOOO 

30. 00 CO 

40.0000 

50.0000 

60 • COO C 

7 C . CUGu 

80.0000 

5C.CCC0 

CHI 

-1.1289 

- .9266 

-.7146 

-.5447 

-.4475 

-.4189 

-.4203 

-.4156 

162MLFF 

.4368 

.5508 

.5621 

.4828 

.346 7 

. 1504 

.0356 

-.0834 

162ULFF 

-.8372 

-.6028 

-.3992 

-.2458 

-.139 2 

-.0604 

.0111 

. 0834 

162*DFF 

.1146 

.1815 

.1718 

.1205 

.0615 

• U 17 2 

-.0032 

C.COOO 

162U0FF 








-.6895 

1 6 2 ML MF 








-.1533 

1 62UL MF 








- • S320 

1 62ML MM 








-.2114 

162ULMM 

- . 8043 

-.6362 

-.4940 

-.3882 

-.3279 

-.3064 

-.2599 

2873 

162MLFM 

.4405 

• 46 00 

.4215 

.3384 

. 2267 

.1020 

- .0224 

-.1122 

1 62ULFM 

-.4992 

-.3290 

-.1971 

-.1013 

-.0326 

. 020 6 

.0 710 

.1215 

162mOFM 

.2716 

.2816 

.2563 

.2169 

.1792 

. 1540 

.1466 

.1576 

162UDFM 

—.6671 

-1 .0976 

-1.7918 

-2.5625 

-1.8555 

-1 . 0666 

-.9217 

-.8630 

1 62 mL F T 

-.8076 

-.9022 

-.8229 

-.0463 

1 • Oil 7 

.6520 

.2505 

.0420 

162ULFT 

- .8260 

-1.0475 

-1 .3002 

-1.4164 

-.5035 

- . 4589 

-.2426 

-.0792 

1 62 MDFT 

-1.3284 

-1.2722 

-1.0997 

-.6630 

-.2993 

-.4432 

-.5436 

-.5751 

162UUFT 








-1.C835 

162MLMT 








.0284 

1 6 2UL inT 








-.5578 

162MLTT 








.4104 

1620LTT 

— . o 33 6 

-.0324 

-.0319 

-.0319 

-.0324 

-. 0329 

-.0328 

-.0334 

16 2ML FP 


0429 

-.0408 

-.0397 

-.0392 

-.0392 

-.0390 

-.0363 

-.0058 

162LLFP 

0^32 

-.0833 

-.0750 

-.0679 

- .0613 

-.0551 

-.0490 

-.0430 

162MDPP 

9015 

.8970 

.8930 

.6896 

.8870 

.8855 

.8851 

. 8857 

1620UFP 


C675 162ML mP 
-.0379 162ULMP 
-.0530 162MDMP 
1.1647 16200 mP 
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1.5-1 


ALPHA - -10 


2 0.0000 

30.0000 

40.00CC 

50.0000 

60.0000 

70.0000 

80.0000 

90*0000 

-.9307 

-.8193 

-.6189 

-.4954 

-.4301 

-.4205 

-.4325 

-.4308 

.3435 

.42 76 

.4341 

.3721 

.2661 

. 1442 

.0248 

-• C635 

-. 7065 

- .5263 

-.3673 

-.2445 

- .1546 

C826 

-.0115 

.0635 

- . U 1 04 

.06 05 

.0724 

.0516 

.0231 

.0028 

-.0045 

0.0000 

-.6806 

-.0507 

-.5458 

-.0304 



-.7394 

- .5771 

-.4461 

-.3594 

-.3194 

-.3161 

- .3243 

-.3208 

.4381 

.4449 

.3989 

.3183 

.2212 

. 1212 

.0256 

- . C4C3 

-.5804 

- .3969 

-.2744 

-. 1870 

- .1233 

- « C699 

-.0161 

.0401 

.2299 

.2429 

.2244 

.1956 

. 1699 

.1533 

.1465 

.1483 

-.4928 

-.6343 

-.7694 

-.9028 

-.8638 

-.7591 

-.8069 

-.8868 

-.1834 

-.1264 

-.0114 

. L906 

. 3968 

. 3631 

.1518 

- • C35 7 

- .4 884 

-.5201 

-.5295 

-.4856 

-.3638 

-.2197 

-.0977 

• C438 

- .8742 

-.7985 

-.7024 

-.5945 

. 5231 

-.5395 

- .5747 

-.5824 

-1.1251 








-.0106 

-.5329 

-.2457 

-.0395 

-.0384 

-.0381 

-.0386 

-.0395 

-.0406 

-.0406 

-.0410 

-.0586 

- .0542 

-.0531 

- • C53 C 

— • C536 

-.0542 

- .0 524 

-.0264 

-.0971 

-.0858 

-.0762 

-.0678 

-.0600 

-.0525 

— .0450 

-.0376 

.846 5 

.8485 

.8497 

.8507 

.651 8 

.6533 

.8555 

.8586 
- . C 79 7 
- . C704 
-.0470 
1.1224 




ALPHA = 

-5 





2 0 .0000 

30.0000 

40.00CC 

50.0000 

60 • 000 C 

7C.CC00 

8C.OOOO 

SO.COOO 

-.931 0 

-.7771 

-.6192 

-.4956 

-.4304 

-.4209 

-.43^8 

-.4309 

.3422 

.4268 

.4336 

.3714 

.2651 

. 1432 

• C238 

- • C645 

- . 7077 

-.5255 

-.3664 

-.2436 

-.1538 

-.0817 

-•0105 

.0645 

- .01 1 9 

.0595 

.0717 

.0510 

• C22 6 

. 0025 

-.0047 

C.CCCO 
-.6844 
— . C63 8 
5471 
-.0562 

-.7272 

-.5721 

-.4455 

-.3608 

-.3210 

-.3170 

-.3242 

-.3198 

.4182 

.4265 

.3836 

. 3060 

.2106 

. 1110 

.0154 

- . C5C6 

-.5387 

- .3809 

-.261 5 

-.1755 

- .1124 

-.0592 

-.0056 

• C505 

.2140 

.2 306 

.2159 

. 1902 

• 1 66 6 

.1520 

.1470 

.1505 

52U3 

- .6888 

-.8604 

-1.0197 

-.9416 

-.7792 

- .8C86 

-.8652 

-.2463 

- . 1992 

-.0736 

.1686 

. 4220 

. 2761 

.1611 

- • C204 

- .52 36 

-.5758 

-.6008 

-.5596 

- • 4 1 4 C 

-.2469 

-.1175 

• C24 3 

-.^171 

- .83 81 

-.7295 

-.5982 

5C90 

-. 5291 

-.5662 

-.5757 

-1.1175 








-.0017 

-.4767 

-.1954 

-.0*91 

- .03b0 

-.0377 

-.0361 

- . C39 C 

-.0401 

-.0403 

- • C404 

- .0560 

-.0536 

-.0525 

-.0524 

- .C53C 

053 5 

-.0516 

- . C256 

-.0963 

- .0650 

-.0756 

-.0672 

- • C 596 

-.0521 

-.0447 

- • C375 

. d46 7 

.8985 

.8496 

.8504 

.8514 

.8529 

.6549 

. 8579 

_ 1 Q C 


-.0785 
— . C6 8 8 
— • C468 

1.1211 


CHI 

-103*LFF 
-103ULFF 
— 1 03toDFF 
- 1 C3U0FF 
-lC3*Lk<F 
-1C3ULWF 
-103*L*W 
— 1 0 3UL kfr 
-103WLFH 
-1C30LFH 
-1 C3wDFi* 
-1C30DFW 
-1C3HLFT 
-103ULFT 
-103WDFT 
-103UDFT 
-1C3NLWT 
- 103UL wT 
-1 03 hL TT 
-1C3ULTT 
-1C3WLFP 
-103ULFP 
-103WDFP 
-1C3UOFP 
- 10 3l*L hP 
-1030L*P 
-1C3HO *P 
-1C3UDWP 


CHI 

-53KLFF 

-53ULFF 

-53W0FF 

-53UUFF 

-53*L*F 

- 5 3 0L hF 

- 53WUnIn 
- 53ULWH 
-53*LFb 
-53LLFW 
-53*DF* 
- 53UDF W 
-53wLFT 
-53ULFT 

- 5 3*0 FT 
-530DFT 
-53*L*T 
-53ULWT 

- 53 WL TT 
-53ULTT 

- 5 3*L FP 
-53ULFP 
-5 3*DFP 
-530UFP 
-53 WO *P 

- 5 3UL *P 

- 53*0 UP 
-53UOWP 
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ALPHA = 0 


2C.CQ00 

30.0000 

40.00CC 

50.0000 

60.0000 

70.0000 

80 .Clou 

90.0000 

CHI 

- .9322 

-.7781 

-.6196 

-.4957 

-.4305 

-.4210 

-.4329 

-.4310 

U3MLFF 

. 241 5 

• 42 o9 

.43 39 

.3714 

.2645 

.1428 

.0234 

- • C64 8 

03 CL 8F 

-.7082 

- .5258 

-.3663 

-.2434 

-.1536 

-.0814 

-.0102 

.0648 

C3MDFF 

-.0127 

.0593 

.0716 

.0509 

.0224 

.0023 

-.0047 

c.ccco 

C30DFF 








-.6871 

03 ML MF 








-.0763 

C3ULMF 








-.5489 

C3MLMM 








-.0815 

0 3UL MM 

-.7173 

-.5686 

-.4460 

-.3628 

-.3231 

-.3184 

-.3248 

-.3194 

C3MLF* 

.3991 

.4095 

.3697 

. 2949 

.2011 

.1017 

.0057 

-.0604 

C3LLFM 

-.5193 

- .3666 

-.2499 

-.1650 

-.1022 

-.0491 

.0044 

. C6C4 

03 MOF In 

. I9d 1 

.21 82 

.2070 

. 1843 

.1630 

.1499 

.1464 

.1517 

C3UUF* 

-.5507 

-.7530 

-.9938 

-1.1693 

-1.0344 

-.8047 

-.8129 

-.8470 

03mLFT 

—.3228 

-.28 27 

-.1472 

.1507 

.4643 

.3936 

.1 723 

- . CC5 1 

C3ULFT 

-.5674 

-.6409 

-.6871 

-.6500 

-.4711 

-.2761 

-.1377 

• CC5 1 

0 3mDFT 

— .9 c9U 

-.8873 

-.7639 

-.6013 

- .4903 

-.5182 

-.5584 

-.5699 

C3LUFT 








-1.1049 

03mLMT 








. CC70 

03LLMT 








-.4481 

C3mLTT 








-.0737 

C3LLTT 

- .0387 

-.0376 

-.0373 

-.0377 

-.0386 

-.0395 

-.0297 

-.0399 

0 3 ML FP 

-.0553 

-.0529 

-.0518 

-.0517 

-.0523 

0527 

-.0509 

-.0248 

03LLFP 

-.0955 

- .0844 

-.0750 

-.0668 

-.0592 

-.C5lv 

-.0446 

-.0374 

C3 mOFP 

.8468 

.8484 

.8494 

.8501 

.8510 

.8523 

.8542 

.8571 

03UUFP 


-.C774 C3wL*P 

-.C673 OSULhP 

C467 03kDI*P 

1.1198 03UCWP 

ALPHA = 5 


20. COOL 

3 C • 0000 

40.0000 

50.0000 

60.0000 

70.0000 

cC .OCCO 

9 C • 0000 

CHI 

-.9341 

-.7743 

-.6201 

-.4957 

-.4305 

-.4206 

-.4329 

-.4309 

5 3 ML FF 

.3417 

.4279 

.4350 

. 3721 

.2653 

.1431 

.0237 

-.0645 

530LFF 

-.7100 

-.5271 

-.3670 

-.2438 

-.1539 

C81 6 

-.0105 

.0645 

5 3 MUFF 

-.0124 

. 06 00 

.0722 

.0512 

.0226 

.0025 

- .0047 

C.CCCO 

53UDFF 








-.6886 

5 3 ML MF 








-.ceei 

5 3LL mF 








-.5510 

53mL MW 








-.1064 

53ULMM 

-.7093 

-.5666 

-.4474 

— • 3656 

-.3258 

-.3204 

-.3259 

-.3198 

5 3h LFw 

. 3 al 6 

.3940 

.3572 

. 2850 

.1926 

.0933 

-.0032 

-.0697 

53LLFM 

-.5022 

-.3540 

-.2396 

-. 1555 

-.0929 

-.0398 

.0138 

.0698 

53mOFM 

.1823 

.2057 

.1980 

.1779 

.1586 

.1470 

.1450 

. 1518 

53UUFW 

-.5837 

-.8289 

-1.1365 

-1.3623 

-1.1392 

-.8342 

-.6201 

-.8314 

5 3 ML FT 

— • 4 l C C 

- .3838 

-.2373 

. 1400 

.5291 

.4164 

.1859 

* C 1 C4 

53ULFT 

-.6170 

-.7177 

-.7929 

-.7622 

- .5351 

-.3074 

- .1388 

- .0143 

53 MDF T 

-1.0324 

-.9492 

-.8081 

-.6024 

-.4650 

- • 5066 

-.5509 

-.5649 

53U0FT 








-1.0867 

53ML MT 








.0158 

53ULMT 








-.4406 

53MLTT 








.0360 

53ULTT 

—.0382 

- .0371 

-.0368 

-.0372 

-.C38C 

-.0390 

-.0392 

-.0393 

53MLFP 

- • 0 546 

-.0522 

-.0511 

-.0510 

-.0515 

-.0520 

-.0501 

-.0241 

53ULFP 

-.0947 

-.0837 

-.0745 

-.0663 

-.0589 

-.0516 

-.0444 

-.0374 

53M0FP 

• 846 7 

.8482 

.8491 

.8497 

.8505 

.8517 

.8534 

.8562 

53UDFP 


-.0763 53KLWP 
- . C658 53ULnP 
-.0468 53hDhP 
1.1182 53L0I.P 
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ALPHA = 10 


2U.C00G 

30 * 00 CO 

40.000C 

50.0000 

60.000C 

7 C . C 0 0 0 

60 . ouoo 

90.0000 

CHI 

-.9367 

-.7611 

-.6707 

-.4956 

-.4302 

-.4206 

-.4326 

-.4301 

103WLPF 

. 3426 

.4297 

.4368 

.3735 

.2664 

. 1441 

.0247 

-.0635 

103ULFP 

-.7130 

-.5253 

-.3685 

-.2449 

-.1549 

082b 

-.0115 

.0635 

1 C3 wDPF 

-.0113 

.0615 

.0735 

.0521 

• C23 1 

.0026 

-.0045 

0.0000 

103UDPF 








-.6889 

1C3*L*F 








-.0992 

103ULWP 








-.5535 

l C3 mL kki 








-.1307 

1 0 3UL tow 

-.7032 

-.5658 

-.4497 

-.3690 

-.3290 

-.3229 

-.3276 

-.3208 

IQ3*LP* 

• 36^5 

.3800 

.3462 

.2765 

.1852 

. C656 

-.0114 

-.0783 

1C3ULFW 

-.4874 

-.3432 

-.2306 

-.1471 

-.0845 

-.0313 

.0224 

.0786 

1 C3 WDFk 

• 1 066 

.1932 

.1886 

.1711 

.1534 

. 1432 

.1425 

.1508 

1C3UDFW 

-.6192 

-.9189 

-1.3197 

-1.6148 

-1.2510 

-.8664 

-.8298 

-.6183 

10 3 J*L FT 

-.3143 

-.5095 

-.3516 

.1441 

.6224 

.4439 

. 2022 

• C264 

1 03 OL FT 

-.6731 

- .80 89 

-.9252 

-.9037 

-.6053 

-.3409 

-.1810 

-.0345 

103wDFT 

-1.1105 

-1.0281 

- .8655 

-.5981 

-.4310 

-.^940 

-.5438 

-.5606 

103CDFT 








-1.0633 

1 C3ML kT 








.0249 

103ULUT 








-.4522 

1 03 hL TT 








.1467 

103ULTT 

-.0377 

-.0366 

-.0363 

-.0 36 7 

-.0375 

-.0385 

-.0386 

-.0387 

1 03 taLFP 

-.053a 

- .0514 

-.0504 

-.0502 

- . C 50 7 

-.0512 

-.0493 

-.0233 

1C3ULFP 

-.0939 

- .0831 

— . 0 74 C 

- • 0660 

-.0586 

-.0515 

-.0444 

-.0374 

103H0FP 

• d466 

. 6480 

.8487 

. 8492 

.8499 

.8509 

• d 526 

.8552 

1C3U0FP 


-.0752 1C3*LWP 
-.0644 103ULWP 
-.0469 10 3wD(nP 
1.1166 1 C3LDhP 

ALPHA = 16 


20.0000 

3U.00C0 

40.0000 

50.0000 

60 • C 000 

70. OuOu 

80.0000 

sc.coou 

OHI 

-.9406 

-.7836 

-.6215 

-.4954 

-.4297 

-.420^ 

-.4323 

-.4306 

1 63 IaLFF 

.3446 

.4329 

.4400 

. 3762 

.2686 

.1463 

.0269 

-.0614 

1630LFF 

-.7183 

-.5332 

-.3698 

-.2471 

-.1570 

— .0848 

-.0136 

.0614 

163WDFF 

- .0086 

.0644 

.0759 

.0538 

.0242 

.0035 

- • C C42 

C.0000 

163UDFF 








-.6875 

1 6 3mL taF 








-.1114 

163UL*F 








-.5565 

162taLfck 








-.1508 

1 63 UL Win 

- . 6v81 

-.5666 

-.4537 

-.3738 

-.3334 

- • 3 2o5 

- .3806 

-.3221 

1 63MLFW 

.34 78 

.3651 

.3348 

.2679 

. 1779 

.0761 

-.0201 

-.0879 

163ULFW 

-.4727 

-.3325 

-.2216 

-.1385 

-.0757 

—.0223 

.0317 

.0883 

163MOFM 

.1477 

.1781 

.1771 

.1622 

.1463 

.1374 

.1381 

.1480 

163UDFH 

- • 6t>40 

-1 .0508 

-1.61 08 

-2.0318 

-1.3796 

-. 5080 

-.8448 

-.8069 

1 63 nLFT 

-.6702 

-.7075 

-.5375 

. 1924 

. 7780 

.4835 

• 2260 

.0467 

163LLFT 

-.7503 

-.9426 

-1.1345 

-1.1289 

-.6548 

-.3846 

-.2100 

-.0604 

163WDFT 

-1.2297 

-1.1544 

-.9602 

-.5760 

-.3758 

-.4774 

-.5356 

-.5566 

163UDFT 








-1.C3C7 

1 63 hL toT 








.0367 

163UL*T 








-.4954 

163*LTT 








.2999 

163ULTT 

-.0371 

-.0360 

-.0357 

-.0360 

-.0369 

-.0376 

- .0379 

-.0381 

163WLFP 

-.0527 

-.0505 

-.0494 

-.0493 

-.0498 

-.0502 

-.0483 

-.0225 

163ULFP 

-.0931 

-.0824 

-.0735 

-.0656 

-.0563 

-.0513 

-.0443 

-.0375 

1 63hDFP 

.8464 

.84 76 

.8481 

.8485 

.8491 

.8500 

.8515 

• 8540 

1 6 3U0FP 


-.0740 163*LhP 

-.0628 163UL t*P 

-.0472 163W0HP 

1.1145 l fc300*iP 
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7X10 

ALPHA = .-10 


20. COOL 

30.0000 

40.0000 

50.CC00 

60.0000 

70.0000 

8C .0000 

90.0000 

CHI 

-1.1129 

-.9093 

-.7039 

-.5446 

— .458 C 

-.4396 

-.4519 

- . 44C3 

-1C4*LFF 

.4294 

.5293 

.5312 

.4527 

.3291 

.1947 

.0673 

-.0298 

-1C4ULFF 

- • a99 0 

-.6676 

-.4672 

-.3150 

-.207 1 

-.1252 

- • C491 

.0298 

-104*DFF 

.0860 

.1539 

.1489 

.1058 

• C56 7 

.0-209 

.0030 

C.CCCO 

-104LDFF 








-.6545 

-104WLHF 








-.0208 

-104ULHF 








-.5698 

-lC4kLkh 








-.0069 

-104ULM* 

-.9765 

- .7615 

-.5768 

-.4484 

- .3844 

-.3738 

-.3851 

-.3871 

-1C4WLF* 

.5157 

.5522 

.51 04 

.4163 

.2984 

.1784 

.0661 

-.0173 

-1C4ULFN 

-.7923 

- .5698 

-.3963 

-.2718 

-.1841 

-.1158 

-.0508 

.0165 

-104HDF* 

.2454 

.2647 

.2335 

.1836 

. 1372 

. 1050 

. C 8 82 

.0840 

-104UDFW 

-.9957 

-1.1313 

-1.1673 

-1.0441 

-.8348 

-. 7314 

-.7523 

-.7852 

-1C4*LFT 

-.1746 

.0105 

.25 13 

.4525 

.4697 

.3104 

.1111 

-.0534 

-1040LFT 

- • 8538 

-.8389 

-.7383 

-.5544 

-.351 7 

1976 

-.0729 

.0585 

-104*UFT 

-.7428 

-.5701 

- .40 62 

-.3109 

-.2118 

-.3602 

-.3511 

-. 2950 

— 1 04UOFT 








-1.C301 

-104WLWT 








-.0422 

-104ULKT 








-.4987 

-1Q41*LTT 








-.1713 

-1C4LLTT 

-.0688 

- .0709 

-.0761 

-.0828 

- .091 C 

-.1001 

-.1081 

-.1143 

-1 C4hLFP 

- . 2 a4 1 

- .2795 

-.2808 

-.2876 

-.2991 

-.3132 

-.3243 

-.3058 

-1C4ULFP 

- .2217 

-.1966 

-.1748 

-.1551 

-.1366 

-.1182 

- • C992 

-.0794 

— 10 4*DFP 

.6975 

.74 54 

.7873 

.8261 

.8645 

.9050 

.9495 

.9591 

— 1 04UDFP 








-.1811 

-104*L*P 








-.5131 

-1C40LHP 








-. 1314 

-1 C4*DkP 








1.2980 

-1C4UU*P 


ALPHA = -5 


2 0. OOUO 

30.0000 

40.0000 

50. 0000 

60.0000 

70.000U 

eo.occc 

9C.COOO 

CHI 

-1.1126 

- .9317 

-.7038 

-.5445 

-.4581 

-.4397 

- .4995 

-.4552 

-54»LFF 

.4286 

.5287 

.5307 

.4521 

.3264 

.1941 

.0666 

- • C3C 5 

-54LLFF 

-.8961 

- .6668 

— .46 64 

-.3143 

-.2065 

-.1246 

-.0484 

.0305 

-54WDFF 

.0850 

.1532 

.1483 

.1053 

.0563 

.0207 

.0029 

0.0000 

-54UUFF 








-.6564 

-54*LfcF 








-.0293 

- 54UL hF 








-.5708 

-54WL WW 








-.0228 

-54ULlnl* 

-.9263 

-.7542 

-.5742 

-.4482 

-.3849 

-.3742 

-.3851 

-.3667 

-54HLF* 

• 5008 

. 5380 

.4987 

.4073 

.2913 

.1720 

.0597 

-.0236 

— 54ULFh 

-.7731 

-.5564 

-.3864 

-.2636 

1767 

-. 1087 

-.0440 

.0232 

- 54WDFH 

.2333 

.2548 

. 2265 

.1793 

.1347 

.1039 

.0881 

.0851 

-54 LOFb 

-l. J703 

-1.2331 

-1.2740 

-1.1166 

-.8616 

-.7377 

-.7503 

-.7749 

-54HLFT 

-.2289 

-.0202 

.2610 

.4936 

.5052 

.3308 

.1296 

-.0329 

-54ULFT 

- *5365 

-.9331 

-.8257 

-.6163 

-.3892 

-.2244 

- .0960 

.0354 

-54HCFT 

7 7 6o 

-.5854 

-.3994 

-.2895 

-.2959 

-.3521 

- .3868 

-.3940 

-54UDFT 








-1.0236 

- 54i*L *T 








-.0253 

-54 UL fcT 








-.4769 

-54*LTT 








-.0997 

-540LTT 

-.0668 

-.0708 

-.0760 

-.0825 

-.0906 

-. 0996 

-.1075 

-.1136 

-54WLFP 

-.2825 

- .2778 

-.2791 

-.2835 

-.2971 

-.3110 

- .3218 

-.3031 

-54ULFP 

-.2202 

-.1952 

-.1 735 

-.1539 

-.1355 

-.1172 

-.0983 

-.0786 

-54*DFP 

.6995 

. 7471 

.7887 

.8272 

.8653 

.9034 

.9496 

.9987 

-54UUFP 


-.1799 -54*LHP 
-.5084 -54ULhP 
-.1299 -54*D*P 
1.2S74 -54U0fcP 
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APPENDIX B 


ALPHA = 0 


20.C00C 

3C.OOCO 

40.0000 

50.0000 

60 • GOO C 

70.0000 

80.000C 

5C.CCCC 

CHI 

-1.1130 

-.9045 

-.7039 

-.5445 

-.4581 

-.4396 

-.4521 

-.4552 

04wLFF 

.4 263 

.5287 

.5307 

.4520 

.3282 

.1939 

• 0664 

-.0307 

C4ULFF 

- • 846 1 

-.6668 

-.4663 

-.3141 

-.2062 

-.1244 

-.0483 

.0307 

04W0FF 

.0647 

.1531 

.1482 

.1052 

• C56 2 

• 0206 

.0028 

C.COOO 

04UDFF 







-.6577 

04HL VtF 








-.0374 

Q4CL hF 








-.5721 

C4WLbl* 








-.0383 

C4ULNW 

-.9501 

-.7484 

-.5727 

-.4487 

-.3859 

-.3750 

-.3656 

-.3868 

C4kL FW 

.4865 

.5246 

.4880 

.3993 

.2849 

. 1662 

.0537 

-.0296 

04ULFM 

-.7569 

-.6445 

-.3775 

-.2562 

- .1695 

-.1022 

-.0376 

• C296 

04*DFW 

.221 l 

.2449 

.2195 

. 1747 

.1319 

.1023 

• 0b76 

.0855 

04UDF* 

-l.lblO 

-1.3568 

-1.4012 

-1.1980 

-.8902 

-.7461 

- . 75 C C 

-. 7660 

C4MLFT 

- .2906 

-.0524 

. 2790 

.5463 

.5452 

. 3522 

.1484 

-.0127 

04ULFT 

- 1 .032 7 

-1 .0427 

-.9255 

-.6835 

-.4280 

-.2516 

- .1190 

• C 12 7 

G4W0FT 

-.8195 

-.6011 

-.3853 

-.2621 

-.2773 

-.3427 

-.3814 

-.3919 

04UDFT 







-1.C158 

C4*L*T 








-.0087 

C40LKT 








-.4648 

04*LTT 








-.0330 

C4ULTT 

- • 0667 

- . U7 C7 

-.0757 

- .0822 

-.0902 

-. 0991 

- • 1 C6 6 

-.1128 

04*LFP 

-.2607 

-.2760 

-.2772 

-.2837 

- .2949 

308b 

-.3191 

- . 3 C C3 

C4ULFP 

-.2167 

- . 1938 

-.1722 

-. 1527 

-.1344 

-.1102 

-.0974 

-.0779 

04hDFP 

. 701b 

.7487 

.7899 

.8281 

.8659 

.9056 

.9493 

.5981 

040DFP 


-•1787 C4*L*P 
-.5036 C4ULKP 
-.1285 C4*DWP 
1.2963 C4UD*P 


ALPHA - 5 


20. OOOO 

30.0000 

40.0000 

50.0000 

60.0000 

70.0000 

6 0 . C C o o 

9 C . 0000 

CHI 

-1.1141 

-.9103 

-.7042 

-.5445 

-.4581 

-.4396 

- .4521 

-.4552 

54kL FF 

• 42 o4 

.52 92 

.5313 

.4524 

.3285 

.1941 

• 0666 

- • C3 C5 

54ULFF 

-. 8992 

- .6675 

-.4667 

3144 

-.2065 

-.1247 

-.0465 

.0305 

54HDFF 

.0649 

.1535 

.1486 

.1054 

.0563 

.0207 

.0029 

C.CCCO 

54U0FF 








-.6584 

54*LnF 








-.0449 

54 UL hF 








-.5736 

54WL 








-.0534 

5 4UL hW 

- .939 1 

-. 7439 

-.5720 

-.4498 

-.3873 

-.3761 

-.3864 

-.3875 

54*L Fw 

. h 72 9 

.5122 

.4782 

.3921 

.2793 

• 1610 

• C 4 6 S 

-.0352 

54ULFW 

-.79 08 

- .5341 

-.3698 

-.2497 

- .1638 

-.0962 

-.0316 

.0356 

54*DFW 

.2090 

. 2351 

.2125 

.1700 

.1288 

.1002 

• C 864 

.0854 

54UDF* 

-L.2714 

-l . 5072 

-1.5515 

-1.2865 

-.919 5 

7563 

-.7513 

-.7584 

54*LFT 

- . 3b 1 1 

- .0653 

.3092 

.6133 

.5900 

.374o 

.1678 

• C C74 

54ULFT 

-1.1457 

-1.1716 

-1.0398 

-.7557 

-.4680 

-.2794 

-.1422 

-.0099 

54bDFT 

- . 6664 

— * 6l 64 

-.3639 

-.2273 

-.2556 

-.3317 

- .2 749 

-.3889 

54UDFT 








-1.C068 

54kLKT 








. C078 

540L fcT 








-.4615 

54taL TT 








.0316 

54LLTT 

- . Obbb 

- .0705 

-.0754 

-.0819 

-.0897 

0965 

-.1061 

-.1120 

54 hL F P 

- • c 7ob 

-. 27 39 

-.2750 

-.2815 

-.2925 

-.3060 

-.3164 

-.2975 

54ULFP 

-.2172 

-. 1925 

- .1710 

-.1516 

-. 1334 

-.1153 

-.0967 

-.0773 

54»iDF P 

. 70.5 4 

. 7501 

.7905 

.8288 

.8662 

.9056 

.9489 

.9971 

54UDFP 


-.1773 54wLhP 
-.4987 54UL bP 
-.1272 54hDWP 
1.2947 54LDWP 
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APPENDIX B 


ALPHA = 10 


2O.C0U0 30.0000 40.0000 50.0000 60.CC0C 7C.0000 80.0CCC SC.CCCO CHI 


-1.1158 -.9114 -.7047 -.5446 -.4580 -.4397 

.4291 .5304 .5325 .4533 .3292 .1947 

-.9012 -.6690 -.4677 -.3151 -.2072 -.1253 

.0857 .1545 .1494 .1060 .0567 .0209 


-.5312 -.7408 -.5724 -.4515 -.3892 -.3777 

.4600 .5008 .4695 .3859 .2746 .1563 

-.7277 -.5252 -.3633 -.2441 -.1584 -.0909 

.1426 .2254 .2055 .1651 .1253 .0977 

-1.4067 -1.6910 -1.7274 -1.3794 -.5475 -.7680 

-.4419 -.1170 .3572 .6971 .6392 .3989 

-1.2797 -1.3244 -1.1704 -.8322 -.5085 -.3061 

-.9203 -.6300 -.3321 -.1839 -.2306 -.3192 


-.0664 -.0702 -.0751 -.0814 -.0892 -.0978 

-.2764 -.2716 -.2727 -.2791 -.2695 -.3033 

-.2158 -.1912 -.1658 -.1505 -.1324 -.1145 

.7051 .7514 .7918 .8292 .8663 .5033 


- . 4 5 2U 

- • 4 63 C 

• 0b73 

- • 0298 

-.0501 

• 02*8 

.0 C30 

c.cooo 


-.6584 


-.0520 


-.5754 


-.0679 

- .2876 

-.3885 

.0435 

-.0404 

- .0 26 2 

• C4 1 2 

. C £47 

. C845 

-.7540 

-.7524 

.I860 

• C277 

- • 1 659 

-.0227 

-.3673 

-.3849 


-.9967 


.0246 


-.4668 


.0964 

-.1054 

-.1111 

-.2136 

-.294o 

-.0*60 

-.0767 

.*461 

.9959 


-.1759 


-.4939 


-.1261 


1.2927 


104*.LFF 
1C4ULFF 
1C4HCFF 
1C4UDFF 
1 C4WL n*F 
1040L8F 
1 C4WLWW 
104ULWW 
1 04 hL FW 
104ULFW 
104h0Fh 
1 C4UDFW 
1C4WLFT 
104LLFT 
1 C4 IkOF T 
104UUFT 
1C48L«T 
104ULWT 
104t»LTT 
1040LTT 
1C4ULFP 
1040LFP 
1048DFP 
1C4U0FP 
1 C4WLHP 
l CALL l*P 
1 C4mOWP 
104UDHP 


ALPHA = 16 


20.0000 30.0000 40.0000 

— 1 . 1 1 66 -.9132 -.7055 

.4306 .5325 .5345 

-.9048 -.6717 -.4697 

.0875 .1564 .1511 


-.9231 -.7388 -.5739 

.4453 .4884 .4604 

-.7145 -.5166 -.3569 

.1626 .2139 .1971 

-1.6116 -1.9668 -1.9741 

-.5351 -.i486 .4496 

-1.4766 -1.5474 -1.3506 

-.9957 -.6408 -.2735 


0661 

-.0698 

-.0746 

2734 

-.2687 

-.2697 

2141 

-.1897 

-.1685 

7069 

.7526 

.7925 


50.0000 60.0000 70.0000 

-.5447 -.4576 -.4394 

.4551 .3307 .1561 

-.3167 -.2083 -.1266 

.1072 .0575 .0214 


-.4543 -.3919 -.3601 

.3797 .2699 .1520 

-.2385 -.1529 -.0653 

.1591 *1208 .0941 

-1.4897 -.5786 -.7832 

.8228 .7031 .4300 

-.9274 - . 5585 -.3440 

-.1187 -.1972 -.3021 


-.0808 -.0885 -.0570 

-.2760 -.2867 -.3000 

-.1494 -.1314 -.1136 

.8295 .8661 . 904o 


80 .0000 

90.0000 

OHI 

-.4516 

-.4550 

1 64wLFF 

.0686 

-.0285 

164ULFF 

-.0504 

.0285 

lo4*0FF 

.0032 

C.COOO 

164UDFF 


-.6574 

l64*LfcF 


-.0596 

1 64UL hF 


-.5777 

!64*Lkk 


-.0845 

164ULHM 

-.3899 

-.3905 

164*LFh 

.0384 

-.0460 

164ULF* 

- .0204 

.0472 

164*0Flrt 

.0819 

. 0827 

164UUFW 

-.7590 

-.7472 

164hLFT 

• 2 128 

.0527 

1640LFT 

-.1954 

-.0609 

164WDFT 

-.3567 

-.2789 

1640DFT 


-.9835 

1 6 4imL WT 


.0435 

1 64UL hT 


-.4848 

1 64I&LTT 


.1778 

1640LTT 

-.1044 

-.1101 

1 64 fcLFP 

-.3102 

-.2912 

164ULFP 

-.0952 

-.0762 

164W0FP 

.9469 

.9942 

164UDFP 


-.1742 

164*LbP 


-.4881 

164UL*P 


-.1251 

164W0WP 


1.2898 

1 6400 fcP 
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TABLE I.- AIRFOIL ORDINATES 



NACA 16-015 

Tip sta.: c = 45.72 cm (18.00 in.) 




Modified NACA 16-017 



L.E. radius: 0.503 cm (0.198 in.) 


L.E. radius: 0.917 cm (0.361 


















TABLE II.- FUSELAGE ORDINATES 


X 


Tr 


X 

R 

cm 

in. 

cm 

in. 

0 

0 

0 



0 

1.70 

.67 

2.18 



.86 

3.38 

1.33 

3.05 



1.20 

6.78 

2.67 

4.24 



1.67 

10.16 

4.00 

5.13 



2.02 

13.54 

5.33 

5.87 



2.31 

20.32 

8.00 

7.01 



2.76 

27.10 

10.67 

7.90 



3.11 

40.64 

16.00 

9.17 



3.61 

54.18 

21.33 

9.40 



3.70 

67.74 

26.67 

10.16 



4.00 
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- Dimensioned drawing of model 





model mounted in Ames 2.13- by 3.05-m (7- by 10-ft) wind tunnel 






' 




(44- by 88-in.) rectangular test section 






by 6- ft) rectangular test section 





Figure 4.- Relative size of test sections with respect to model. Langley full-scale tunnel 
has been omitted because of its vastly greater size. (See fig. 5 for relative size of 
Ames full-scale tunnel and fig. 6 for relative size of Langley full-scale tunnel.) 





A- 72-1215 


(a) Overall view. 


Figure 5.- Fan-in-wing model mounted in Ames 12.2- by 24.4-m (40- bv 80- ft) tunnel. 



(b) View from "below" model. 
Figure 5.- Continued. 


A- 72- 1216 



A- 72-1218 


(c) View from "above” model. 


Figure 5.- Concluded. 


L-72-2844 

(a) Near view. 

Figure 6.- Fan-in-wing model mounted in 9.14- by 18.3-m (30- by 60- ft) 

Langley full-scale tunnel. 
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(a) Front view. 

Figure 7.- Fan-in- wing model mounted in 2.13- by 3.05-m (7- by 10-ft) insert in Langley full-scale tunnel. 
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Figure 7.- Continued. 



re 7.- Concluded. 
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Figure 8.- Uncorrected data for the model with the fans covered in several different test sections. The drag of the 
model at zero angle of attack has been removed from the data. The solid curve is a least-squares quartic faired 
through the data from the 9.14- by 18.3-m (30- by 60-ft) tunnel. 
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Figure 8.- Continued. 




Figure 8.- Concluded. 

















Figure 10.- Corrected data for the model with the fans covered in several different test sections. The drag of the 
model at zero angle of attack has been removed from the data. The solid curve is a least- squares quartic faired 
through the data from the 9.14- by 18.3-m (30- by 60-ft) tunnel. 
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Figure 10.- Continued. 
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v/v, 

(a) a - -5°. 

Figure 11.- Uncorrected values of lift coefficient as a function of v/Vj. The solid curve represents the sum of the 
vertical component of the static thrust and the lift of the wing expressed in coefficient form. It is assumed that 
there is no interference between the wing and the fans. 




Figure 11.- Continued. 








Figure 11.- Continued. 
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Figure 11.- Concluded. 
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C_) 


A/A 


(a) a = -5°. 

Figure 12.- Uncorrected values of drag coefficient as a function of v/Vj. The solid curve represents the sum of the 
momentum-theory value of fan drag and the drag of the wing expressed in coefficient form. It is assumed that 
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Figure 12.- Continued. 





Figure 12.- Continued. 





Figure 12.- Continued. 





Figure 12.- Concluded. 




Figure 13.- Uncorrected values of tail normal-force coefficient as a function of V/Vi. 
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Figure 13.- Continued. 
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Figure 13.- Continued. 
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Figure 13.- Concluded. 
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Figure 14.- Uncorrected values of ratios of lift to static thrust as a function of V/Vj . The curve is the sum of the 
vertical component of static thrust and the lift of the wing. It is assumed that there is no interference between 
the fans and the wing. 




Figure 14.- Continued. 
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Figure 14.- Continued. 




14.- Continued. 





Figure 14.- Concluded. 





Figure 15.- Continued. 




Figure 15.- Continue' 




Figure 15.- Concluded. 





s! 

S 


CO 

-J 

LU 

z 

z 

ZD 


cr 

»— 

u 

LU 

az 


LU 

d 


a 

& 


r: — n - e 
i z i z I ' 

SIT 5jT St 

“CD “CD • 

CVI CO CVJ CD _ 
XX xx XX 
CVJ 3“ CM :T 0J3* 


* CO 


CD 

z 

a: 

i— 

u 

ck 

SY 

pi2 

SfS 


— 1 —3“ CVJ w 

• w • • « 

— — CM 


d 6c 


CE 

—I 

a_ 


cnu, 

CD O 
tO 

XX 

:r o 
— cn 


E-' 
I t— 
3“U. 

a 

,_CD 

Xx 

CM O 
* 3“ 
CM — 


S 


(T) — 


o □ o < A a 



121 


Figure 16.- Uncorrected values of external drag-lift ratios as a function of Vj Vj. The curve is calculated from the 
momentum-theory values of the fan forces and the lift and drag of the wing. It is assumed that there is no inter- 
ference between the fans and the wing. 
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Figure 16.- Continued. 




Figure 16 .- Continued. 




Figure 16.- Continued. 
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Figure 17.- Uncorrected values of equivalent lift-drag ratios as a function of V/ Vj. The curve is calculated from 
the momentum-theory values of the fan lift, drag, and shaft power together with the lift and drag of the wing. It 
is assumed that there is no interference between the fans and the wing. 
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Figure 17.- Continued. 




Figure 17.- Continued. 
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Figure 17.- Continued. 





Figure 17.- Concluded. 
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100 1000 10000 100 1000 10000 
DISK LORDING. N/SQ M DISK LORDING. N/SQ N 

Figure 18.- Relationship between the present model and test sections and the model sizes determined by reference 7 

to have negligible wall effects. 
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Figure 19.- Influence of the ratio of fan area to wing area on "fan-induced" lift (from ref. 19). The solid symbols 
show the values measured during the present tests in the 1.12- by 2.24-m (44- by 88-in.) flat-oval test section 
a = 0°; V/V; = 0.4. 
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V, KNOTS 


Figure 20.- Comparison between flight test and uncorrected wind-tunnel test (Ames 
test 388) measurements of the elevator deflection required to trim the YOV-10 air- 
craft when fitted with a rotating- cylinder flap (ref. 20). Flap set at 60°/30°; c.g. at 
0.219c. 
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Figure 22.- Location of the YOV-10 aircraft with relation to the model sizes determined by reference 7 to have neg- 
ligible wall effects. 






Figure 23.- Wall-induced angle of attack Aa in the uncorrected data from the fan-in- wing model as a 

function of v/Vi in several different wind tunnels. 



Figure 23.- Continued. 







Figure 23.- Continued. 





Figure 23.- Continued. 







Figure 23.- Concluded. 
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function of V/Vi in several different wind tunnels. 






UJ 

UJ 

UJ 


□L 

— 1 

-J 

— 1 


> 

LD 

LD 

LD 


O 

Z 

Z 

Z 



cc 

tr 

CC 


l— 

1— 

t— 

L— 


cn 

u 

u 

u 


_j 

UJ 

UJ 

UJ 


Ll_ 

cn 

tn 

cn 


c — 

z — 

Z 

z — 

CO 

i z. 

l z 

1 * — ' 

1 

1 

:r — - 

rr — • 

cn »— 

10 U. 

UJ 

OJ l 

CU 1 

CD Lt. 

O l 

2 

' CD 

’ CD 

* 1 

• o 


OJ CD 

CU CD 

— LD 

cn ^ 

Z 

X X 

X K 

X X. 

x x 

m? 

c\j rr 

oj rr 

OJ IT 

m r- 

j — 

— • IT 

— • ZT 

OJ w 

*— « w 


1 w-» 

' w 

% 

* 


— 1 



CVJ 


O □ o <3 


O 0 

O 


$ 

O o 


-3 A 


Figure 24.- Continued. 









Figure 24.- Continued. 










Figure 25.- Wall-induced tail incidence Ai-p in the uncorrected data from the fan-in-wing model as 

a function of V/Vi in several different wind tunnels. 




Figure 25.- Continued. 






Figure 25.- Continued. 






Figure 25.- Continued. 




Figure 25.- Concluded. 
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Figure 26.- Wall-induced values of Q T /q c in the uncorrected data from the fan-in-wing model as 

a function of V/V; in several different wind tunnels. 





Figure 26.- Continued. 






Figure 26.- Continued. 
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Figure 26.- Continued. 



Figure 26.- Concluded. 
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Figure 27.- Wall-induced fan incidence Aip in the uncorrected data from the fan-in-wing model as 

a function of v/v, in several different wind tunnels. 













Figure 27.- Concluded. 
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Figure 28.- Wall-induced values of qp/q c in the uncorrected data from the fan-in- wing model as 

a function of V/Vi in several different wind tunnels. 
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Figure 28.- Continued. 
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Figure 29.- Corrected values of lift coefficient as a function of v/Vj. The solid curve 
represents the sum of the vertical component of the static thrust and the lift of the wing 
expressed in coefficient form. It is assumed that there is no interference between the 
wing and the fans. 
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Figure 29.- Continued. 
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Figure 29.- Continued. 





Figure 29.- Continued. 
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Figure 29.- Concluded. 
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Figure 30.- Continued. 





Figure 30.- Continued. 






Figure 30.- Concluded. 
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Figure 31.- Corrected values of tail normal-force coefficient as a function of V/Vj. 
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Figure 31.- Continued. 






UJ 

LU 

UJ 


cr 

— 1 

-J 

_J 

cr 

:> 

CD 

CD 

CD 

:> 

o 

z 

z 

z 

CD 


cr 

a: 

a: 


y— 

h— 

V— 

t— 

h— 

cr 

<_) 

CJ 

CJ 

cr 

_j 

UJ 

UJ 

UJ 

_j 

u_ 

or 

ct 

CK 

Ll 



Z — 

z — 

z 

Z — 

Z~ 

CQ 

l Z 

l z 

1 — 

i y— 

1 f— * 

_J 

zr — 

zr — 

OOh- 

LO u_ 

CO U- 

UJ 

CM 1 

CM 1 

CO u. 

O 1 

CD O 

• CO 

* 00 

* 1 

• o 

«C— . 

C\J CO 

CM CD 

—•to 

CO — 

— CO 

z 

X x 

x X 

XX 

XX 

XX 

ZD 

CMrr 

cm zr 

cm zr 

co r- 

z r cd 

h— 

— • zr 

-« =r 

CM — 

— • • — 

— co 








— 


•— < 

CM 

CD 


o □ O < J 


Figure 31.- Continued. 
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Figure 31.- Concluded. 
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Figure 32.- Corrected values of ratios of lift to static thrust as a function of V/Vj . The curve is the sum of the 
vertical component of static thrust and the lift of the wing. It is assumed that there is no interference between 
the fans and the wing. 
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Figure 32.- Continued. 
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Figure 32.- Continued. 
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Figure 33.- Concluded. 
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Figure 34.- Continued. 
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Figure 35.- Continued. 
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Figure 40 .- Continued. 
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Figure 40.- Continued. 
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